CN104678434A - Method for predicting storage layer crack development parameters - Google Patents
Method for predicting storage layer crack development parameters Download PDFInfo
- Publication number
- CN104678434A CN104678434A CN201310631520.8A CN201310631520A CN104678434A CN 104678434 A CN104678434 A CN 104678434A CN 201310631520 A CN201310631520 A CN 201310631520A CN 104678434 A CN104678434 A CN 104678434A
- Authority
- CN
- China
- Prior art keywords
- value
- zone
- interest
- data
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Abstract
The invention provides a method for predicting storage layer crack development parameters, and belongs to the field of petroleum geophysics exploration. The method comprises the following steps that (1) earthquake data is input; (2) earthquake data with different azimuth angles is used for extracting a target layer two-way travel time value, the normalization processing is carried out, and the two-way travel time value of each target layer after the normalization processing is obtained; (3) an earthquake rose figure and an integral crack development intensity plane figure are respectively generated on each point on a set lattice, the two figures are laminated and are then displayed, and a crack parameter plane figure is obtained. When the method provided by the invention is used, a crack type storage layer can be accurately predicted.
Description
Technical field
The invention belongs to field of petroleum geophysical exploration, be specifically related to a kind of method that Predicating Reservoir Fractures grows parameter, belonging to seismic data interpretation category, is that a kind of two way travel time value of the zone of interest of each position angle geological data that utilizes is to calculate the method for the fracture development parameter of reservoir.
Background technology
The carbonate fracture type reservoir relative in southern china exploratory area is grown, and is embodied in the Sichuan Basin, bores and has met a large amount of Fractured Gas Reservoirs, and obtain industrial gas.These slit formation carbonate reservoirs are formed under the effect being subject to tectonic stress substantially, have very strong nonuniformity, and in uneven distribution in plane, exploratory development difficulty is very large, and how identifying crack elimination is one of emphasis of Some Comments On Geophysical Work.So, the prediction of the type reservoir and technology exploration tool are had very important significance.
Find that underground zone of interest is under the state of growing crack according to research, P-wave amplitude, two way travel time value are by having the feature with azimuthal variation during this crack body.Such as, Xiangyang Li " the Fracture detection using azimuthal variation of P-wave moveout from orthogonal seismic survey lines " that " Geophysics " deliver in 1999.Research data shows: 1. the reflected P ripple two way travel time of crack body top layer is not with azimuthal variation, and there is faint change at amplitude position angle, not obvious to anisotropy reflection.2. the reflected P wave table of crack body bottom reveals very strong azimuthal anisotropy, P wave reflection amplitude and two way travel time relevant with survey line and fracture azimuth, when survey line is parallel with crack, amplitude is the strongest, two way travel time is the shortest; Along with the increase of survey line and crack angle, amplitude weakens gradually, two way travel time is elongated gradually; Time vertical with fracture orientation to survey line, amplitude is the most weak, and two way travel time is the longest.As can be seen here, P ripple is by showing very strong azimuthal anisotropy feature during the body of crack, the influence degree of crack to P wave response depends on that the scale of Fracture System is necessarily opened in angle between fracture azimuth and survey line and underground, and this phenomenon is especially larger on impact during the walking of bottom.
The technical method of existing typical crack type reservoir prediction, mainly relevant with poststack and prestack correlation technique.As coherent technique, curvature, P ripple anisotropic analysis technology, the calculating of the road collection time difference etc. realize prediction.Some patents as invention " crack prediction method and device " (application number: 201010205983.4) patent discloses a kind of crack prediction method and device, be mainly pickup center with zone of interest, utilize picking window to obtain the reflection amplitude of each seismic trace; The position angle of acquisition and reflection amplitude is utilized to carry out ellipse fitting, to determine fracture orientation and fracture density; The defining method that patent discloses a kind of Reservoir Fracture of invention " a kind of defining method of Reservoir Fracture " (application number: 201210242061.X), comprise: according to actual seismic data type, the analysis of multiple method FRACTURE PREDICTION is carried out to actual seismic data, obtain crack attribute and the distribution of reservoir; Utilize actual geologic parameter just drilling and FMI data, comparison Reservoir Fracture distributes, and obtains matching degree and the goodness of fit; From multiple FRACTURE PREDICTION analytical approach, a kind of best FRAC method is preferably chosen according to matching degree and the goodness of fit.According to the FRAC method optimized, explain actual seismic data, determine Reservoir Fracture.But the technical method of these routines implements difficulty quite greatly, operates and is subject to the impact of various unfavorable factor, mainly there is following Railway Project:
(1) well data is utilized to carry out the constraint of FRACTURE PREDICTION, often bad in the regional prediction effect away from well.
(2) utilize the attribute such as amplitude, frequency to carry out ellipse fitting and calculate fracture azimuth and intensity, predict the outcome often widely different, be difficult to the judgement these technical methods being accepted or rejected and affects achievement.
(3) to the evaluation of miniature fracture development parameter, convenient technical process is often difficult to accurate description.
Summary of the invention
The object of the invention is to solve the difficult problem existed in above-mentioned prior art, provide a kind of Predicating Reservoir Fractures to grow the method for parameter, adopt relatively simple techniqueflow to solve the identification of fracture type reservoir, more effectively the type reservoir is predicted.
The present invention is achieved by the following technical solutions:
Predicating Reservoir Fractures grows a method for parameter, comprising:
(1) geological data is inputted;
(2) utilize the geological data of different orientations to extract zone of interest two way travel time value, and be normalized, obtain the two way travel time value of each zone of interest after normalized;
(3) earthquake rose diagram and comprehensive fracture development strength plane figure are generated respectively to each point on setting grid, (be calculate through ellipse fitting the earthquake rose diagram formed on each net point on the one hand, utilize the coefficient of dispersion value that net point calculates on the other hand, comprehensive fracture development strength plane figure is obtained again after interpolation calculation) show after superimposed for this two figure, obtain fracture parameters planimetric map.
In described step (1), the geological data of input is the CDP road collection data formed after normal moveout correction.
Described step (2) comprises the following steps:
(21), according to seismic data and actual conditions, to processing the superposition and migration processing that Hou CDP road collection data carry out a point position angle, ranges of incidence angles limits, each point of azimuthal displacement data volume is obtained;
(22), according to the actual conditions of seismic data and the reflectance signature of zone of interest, use identical interpretation parameters to carry out pickup to each point of azimuthal displacement data volume to zone of interest to explain, obtain the two way travel time value of the zone of interest position of each bearing data body, and gridding interpolation and the disposal of gentle filter are carried out to the two way travel time value of the zone of interest position of picking up each bearing data body obtained, obtain the zone of interest two way travel time data after the process of each bearing data body.
Gridding interpolation in described step (22) and the disposal of gentle filter specific as follows:
By utilizing well-log information and each azimuthal displacement data volume to carry out accurate calibration to aboveground crack elimination, and carry out choosing of zone of interest.Zone of interest selection principle is choose for the bottom of crack elimination section.And determine the reflectance signature of zone of interest, as crest, trough or zero phase etc., these easily carry out the reflectance signature of tracing of horizons.Wherein, the process such as layer position gridding interpolation, smothing filtering after explanation is according to being inserted into 1 line X1 road in gridding interpolation algorithm by the two way travel time value of the layer position of explanation grid (usually explaining that grid is greater than 1 line X1 road), and the filtering parameter smoothing filtering certain to the layer position two way travel time data setting after gridding calculating, make it eliminate some and walk duration catastrophe point, closer to the explanation of actual purpose layer.Gridding interpolation algorithm can have least square method, method of weighted mean, (general) Kriging technique etc., uses these computing method should determine according to interpretation horizon concrete condition.
(23), the two way travel time value obtaining each zone of interest after normalized is normalized to the zone of interest two way travel time data after the process of each bearing data body described.Specific as follows:
According to the maximal value of zone of interest two way travel time and the actual conditions of minimum value of each bearing data body, when each position angle is walked by unification, numerical value normalization is in relatively little n ~ m codomain.Normalized that is be utilize add, subtract, the combinatorial operation such as multiplication and division, make two way travel time value all have diminishing and making that the maximum two way travel time value on each point is relative with the multiple of minimum two way travel time value to be increased to a certain extent, be beneficial to carry out ellipse fitting calculating.
Described step (3) comprises the following steps:
(31), the mesh parameter of setting grid is set according to actual conditions; Described mesh parameter comprises mesh spacing and grid number;
(32), extract the two way travel time value setting described each zone of interest after normalized on net point, and coefficient of dispersion calculating is carried out to this point; Using the coefficient of dispersion on this aspect as the comprehensive fracture development intensive parameter of this point;
(33), calculate fracture strike on setting net point and the crack relative intensity of this trend, obtain earthquake rose diagram;
(34), utilize the value of the coefficient of dispersion on setting net point after interpolation calculation, generate comprehensive fracture development strength plane figure;
(35), by described comprehensive fracture development strength plane figure and described earthquake rose diagram be superimposed on a sheet of planar figure and show, thus obtain fracture parameters planimetric map.
That the matching of employing three null ellipse is carried out calculating in described step (33), specific as follows:
Oval short-axis direction is fracture strike, and the ratio of its major axis and minor axis is as the crack relative intensity of this trend; These parameters form the earthquake rose diagram on this net point, and the earthquake rose diagram comprehensively on each net point obtains the distribution situation of fracture strike in plane on each net point and crack relative intensity.
Described step (34) is achieved in that
By the value of the coefficient of dispersion on described setting net point after gridding interpolation calculation, form new grid and interpolation point, and isoplethes drawing is carried out to these data points or forms color filling display, thus obtain comprehensive fracture development strength plane distribution plan.
Described step (35) is achieved in that
A sheet of planar figure arranges two Z values, the value of comprehensive fracture development intensity in a corresponding comprehensive fracture development strength plane figure of Z value, the scope adjusting this value can control the display effect of comprehensive fracture development intensity; The value of crack relative intensity in the corresponding earthquake rose diagram of another Z value, the scope adjusting this value can control the display effect of crack relative intensity value in each earthquake rose diagram; Wherein, the X on superimposed planimetric map, Y-coordinate with in comprehensive fracture development strength plane figure and described earthquake rose diagram two figure X, Y-coordinate be the same.
Compared with prior art, the invention has the beneficial effects as follows: the crack elimination distributed areas that the technology of the present invention is predicted than other method are more clear, accurate and understand, relative to traditional crack elimination Forecasting Methodology, can detect by fracture reservoir exactly.Explore certain block for south to illustrate, bore the slit formation carbonate reservoir of meeting when drilling well and can obtain middle and high product industrial gas after tested.Known to being disclosed by the known drilling well in study area, the achievement of invention and drilling data match result better, can go out crack elimination by Accurate Prediction.
Accompanying drawing explanation
Fig. 1 is the step block diagram of the inventive method;
Fig. 2 is the two way travel time value of zone of interest on each well point of example one of the present invention;
Fig. 3 is the situation of inventive method on each well point of example one of the present invention and conventional method and aboveground average crack density value, formation testing gas data;
Fig. 4 is example two survey line zone of interest two way travel time value of the present invention and normalized value thereof;
Fig. 5 is example two survey line fitted ellipse of the present invention and fracture strike relation schematic diagram;
Fig. 6 is the crack vector rose diagram of example two of the present invention design;
Fig. 7 is zone of interest two way travel time value and normalized value thereof on example three survey line of the present invention;
Fig. 8 is example three fitted ellipse of the present invention and known fracture strike schematic diagram.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail:
The present invention is that a kind of coefficient of dispersion of the zone of interest two way travel time value calculating different orientations geological data and ellipse fitting of utilizing calculates thus realize the technology of the fracture development parameter of predicting reservoir, can the distribution situation exactly in fracture type reservoir plane effectively predict.
Crack elimination forecasting techniques provided by the invention comprises: obtain a point position angle superposition, the rear geological data of skew, and carried out zone of interest explanation again after carrying out well composite seismogram calibration to these geological datas respectively; Then set mesh parameter, utilize each position angle zone of interest two way travel time value of grid each point in plane to carry out coefficient of dispersion and the ellipse fitting calculating of this point; The coefficient of dispersion data of recycling setting grid each point, calculate after interpolation through gridding and generate isoline planimetric map, thus obtain comprehensive fracture development strength plane distributed data; Secondly, each earthquake rose diagram plane set on net point is calculated through ellipse fitting; Both superimposed displays comprehensive analysis, thus realize the prediction to the fracture development parameter of reservoir in plane.Coefficient of dispersion is also called the coefficient of variation, is the conventional statistical indicator in the middle of statistics, is mainly used in the dispersion degree of the variable series comparing varying level and the representativeness of average.
The present invention is the developmental state utilizing zone of interest two way travel time value fracture detection type reservoir, can cause anisotropy but not reservoir then can not produce anisotropic feature by azimuthal geological data after usual crack elimination oily to the geological data of different orientations.Concrete upper can show zone of interest two way travel time value generation difference on seismic location angular data body, this walk time coefficient of dispersion and the comprehensive fracture development intensity of reservoir be often proportionate.So the coefficient of dispersion value calculating the two way travel time of zone of interest just can the comprehensive fracture development intensity of Accurate Prediction crack elimination.The crack relative intensity value in different fracture azimuth on different net point and this direction is calculated again by ellipse fitting.Again two kinds of achievements are carried out to the superimposed display in plane, so just can judge plane distribution and the fracture azimuth of crack elimination exactly.
The present invention is a kind of technology carrying out fracture detection type reservoir development parameter utilizing the zone of interest two way travel time value of different orientations, can reach and predict the fracture parameters of reservoir accurately and efficiently.By dividing superposition, the migration processing at position angle and incident angle, the response characteristic when bottom utilizing well data to extract the corresponding same fracture reservoir of different orientations geological data is walked, and according to these features, each geological data is carried out to the explanation of zone of interest; Secondly, carry out coefficient of dispersion to the different orientations zone of interest two way travel time value on setting net point and calculate and interpolation generation isoline planimetric map, this planimetric map is exactly the comprehensive fracture development intensity map of fracture reservoir.Secondly, the two way travel time value on grid each point after normalization is utilized to carry out ellipse fitting calculating, obtain the earthquake rose diagram of fracture azimuth on grid each point and this direction relative intensity value, the earthquake rose diagrams on these setting net points just form a tension gash and grow direction and this direction relative intensity planimetric map.Be presented on a sheet of planar figure by superimposed for these two planimetric maps, the fracture reservoir that just can dope different fracture development intensity and fracture strike grows the region roughly distributed, and can increase substantially the exploration economic benefit of crack elimination.
As shown in Figure 1, the embodiment of the inventive method is:
1. geological data input.
2. to utilizing the geological data of different orientations to extract zone of interest two way travel time value and normalized after step 1.
3. carry out earthquake rose diagram and coefficient of dispersion respectively to each point on the setting grid after step 2 to calculate, comprehensive and complete fracture parameters planimetric map.
Wherein, geological data inputs, and comprises the following steps:
I, pre-service is carried out to the longitudinal wave earthquake data of field acquisition.Pre-service is mainly for the common seismic relative amplitude preserved processing work in early stage, and road collection forms normal moveout correction CDP road collection data again after normal moveout correction.
If the single bin superposition of ii can not meet the requirement of seismic data signal to noise ratio (S/N ratio), the process of extended surface unit can be carried out the single bin before superposition, form grand bin, to improve stacking fold, increase the signal to noise ratio (S/N ratio) of pre-stack data.Grand bin extracts the process of the bin that is exactly that to be pre-stack seismic road collection on several adjacent depth points by the prestack road collection data pick-up on single depth point be.The distance how many (i.e. the size of grand bin) of neighboring track collection, can be selected according to actual conditions.
Wherein, utilize the geological data of different orientations to extract zone of interest two way travel time value and normalized, comprise the following steps:
I, to process after normal moveout correction road collection data carry out a point position angle, ranges of incidence angles restriction superposition and migration processing.The division of position angle number is moved towards situation according to seismic data and study area fracture development and is determined, and can carry out decile or not decile position angle division in principle, limits ranges of incidence angles aspect and will ensure that each the azimuthal maximum offset divided is consistent or close.Setting recording geometry direction in process with direct north is 0 °, and clockwise direction rotates, and counts 360 °.Based on symmetry principle, 360 ° of orientation of the geological data of field acquisition shot point-geophone station are converted into 180 ° of position angles, concerning some azimuth coverages, calculate its central angle, central angle then represents this azimuth coverage of division.Computing formula is as follows:
θ in formula (1)
ifor the central angle of i-th superposition of data of design,
for i-th minimum position angle of superposition of data of design,
for i-th maximum position angle of superposition of data of design, i>=2.
It is the scope will getting certain offset distance that ranges of incidence angles limits principle, and the size of incident angle should according to the principle of AVO, and for zone of interest, the incident angle of maximum offset can not be greater than 30 °.Relevant incidence angle computing formula is as follows:
θ=arctgD／2h (2)
In formula (2), D is the earthquake data offset distance of design, and θ is the incident angle of design, and h is prediction zone of interest buried depth.
The stacking fold of each position angle superposition of data after ii, point position angle and ranges of incidence angles limit should be unanimous on the whole, to ensure the noise when balancing energy of section.The setting range of incident angle is-be not more than 30 ° in principle for its maximum incident angle of zone of interest, and minimum incident angle should be greater than 5 °.Azimuthal number should be more than or equal to 3 in principle, and concrete visual actual conditions design accordingly.
Iii, utilizing well-log information and each azimuthal displacement data volume to carry out accurate calibration to aboveground crack elimination, and carry out choosing of zone of interest, zone of interest selection principle is choose for the bottom of crack elimination section.And determine the reflectance signature of zone of interest, as crest, trough or zero phase etc., these easily carry out the reflectance signature of tracing of horizons.
The reflectance signature of the zone of interest that iv, utilization are chosen and identical interpretation parameters carry out horizon picking explanation respectively to each azimuthal displacement data volume, interpretation parameters comprises start-finish line, Taoist monastic name and setting grid, reflectance signature, and carries out process such as gridding calculating interpolation, smothing filtering etc. to zone of interest position.Gridding calculates interpolation, smothing filtering etc. and is treated to the two way travel time value of the layer position by explaining grid (usually explaining that grid is greater than 1 line X1 road) according to being inserted into 1 line X1 road in gridding interpolation algorithm, and the filtering parameter smoothing filtering certain to the layer position two way travel time data setting after gridding calculating, make it eliminate some and walk duration catastrophe point, closer to the explanation of actual purpose layer.Gridding interpolation algorithm can have least square method, method of weighted mean, universal Kriging etc., uses these computing method should determine according to actual conditions.
V, each zone of interest two way travel time data to be normalized, by its numerical value normalization to n ~ m codomain.Normalized that is be utilize add, subtract, the combinatorial operation such as multiplication and division, two way travel time value is made all to have diminishing and making that the maximum two way travel time value on each point is relative with the multiple of minimum two way travel time value to be increased to a certain extent, be beneficial to carry out ellipse fitting calculating, each azimuthal walk duration be the corresponding oval position angle in coordinate system-walk duration of calculating.X in coordinate system, Y-axis are Value Data when walking after representing normalization, also represent azimuth value.The setting of azimuth value scope refer to when Y-axis positive axis (azimuth value is 0 °) revolve in the direction of the clock turn around time, azimuth value scope is 0 °-360 °.That is corresponding two data values of each data point in this coordinate system-be walk duration size and azimuth value after normalization respectively, length value between this point and initial point is walk duration after normalization, and is azimuth value with the line of initial point and the size of Y-axis positive axis angle.Again the ellipse fitting that each point on this coordinate is correlated with is calculated.Specifically as shown in Figure 5.
Wherein, on setting grid, each point carries out earthquake rose diagram and coefficient of dispersion respectively and calculates, comprehensive and complete fracture parameters planimetric map, comprises the following steps:
I, arrange mesh parameter according to actual conditions, this grid is called setting grid.Mesh parameter comprises mesh spacing and grid number etc., the large I of mesh parameter needs situation to determine according to the grid distribution of the crack point that will predict and precision, usually the mesh spacing of setting is larger in principle, then the precision predicted is relatively low, can lose some characteristic informations; Mesh spacing is less, then the precision predicted is relatively high, and institute's result of painting is more detailed.The size of setting grid should need situation according to actual and FRACTURE PREDICTION and determine, and the grid of setting is normally well-behaved.
The two way travel time value of each zone of interest on ii, extraction setting net point after normalized
and coefficient of dispersion calculating is carried out to this point.Coefficient of dispersion on this aspect can be used as the comprehensive fracture development intensive parameter of this point.
Coefficient of dispersion computing formula in previous step (32) is as follows:
V in formula (3) ~ (5)
irepresent i-th some coefficient of dispersion in aspect, σ
irepresent i-th some standard deviation,
represent the arithmetic mean of each position angle zone of interest two way travel time value of i-th point,
represent a jth position angle zone of interest two way travel time value of i-th point, k is a point position angle number, j≤k.
Iii, to the trend in crack on setting net point and the crack relative intensity of this trend, adopt three null ellipse matchings to calculate, oval short-axis direction is generally the trend in crack, and the ratio of its major axis and minor axis can be used as this crack, direction relative intensity value.These parameters can form the content of earthquake rose diagram, and this graphic designs should have vector characteristic.Fracture strike polar plot on each setting net point comprehensive thus the earthquake rose diagram of the fracture strike obtained in plane on each net point and relative intensity.In principle, fracture azimuth can be understood as from the north pointer direction line of certain point on crack, according to the horizontal angle between clockwise direction to fracture strike line, position angle and the relative intensity of the fracture strike of the upper prediction of earthquake rose diagram design are all symmetrical in 180 °, and add using this position angle, subtract 5 ° as the fracture strike error ranges predicted.The dial plate of earthquake rose diagram also adopts the span of 0 ° ~ 360 °, dial plate has two parameter coordinate axis, is circular angle coordinate (position angle in reflection crack) and linear relative intensity value coordinate (relative intensity value in reflection crack) respectively.
Iv, on the setting net point in plane coefficient of dispersion value generate isoline planimetric map or color filling display.The method that computer drawing isogram is conventional has triangulation and grid method.Be specially the coefficient of dispersion value on setting net point after gridding interpolation calculation, form new grid and interpolation point, and isoplethes drawing is carried out to these data points or forms color filling display, thus obtain comprehensive fracture development strength plane distribution plan.
V, the fracture strike on comprehensive fracture development strength plane figure and setting net point and the earthquake rose diagram of crack relative intensity to be superimposed on a sheet of planar figure, thus to complete the fracture parameters planimetric map of reservoir.The concrete grammar of superimposed display is, a sheet of planar figure arranges two Z values, and a Z value is comprehensive fracture development intensity level, adjusts the display effect of the comprehensive fracture development intensity of this value range-controllable system; Another Z value is crack relative intensity value, adjusts the display effect of each rose diagram crack, crack relative intensity value on this value range-controllable grid.Wherein, X, the Y-coordinate of the X on superimposed planimetric map, Y-coordinate and these two planimetric maps are the same.
Example of the present invention is:
According to this invention techniqueflow (Fig. 1), formulate job step, example one carries out fracture development parameter prediction to the Daanzhai member shelly limestone in certain three-dimensional work area.
To step 1. in pre-service first pre-service is carried out to the three-dimensional P wave data of field acquisition, form the Hong Mianyuandao collection after normal moveout correction, bin size expands 40mX40m to by 20mX20m, and after grand bin process, the signal to noise ratio (S/N ratio) of each azimuthal stacked section is enhanced.
To step 2. in Hong Mianyuandao collection data carry out a point position angle overlap-add procedure, central angle be respectively 15 °, 60 °, 120 °, 165 ° amount to four position angle seismic data volume-azimuth coverages be followed successively by 0 ° ~ 30 °, 30 ° ~ 90 °, 90 ° ~ 150 °, 150 ° ~ 180 °, the zone of interest buried depth 2400m of design, the geophone offset scope of superposition is defined as 420m ~ 2500m, in the ranges of incidence angles allowed, respectively migration processing is carried out to each position angle superposition of data, reflection wave groups is playbacked.
Utilize well-log information carry out synthetic record and demarcate respectively to each azimuthal displacement data, determine the crack section of reservoir in each azimuthal displacement data, and zone of interest is picked up.That chooses carbonate fracture Reservoir Section zone of interest in example is reflected into zero phase, and this zone of interest is the bottom reflection of carbonate fracture section, so will make an explanation along zero phase.Layer position interpretation parameters is for from 1880 Xian Ji 300 roads, terminate to 3600 lines and 1800 roads, actual conditions setting according to seismic data explains that grid is 10 line X10 roads, afterwards interpolation, round and smooth process are carried out to layer position, form the layer bit data that grid is 1 line X1 road, and in each azimuthal displacement data, use identical interpretation parameters to carry out whole district's explanation respectively to this zone of interest.
To step 3. in plane on setting grid on each point extract each position angle zone of interest two way travel time value respectively
utilize coefficient of dispersion to calculate its comprehensive fracture development intensity and utilize ellipse fitting to calculate the earthquake rose diagram set on net point.And gridding interpolation calculation is done to the comprehensive intensity value on setting net point, obtain the comprehensive fracture development strength plane distribution plan in this district.Earthquake rose diagram on comprehensive fracture development strength plane figure and setting net point is carried out superimposed display, thus obtains the flat distribution map of work area Daanzhai member fracture reservoir.The setting mesh spacing of example setting is 5 road X5 lines, and new grid when calculating comprehensive fracture development intensity is 2 road X2 lines.
From Fig. 2, produce difference when crack makes the zone of interest of bottom walk, thus produce discrete phenomena.In Fig. 3, it is the interpretation of result of this invention and original routine techniques comparison, production capacity and the average crack density value of the value using conventional P ripple anisotropy technology to predict and the fracture reservoir on well point are corresponding bad, and the production capacity of well point, average crack density value and coefficient of dispersion value anastomose property are good, illustrate that this invention technical method is effective.
Example two utilizes the two way travel time difference of two-dimentional survey line to come Prediction of fracture direction and relative intensity value to certain laboratory data of employing.
The position angle of survey line 1 is 15 °, and the position angle of survey line 2 is 60 °, and survey line 4 position angle is 150 °.This model is made up of a Fractured gassiness layer of sand and a position angle isotropy shale superstratum, and two way travel time value and the normalized value thereof of Fig. 4 to be the offset distance of these three lines be 2000m, visible difference is obvious.Fracture strike (short-axis direction) that three null ellipse the Fitting Calculation draw is good with actual fracture strike consistance in order to utilize normalized value to carry out for Fig. 5, the fracture strike that Fig. 6 is the oval prediction of utilization and the earthquake rose diagram designed by relative intensity value, visible earthquake rose diagram is clear.
Example three is to adopting the real data in oil field district, the North Sea to utilize the two way travel time difference (as shown in Figure 7) of two-dimentional survey line (survey line 1, survey line 3, survey line 4) to come Prediction of fracture direction.The minor axis of the ellipse in Fig. 8 is consistent with actual fracture strike.
Technique scheme is one embodiment of the present invention, for those skilled in the art, on the basis that the invention discloses application process and principle, be easy to make various types of improvement or distortion, and the method be not limited only to described by the above-mentioned embodiment of the present invention, therefore previously described mode is just preferred, and does not have restrictive meaning.
Claims (7)
1. Predicating Reservoir Fractures grows a method for parameter, it is characterized in that: described method comprises:
(1) geological data is inputted;
(2) utilize the geological data of different orientations to extract zone of interest two way travel time value, and be normalized, obtain the two way travel time value of each zone of interest after normalized;
(3) earthquake rose diagram and comprehensive fracture development strength plane figure are generated respectively to each point on setting grid, show after superimposed for this two figure, obtain fracture parameters planimetric map.
2. Predicating Reservoir Fractures according to claim 1 grows the method for parameter, it is characterized in that: in described step (1), the geological data of input is the CDP road collection data formed after normal moveout correction.
3. Predicating Reservoir Fractures according to claim 2 grows the method for parameter, it is characterized in that: described step (2) comprises the following steps:
(21), according to seismic data and actual conditions, to processing the superposition and migration processing that Hou CDP road collection data carry out a point position angle, ranges of incidence angles limits, each point of azimuthal displacement data volume is obtained;
(22), according to the actual conditions of seismic data and the reflectance signature of zone of interest, use identical interpretation parameters to carry out pickup to each point of azimuthal displacement data volume to zone of interest to explain, obtain the two way travel time value of the zone of interest position of each bearing data body, and gridding interpolation and the disposal of gentle filter are carried out to the two way travel time value of the zone of interest position of picking up each bearing data body obtained, obtain the zone of interest two way travel time data after the process of each bearing data body;
(23), the two way travel time value obtaining each zone of interest after normalized is normalized to the zone of interest two way travel time data after the process of each bearing data body described.
4. Predicating Reservoir Fractures according to claim 3 grows the method for parameter, it is characterized in that: described step (3) comprises the following steps:
(31), the mesh parameter of setting grid is set according to actual conditions; Described mesh parameter comprises mesh spacing and grid number;
(32), extract the two way travel time value setting described each zone of interest after normalized on net point, and coefficient of dispersion calculating is carried out to this point; Using the coefficient of dispersion on this aspect as the comprehensive fracture development intensive parameter of this point;
(33), calculate fracture strike on setting net point and the crack relative intensity of this trend, obtain earthquake rose diagram;
(34), utilize the value of the coefficient of dispersion on setting net point after interpolation calculation, generate comprehensive fracture development strength plane figure;
(35), by described comprehensive fracture development strength plane figure and described earthquake rose diagram be superimposed on a sheet of planar figure and show, thus obtain fracture parameters planimetric map.
5. Predicating Reservoir Fractures according to claim 4 grows the method for parameter, it is characterized in that: be that the matching of employing three null ellipse is carried out calculating in described step (33), specific as follows:
Oval short-axis direction is fracture strike, and the ratio of its major axis and minor axis is as the crack relative intensity of this trend; These parameters form the earthquake rose diagram on this net point, and the earthquake rose diagram comprehensively on each net point obtains the distribution situation of fracture strike in plane on each net point and crack relative intensity.
6. Predicating Reservoir Fractures according to claim 5 grows the method for parameter, it is characterized in that: described step (34) is achieved in that
By the value of the coefficient of dispersion on described setting net point after gridding interpolation calculation, form new grid and interpolation point, and isoplethes drawing is carried out to these data points or forms color filling display, thus obtain comprehensive fracture development strength plane distribution plan.
7. Predicating Reservoir Fractures according to claim 6 grows the method for parameter, it is characterized in that: described step (35) is achieved in that
A sheet of planar figure arranges two Z values, the value of comprehensive fracture development intensity in a corresponding comprehensive fracture development strength plane figure of Z value, the scope adjusting this value can control the display effect of comprehensive fracture development intensity; The value of crack relative intensity in the corresponding earthquake rose diagram of another Z value, the scope adjusting this value can control the display effect of crack relative intensity value in each earthquake rose diagram; Wherein, the X on superimposed planimetric map, Y-coordinate are the same with the X in comprehensive fracture development strength plane figure and described earthquake rose diagram two figure, Y-coordinate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310631520.8A CN104678434B (en) | 2013-12-02 | 2013-12-02 | Method for predicting storage layer crack development parameters |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310631520.8A CN104678434B (en) | 2013-12-02 | 2013-12-02 | Method for predicting storage layer crack development parameters |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104678434A true CN104678434A (en) | 2015-06-03 |
| CN104678434B CN104678434B (en) | 2017-05-24 |
Family
ID=53313754
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310631520.8A Active CN104678434B (en) | 2013-12-02 | 2013-12-02 | Method for predicting storage layer crack development parameters |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104678434B (en) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105332685A (en) * | 2015-11-13 | 2016-02-17 | 西南石油大学 | Propping agent multistage paving method for improving coal bed complicated crack supporting effect |
| CN106154326A (en) * | 2016-08-18 | 2016-11-23 | 中国石油天然气集团公司 | The method and device that a kind of vertical buckling fold fracture spacing is evaluated |
| CN107728214A (en) * | 2016-08-11 | 2018-02-23 | 中国石油化工股份有限公司 | A kind of crack prediction method |
| CN107817517A (en) * | 2016-09-12 | 2018-03-20 | 中国石油化工股份有限公司 | Trace gather preprocess method and device for azimuthal AVO inverting |
| CN108761520A (en) * | 2018-05-08 | 2018-11-06 | 中石化石油工程技术服务有限公司 | A kind of method of seismic acquisition parameter design |
| CN109425900A (en) * | 2017-09-05 | 2019-03-05 | 中国石油化工股份有限公司 | A kind of Seismic Reservoir Prediction method |
| CN110390702A (en) * | 2019-07-09 | 2019-10-29 | 中国石油化工股份有限公司 | The change face element plane interpolation method of analysis is moved towards based on geology |
| CN111399049A (en) * | 2020-04-29 | 2020-07-10 | 西南石油大学 | Crack strength prediction method based on data volume dimensionality reduction and discrete coefficient calculation |
| CN111506861A (en) * | 2020-04-24 | 2020-08-07 | 成都捷科思石油天然气技术发展有限公司 | Method for calculating crack strength of favorable region of target layer |
| CN111505715A (en) * | 2020-04-26 | 2020-08-07 | 成都捷科思石油天然气技术发展有限公司 | Method for calculating crack parameters based on central incidence angle of depth domain |
| CN112198551A (en) * | 2020-09-12 | 2021-01-08 | 北京恺标技术发展有限公司 | Universal novel method for quantitatively evaluating reservoir fracture strength |
| CN113970785A (en) * | 2020-07-22 | 2022-01-25 | 中国石油化工股份有限公司 | Method and system for predicting development of underground river crack, storage medium and electronic equipment |
| CN116520419A (en) * | 2023-05-05 | 2023-08-01 | 中国矿业大学(北京) | Hot fluid crack channel identification method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102012521A (en) * | 2010-10-27 | 2011-04-13 | 中国石油化工股份有限公司 | Method for detecting pre-stack cracks in seismic reservoir prediction |
| CN102213769A (en) * | 2010-04-07 | 2011-10-12 | 中国石油天然气集团公司 | Method for determining anisotropic parameters by utilizing data of three-dimensional VSP (Vertical Seismic Profile) |
| CN102540250A (en) * | 2010-12-08 | 2012-07-04 | 同济大学 | Azimuth fidelity angle domain imaging-based fractured oil and gas reservoir seismic exploration method |
-
2013
- 2013-12-02 CN CN201310631520.8A patent/CN104678434B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102213769A (en) * | 2010-04-07 | 2011-10-12 | 中国石油天然气集团公司 | Method for determining anisotropic parameters by utilizing data of three-dimensional VSP (Vertical Seismic Profile) |
| CN102012521A (en) * | 2010-10-27 | 2011-04-13 | 中国石油化工股份有限公司 | Method for detecting pre-stack cracks in seismic reservoir prediction |
| CN102540250A (en) * | 2010-12-08 | 2012-07-04 | 同济大学 | Azimuth fidelity angle domain imaging-based fractured oil and gas reservoir seismic exploration method |
Non-Patent Citations (3)
| Title |
|---|
| XIANG-YANG LI: "Fracture detection using azimuthal variation of P-wave moveout from orthogonal seismic survey lines", 《GEOPHYSICS》 * |
| 孙炜等: "P波各向异性裂缝预测技术在ZY区碳酸盐岩储层中的应用", 《石油与天然气地质》 * |
| 董渊等: "利用P波层间时差确定裂缝性地层的各向异性参数", 《石油地球物理勘探》 * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105332685A (en) * | 2015-11-13 | 2016-02-17 | 西南石油大学 | Propping agent multistage paving method for improving coal bed complicated crack supporting effect |
| CN105332685B (en) * | 2015-11-13 | 2018-05-25 | 西南石油大学 | A kind of proppant multistage spread method for improving coal seam complex fracture support effect |
| CN107728214B (en) * | 2016-08-11 | 2019-04-09 | 中国石油化工股份有限公司 | A kind of crack prediction method |
| CN107728214A (en) * | 2016-08-11 | 2018-02-23 | 中国石油化工股份有限公司 | A kind of crack prediction method |
| CN106154326A (en) * | 2016-08-18 | 2016-11-23 | 中国石油天然气集团公司 | The method and device that a kind of vertical buckling fold fracture spacing is evaluated |
| CN107817517B (en) * | 2016-09-12 | 2019-10-08 | 中国石油化工股份有限公司 | Trace gather preprocess method and device for azimuthal AVO inverting |
| CN107817517A (en) * | 2016-09-12 | 2018-03-20 | 中国石油化工股份有限公司 | Trace gather preprocess method and device for azimuthal AVO inverting |
| CN109425900A (en) * | 2017-09-05 | 2019-03-05 | 中国石油化工股份有限公司 | A kind of Seismic Reservoir Prediction method |
| CN108761520A (en) * | 2018-05-08 | 2018-11-06 | 中石化石油工程技术服务有限公司 | A kind of method of seismic acquisition parameter design |
| CN110390702A (en) * | 2019-07-09 | 2019-10-29 | 中国石油化工股份有限公司 | The change face element plane interpolation method of analysis is moved towards based on geology |
| CN110390702B (en) * | 2019-07-09 | 2022-11-08 | 中国石油化工股份有限公司 | Variable surface element plane interpolation method based on geological trend analysis |
| CN111506861A (en) * | 2020-04-24 | 2020-08-07 | 成都捷科思石油天然气技术发展有限公司 | Method for calculating crack strength of favorable region of target layer |
| CN111505715A (en) * | 2020-04-26 | 2020-08-07 | 成都捷科思石油天然气技术发展有限公司 | Method for calculating crack parameters based on central incidence angle of depth domain |
| CN111399049A (en) * | 2020-04-29 | 2020-07-10 | 西南石油大学 | Crack strength prediction method based on data volume dimensionality reduction and discrete coefficient calculation |
| CN113970785A (en) * | 2020-07-22 | 2022-01-25 | 中国石油化工股份有限公司 | Method and system for predicting development of underground river crack, storage medium and electronic equipment |
| CN112198551A (en) * | 2020-09-12 | 2021-01-08 | 北京恺标技术发展有限公司 | Universal novel method for quantitatively evaluating reservoir fracture strength |
| CN116520419A (en) * | 2023-05-05 | 2023-08-01 | 中国矿业大学(北京) | Hot fluid crack channel identification method |
| CN116520419B (en) * | 2023-05-05 | 2024-01-30 | 中国矿业大学(北京) | Hot fluid crack channel identification method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104678434B (en) | 2017-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104678434B (en) | Method for predicting storage layer crack development parameters | |
| CN103454685B (en) | Study of The Impedence Inversion Restrained By Well Log is utilized to predict the method and apparatus of sand thickness | |
| Wang et al. | Current developments on micro-seismic data processing | |
| CN102540250B (en) | Azimuth fidelity angle domain imaging-based fractured oil and gas reservoir seismic exploration method | |
| CN101551466B (en) | Method for improving prediction precision of oil and gas reservoir by using seismic attribute related to offset distance | |
| CN109425896A (en) | Dolomite oil and gas reservoir distribution forecasting method and device | |
| CN102466815B (en) | Triassic clastic rock petroleum reservoir identification method | |
| CN104853822A (en) | Method for evaluating shale gas reservoir and searching sweet spot region | |
| CN104636980B (en) | Collect the geophysics characterizing method of condition for channel reservoir type oil gas | |
| CN106154323A (en) | The thin method for predicting reservoir of phased stochastic inverse that frequency processes is opened up based on earthquake | |
| CN105510993A (en) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method | |
| CN102213769A (en) | Method for determining anisotropic parameters by utilizing data of three-dimensional VSP (Vertical Seismic Profile) | |
| CN112883564B (en) | Water body temperature prediction method and prediction system based on random forest | |
| CN103675907A (en) | AVO inversion hydrocarbon detection method based on petrographic constraints | |
| CN104297784A (en) | Primary wave azimuthal anisotropy based fracture predicting method | |
| CN103257363B (en) | A kind of method of fracture dip in Underground fracture-type reservoir | |
| CN103728659A (en) | Method for improving underground karst detecting precision | |
| CN104675392A (en) | Reservoir lithology identification method based on pre-stack multi-parameter dimensionality reduction | |
| CN102375154A (en) | Wide azimuth three-dimensional earthquake-based fracture parameter determining method | |
| CN104155701A (en) | Multi-scale fracture predicting method by adopting pre-stacking seismic data and well information | |
| CN103869362A (en) | Method and equipment for obtaining body curvature | |
| CN106405639B (en) | A kind of inversion method of pre-stack seismic reservoir lithology parameter | |
| CN102879815B (en) | Based on the structural earthquake attributes extraction method of spatial mode gray scale adjoint matrix | |
| CN107942405A (en) | The method for predicting thin sand-mud interbed sand body cumulative thickness | |
| CN105445787A (en) | Crack prediction method for preferred orientation daughter coherence |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |