Summary of the invention
For defect of the prior art, the invention solves the problem of Reservoir Fracture density calibration in well logging.
In order to solve above technical matters, a kind of Reservoir Fracture well logging scaling method, specifically comprises:
Receive well-log information and seismic data in work area, obtain well logging FRACTURE CHARACTERISTICS parameter, according to described forecast for seismic data fracture density according to described interpretation of logging data;
Described well logging FRACTURE CHARACTERISTICS parameter and described earthquake prediction fracture density are carried out crossing and obtains two dimension and to cross figure;
Cross in figure in described two dimension, by well logging FRACTURE CHARACTERISTICS parameter level interval division earthquake prediction fracture density grade interval, according to described well logging crack grade interval, described earthquake prediction fracture density is demarcated.
Preferably, described well logging FRACTURE CHARACTERISTICS parameter comprises: one or more in depth side direction resistivity, depth side direction resistivity ratio, factor of porosity, fracture porosity.
Preferably, described according to described interpretation of logging data acquisition well logging FRACTURE CHARACTERISTICS parameter, step also comprises: carry out analysis interpretation according to micro resistor FMI in described well-log information, obtains fracture porosity.
Preferably, described described well logging FRACTURE CHARACTERISTICS parameter and described earthquake prediction fracture density are carried out crossing obtain two dimension and to cross figure, step also comprises: carrying out crossing with described earthquake prediction fracture density respectively according to described depth side direction resistivity, depth side direction resistivity ratio, factor of porosity, fracture porosity obtains multiple two dimension and to cross figure.
Preferably, describedly cross in figure in described two dimension, interval by well logging FRACTURE CHARACTERISTICS parameter level, divide earthquake prediction fracture density grade interval. step also comprises: choose two dimension to cross figure according to the cross data point concentration class of figure of described multiple two dimension; Cross in figure in described two dimension, by well logging FRACTURE CHARACTERISTICS parameter level interval division earthquake prediction fracture density grade interval.
Preferably, describedly cross in figure in described two dimension, also comprised by the fracture density grade interval step of well logging FRACTURE CHARACTERISTICS parameter level interval division earthquake prediction: cross in figure in described multiple two dimension, character numerical value corresponding to the grade of well logging crack is interval, the fracture density value obtained by earthquake prediction by crossing is divided into different interval, obtains the distribution grade of different interval Seismic Fracture density.
Preferably, describedly obtain fracture density step according to described seismic data and also comprise: carry out anisotropic analysis according to seismic data and obtain fracture density.
Preferably, describedly to be comprised by seismic data: prestack road collection data.
Preferably, described according to described well logging crack grade interval, demarcating steps is carried out to described earthquake prediction fracture density after also comprise: carry out well logging FRACTURE CHARACTERISTICS analysis interpretation according to research purpose interval.
Compared with prior art, above-mentioned embodiment of the present invention has the following advantages: for the reservoir of fracture development, define the technical scheme that a set of crack is quantitatively carved, by demarcating the well logging of earthquake prediction fracture density, realize the division of crack rank, be conducive to identifying effective fracture in real work, decrease further work the risk and cost, thus reduce the Primary Stage Data processing cost of opposite joint hole type Reservoir Fracture research, obtain higher income.
Embodiment
Below in conjunction with accompanying drawing, the embodiment of the present invention is described.
Embodiment one
Fig. 1 is a kind of schematic diagram demarcating the method in logging reservoir crack of the present invention, and the method comprises the following steps:
S101: gather well-log information and seismic data;
Receive well-log information and seismic data in work area in this step, obtain well logging FRACTURE CHARACTERISTICS parameter, according to described forecast for seismic data fracture density according to described interpretation of logging data; Specifically realize by following steps:
First, the well logging of same work area is carried out the collection of two class data, and a class is well-log information, and a class is seismic data, and wherein, well-log information directly obtains by well logging exploration, and seismic data realizes by the other seismic data acquisition of well.By well logging interpretation, well logging FRACTURE CHARACTERISTICS parameter is obtained for well-log information, FRACTURE CHARACTERISTICS parameter of logging well under normal circumstances comprises: the parameters such as depth side direction resistivity, depth side direction resistivity ratio, micro resistor FMI(FormationMicroScanner Image is carried out according to described well-log information, FMI) well logging interpretation acquisition pores porosity, fracture porosity, carries out the extraction of above well logging FRACTURE CHARACTERISTICS parameter usually by reservoir prediction and interpretation software LD-EPS reservoir V5.2.Wherein, described fracture density specifically comprises: earthquake ellipticity.Seismic data can specifically can comprise: prestack road collection data and poststack road collection data, under normal circumstances, prestack road collection data can obtain and more retain more original seismic event information, so can above processing procedure by pre-stack seismic inversion and layer description software LD-GMAX
tMreservoir V2.0 realizes, and carry out prestack azimuthal anisotropy analysis according to described seismic wave parameter and obtain described multiple effective fracture density, above processing procedure is by fracture reservoir comprehensive description software LD-FRS
tMfracture V3.2 realizes.
S102: carry out data exchange;
In this step: described well logging FRACTURE CHARACTERISTICS parameter and described earthquake prediction fracture density being carried out crossing obtains two dimension and to cross figure, as shown in Figure 5, horizontal ordinate is chosen for earthquake ellipticity (obtaining from seismic data) both fracture densities, ordinate is chosen for fracture porosity (obtaining from well-log information), in above-mentioned coordinate, two class data are crossed, obtain two dimension to cross figure, simultaneously also can according to described depth side direction resistivity in above step, depth side direction resistivity ratio, factor of porosity, multiple and described earthquake prediction fracture density in fracture porosity is carried out crossing and is obtained multiple two dimension and to cross figure, in other words, fracture density according to obtaining in the parameter of number of different types and seismic data crosses, obtain multiple two dimension to cross figure, different fracture parameters and the corresponding relation of fracture density can be reflected in its each figure
S103: Reservoir Fracture well logging is demarcated;
In above step: cross in figure in described two dimension, by well logging FRACTURE CHARACTERISTICS parameter level interval division earthquake prediction fracture density grade interval, according to described well logging crack grade interval, described earthquake prediction fracture density is demarcated.As shown in Figure 5, fracture density can be divided into I, II, III level, wherein the degree of porosity interval of I level correspondence is 0.320 ~ 0.640, the degree of porosity interval of II level correspondence is 0.080 ~ 0.32, degree of porosity interval corresponding to III level is 0 ~ 0.320, and can according to divide after degree of porosity interval 0.320 ~ 0.640 by well-log information as Fig. 3 is demarcated as I level crack section, simultaneously as described in step S102 by depth side direction resistivity, depth side direction resistivity ratio, factor of porosity, multiple parameter such as fracture porosity and described fracture density are carried out crossing and are obtained multiple two dimension and to cross figure, therefore, also can comprise in above step: choose two dimension to cross figure according to the cross data point concentration class of figure of described multiple two dimension, cross in figure in described two dimension, by well logging FRACTURE CHARACTERISTICS parameter level interval division earthquake prediction fracture density grade interval.Describedly choose the two dimension figure step that crosses also comprise according to the cross data point concentration class of figure of described multiple two dimension: cross in figure in described multiple two dimension, character numerical value corresponding to the grade of well logging crack is interval, the fracture density value obtained by earthquake prediction by crossing is divided into different interval, obtains the distribution grade of different interval Seismic Fracture density.Cross in figure in described multiple two dimension, the data point of figure of described two dimension being crossed is divided into multiple independently data area, cross in described multiple two dimension and compare number of data points in corresponding region between figure and obtain the maximum data point region, the two dimension corresponding according to the described the maximum data point region figure that crosses is chosen for two dimension and crosses figure.Therefore carrying out data exchange, according to crossing in result determining hole porosity grade interval process, the mode of multiple well logging FRACTURE CHARACTERISTICS parameter can adopted, can choose and determine that two dimension crosses figure according to described multiple two dimension figure that crosses; From multiple two dimension crosses figure, wherein choose one can reflect that the two dimension of geological condition crosses figure, aggregation extent by its data is chosen, its data are more assembled, and illustrate that the data that data reflect are corresponding more obvious, the relation that more can clearly state data is clearer and more definite; Described multiple two dimension crosses in figure, can be multiple independently data areas by diagram root, cross in described multiple two dimension and compare data volume in corresponding region between figure and obtain the maximum data amount region, the two dimension corresponding according to the described the maximum data amount region figure that crosses is defined as two dimension and crosses figure, the figure thus the good two dimension of selection data gathering crosses, can improve the accuracy of Fracture evaluation, simultaneously the method achieves and chooses according to strata condition the mode that comparatively " sensitivity " parameter carries out assessing.
Also specifically comprise after above step: carry out well logging FRACTURE CHARACTERISTICS analysis interpretation according to research purpose interval.
Embodiment two
Embodiment when to be that the present invention is a kind of determine that the method for effective Reservoir Fracture detects for Tarim Oilfield block carbonate rock fractured cave: specifically comprise the following steps:
S201: obtain multiple significant depth; First, well-log information is analyzed, the rank in crack is determined according to conventional logging and FMI logging technology, as shown in Figure 2, from FMI image, obvious identifiable design, its crack are intensive or to grow good penetration of fracture position be 21, correspond in the log data of the same degree of depth, thus to realize well logging be 0.65% with fracture porosity in the degree of depth, and determine that crack Grade I level net porosity is 0.65% according to this parameter, restrictive pore degree scope is for being greater than 0.3%.Or can passing hole porosity PHIE (%) and dark direction finding resistivity RLLD (ohmm) data carry out crossing divide crack rank (Fig. 4) by setting up two-dimensional coordinate system, be I, II, III level according to data gathering situation by crack partition of the level.
S202: obtain effective fracture density scope;
On the basis that a large amount of drilling well is demarcated and added up, demarcate the fracture density magnitude range that I, II, III tri-class fracture reservoir are corresponding, as shown in Figure 5, for according to demarcate I, II, III tri-class crack degree of depth respective depth from geological data stratum in fracture density, thus form crossing of data, and the corresponding relation determining well logging crack rank and seismic data crack rank is as shown in table 1:
Table 1
S203: determine effective fracture parameters;
When gathering Reservoir Fracture density for when both ellipticity is 1.3, when falling into the fracture density scope >1.22 of the described degree of depth, determine that effective fracture porosity is 0.3% according to described effective well logging FRACTURE CHARACTERISTICS parameter area fracture porosity >=0.3%, thus determine by effective fracture porosity value, and confirm that the fracture pore in its respective depth stratum is active porosity, have and estimate oil reserve.
In order to quantitative description crack must utilize well logging, earthquake two kinds of information, by correlation analysis and multiple intersection, determine the dimension of fracture reservoir prediction.On this basis, divide I, II, III grade of crack, demarcated by well between the fracture density (ellipticity) that well logging different stage crack and earthquake prediction are obtained and just set up relation, return well logging and can realize three-dimensional objective interval different stage fracture reservoir thickness quantitative description.The method feature is: all well layer positions are carried out analyzing, arrange, comprehensively, namely from well, returned well through earthquake.The I level crack section predicted, necessarily has the I level crack of well logging interpretation, but also may grow II grade or III grade of crack.Prediction is in the past instruction " having crack " just, is at most that the power of fracture density is qualitative.And predict III grade of crack section, then generally there is not I, II grade of crack.Quantitative forecast precision is higher, but is difficult to differentiate thinner crack individual layer as well logging, and does not generally comprise low angle crack.
Through the above description of the embodiments, those skilled in the art can be well understood to the present invention can by hardware implementing, and the mode that also can add necessary general hardware platform by software realizes.Based on such understanding, technical scheme of the present invention can embody with the form of software product, it (can be CD-ROM that this software product can be stored in a non-volatile memory medium, USB flash disk, portable hard drive etc.) in, comprise some instructions and perform method described in each embodiment of the present invention in order to make a computer equipment (can be personal computer, server, or the network equipment etc.).
It will be appreciated by those skilled in the art that accompanying drawing is the schematic diagram of a preferred embodiment, the module in accompanying drawing or flow process might not be that enforcement the present invention is necessary.
It will be appreciated by those skilled in the art that the module in the device in embodiment can carry out being distributed in the device of embodiment according to embodiment description, also can carry out respective change and be arranged in the one or more devices being different from the present embodiment.The module of above-described embodiment can merge into a module, also can split into multiple submodule further.
The invention described above embodiment sequence number, just to describing, does not represent the quality of embodiment.
Be only several specific embodiment of the present invention above, but the present invention is not limited thereto, the changes that any person skilled in the art can think of all should fall into protection scope of the present invention.