CN103376463A - Inversion modeling method based on fault control - Google Patents
Inversion modeling method based on fault control Download PDFInfo
- Publication number
- CN103376463A CN103376463A CN201210110254XA CN201210110254A CN103376463A CN 103376463 A CN103376463 A CN 103376463A CN 201210110254X A CN201210110254X A CN 201210110254XA CN 201210110254 A CN201210110254 A CN 201210110254A CN 103376463 A CN103376463 A CN 103376463A
- Authority
- CN
- China
- Prior art keywords
- line
- zone
- special processing
- grid
- rent
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses an inversion modeling method based on fault control, and belongs to the field of oil gas and coal bed gas earthquake exploration and development. The inversion modeling method includes (1) meshing a seismic section; (2) defining a control line; (3) finding an intersection of a fault and the control line; (4) making a vertical line; (5) defining mesh special processing areas; (6) defining mesh properties of the mesh special processing areas; (7) carrying out interpolation calculation on the split mesh special processing areas by utilizing a wave impedance curve of well log data; and (8) carrying out vertical interpolation calculation on the non-mesh special processing areas by utilizing the wave impedance curve of well log data. By utilizing the inversion modeling method based on fault control, a more elaborating initial model is obtained, earthquake model data, which is more reliable and closer to a real earthquake structure, is provided for a variety of inversion methods, solution multiplicity of inversion is reduced, and mathematic and physical method employed in an inversion stage are guaranteed to have geological meanings.
Description
Technical field
The invention belongs to oil gas and coalbed gas seismic exploration and development field, be specifically related to a kind of inverting modeling method based on faults control.
Background technology
Modeling in the seismic inversion is the wide band surge impedance model that comprehensive earthquake, decryption and well-log information are set up.When modeling the decryption that relies on is had relatively high expectations, the personnel that generally explain are familiar with to set up model according to self to the structure in somewhere, and are just drilling, but it is just drilling the result and seismic data has error.Need to reduce gradually error by constantly revising model.Finally obtain satisfied result.
The impedance initial value model that utilizes the mode of interpolation to carry out the overall situation on the basis of setting up the other primary earthquake wave impedance of well among the commercial Inversion Software Strata is set up.Concrete grammar is as follows: at first the correlativity by seismic trace near well and composite traces stretches longitudinally to logging trace and compresses, and just can obtain the impedance initial value of well when related coefficient reaches certain standard; Under the constraint of seismic horizon, choose suitable interpolation method the impedance initial value of well is carried out interpolation and extrapolation, set up overall impedance initial value model; Then take this model as initial point, adopt corresponding inversion method to carry out seismic inversion.
Utilize the coordinate information of geological data among the commercial Inversion Software Jason, fully the data such as combined structure, geology, deposition and well logging form the closed model of a parameterized time, Depth Domain.This model is the basis of reservoir, hydrocarbon-bearing pool quantitative description.This parameterized model has comprised layer a position, tomography, formation contact (integration/unconformability, river course/reef etc.), logging trace and based on the weight coefficient of layer position.Utilize this closed model to create the three-dimensional properties model (such as in length and breadth wave impedance, factor of porosity, density etc.) of different logging traces.
Tomography is significant in fields such as Solid Mineral exploration, oil-gas exploration, earthquake forecasts, but all fail to take full advantage of seismic profile interpretation information in the modeling method in the above-mentioned commercial Inversion Software, the model interrupting layer information response that they are set up is not obvious, and the zone of interest continuity is bad.
Summary of the invention
The object of the invention is to solve a difficult problem that exists in the above-mentioned prior art, a kind of inverting modeling method based on faults control is provided, be used for generating meticulousr initial model, for various inversion methods provide the reliable seismic model data of more pressing close to real construction that reach, reduce the multi-solution of inverting, the mathematics that the assurance inversion stage adopts, the method for physics have geological Significance.
The present invention is achieved by the following technical solutions:
A kind of inverting modeling method based on faults control said method comprising the steps of:
(1) gridding seismic section: the coordinate that obtains seismic horizon and line of rent by structure elucidation from original earthquake data, then seismic section is counted seismic section after grid division obtains gridding by number of samples and road, the ordinate of grid is number of samples, and horizontal ordinate is the road number of seismic trace; On the seismic section after the described gridding, the corresponding grid data of a sampling point on each seismic trace;
(2) definition control line: according to the coboundary of the zone of interest in the seismic section after the described gridding and the control line of two levels of lower boundary definition, what be positioned at the top is upper control line, control line under being below being positioned at; Described zone of interest refers to the reservoir of oily, by two zones that seismic horizon forms in the seismic section after the gridding, the seismic horizon that is positioned at the top is called as the coboundary of this zone of interest, and the seismic horizon that is positioned at the below is called as the lower boundary of this zone of interest; Described upper control line is positioned at the top of described coboundary, and described lower control line is positioned at the below of described lower boundary;
(3) ask for the intersection point of line of rent and control line: from the seismic section after the described gridding, find each line of rent respectively with the intersection point of upper control line and lower control line;
(4) do perpendicular line: in zone of interest, cross each line of rent and do the boost line vertical with upper control line with the intersection point of upper control line, cross each line of rent and do the boost line vertical with lower control line with the intersection point of lower control line;
(5) definition grid special processing zone: in zone of interest, the zone that is surrounded by boost line, line of rent and seismic horizon and be exactly grid special processing zone by the zone that line of rent and seismic horizon surround; In zone of interest, the zone beyond the described grid special processing zone just is non-grid special processing zone;
(6) grid property in definition grid special processing zone: in the grid special processing zone that is obtained by step (5), from line of rent begin to do one group parallel with line of rent and the calibration line of width interval arranged each other; Described calibration line is the identical curve of straight line or shape with described line of rent;
(7) utilize the wave impedance curve in the log data that the interpolation calculation result that interpolation calculation obtains grid special processing zone is carried out in described grid special processing zone;
(8) utilize the wave impedance curve in the log data that the interpolation calculation result that vertical interpolation calculation obtains non-grid special processing zone is carried out in described non-grid special processing zone;
The interpolation calculation result in the interpolation calculation result in described grid special processing zone and non-grid special processing zone has consisted of inverse model.
Described step (2) is achieved in that the seismic section after the computer screen coordinate system with standard defines described gridding, and namely the upper left corner of the seismic section after the gridding is as initial point; The value of the ordinate of upper control line is less than the minimum value of the value of the ordinate of described coboundary, and the value of the ordinate of lower control line is greater than the maximal value of the value of the ordinate of lower boundary; The seismic horizon data recording the position of each sampling point on the seismic section after the gridding on this seismic horizon, earthquake layer bit data is done traversal can obtain this seismic horizon with respect to maximal value and the minimum value of computer screen coordinate system; When the seismic section after the described gridding has a plurality of zone of interest, set up successively upper control line and the lower control line of each zone of interest.
In the described step (4), if the point that line of rent and control line, then utilize the line of rent least significant end without natural intersection point therewith any point on the line of rent do connecting line, and this connecting line is prolonged until it is crossing with control line; Described any point is better close to more from the least significant end of line of rent.
Described step (7) be achieved in that will be corresponding with calibration line with the calibration line direction in the grid special processing zone grid data sort out and to calculate, namely utilize computerese will the grid data corresponding with calibration line to classify as the array of different length, adopt interpolation algorithms to carry out the interpolation calculation result that interpolation calculation obtains grid special processing zone to these arrays.
Described step (8) is achieved in that perpendicular to the direction of computer screen horizontal ordinate vertical grid data is sorted out and calculates, namely utilize computerese that vertical grid data is classified as different arrays, adopt interpolation algorithm to carry out the interpolation calculation result that interpolation calculation obtains non-grid special processing zone to these arrays.
Described interpolation algorithm adopts association's Kriging method or Based on Interpolating Spline or anti-distance weighted algorithm or nature interpolation algorithm at zero point.
Compare with existing inverting modeling technique, the invention has the beneficial effects as follows:
(1) form of the inventive method is simple, utilized innovatively grid data the complicated geological structure that is cut by line of rent and seismic horizon to be made the definition of inclination angle tendency, and these information recording /s can be become file, make it when carrying out modeling, can embody the continuity that fault information can well guarantee again zone of interest.
(2) the inventive method utilization is stored in geological data and the decryption (comprising line of rent and seismic horizon) of the SEGY standard format in the computing machine, in conjunction with Geometry method current seismic section being done gridding processes, only need to set up preferably by eight steps the model of a certain seismic section, and so that model can significantly embody fault information.The method is beneficial to the realization of computerized algorithm more, and can guarantee to reach good actual effect according to the seismic inversion that this model carries out.
Description of drawings
Fig. 1-the 1st the present invention is based on the synoptic diagram of line of rent, seismic horizon and the well of the original earthquake data among the inverting modeling method embodiment of faults control.
Fig. 1-2 the present invention is based on the synoptic diagram that among the inverting modeling method embodiment of faults control original earthquake data is converted to grid data.
Fig. 2 is the synoptic diagram that the present invention is based on the upper and lower control line among the inverting modeling method embodiment of faults control.
Fig. 3 is the synoptic diagram of doing perpendicular line that the present invention is based among the inverting modeling method embodiment of faults control.
Fig. 4 is the synoptic diagram that the present invention is based on the definition special processing zone among the inverting modeling method embodiment of faults control.
Fig. 5 is the synoptic diagram that the present invention is based on the definition special processing REGION INTERPOLATION mode among the inverting modeling method embodiment of faults control.
Fig. 6 is the step block diagram that the present invention is based on the inverting modeling method of faults control.
Fig. 7-the 1st, the as a result figure after the interpolation processing is carried out in the special processing zone.
Fig. 7-the 2nd, the as a result figure after the interpolation processing is carried out in non-special processing zone.
Fig. 8-the 1st, the model that utilizes prior art to obtain.
Fig. 8-the 2nd, the model that utilizes the inventive method to obtain.
Embodiment
Below in conjunction with accompanying drawing the present invention is described in further detail:
A kind of inverting modeling method based on faults control, its utilization are stored in geological data and the decryption (comprising line of rent and seismic horizon) of the SEGY standard format in the computing machine and process, and as shown in Figure 6, said method comprising the steps of:
(1) gridding seismic section
(" structure elucidation " Fig. 6 is exactly in order to obtain the coordinate of seismic horizon and line of rent at geological data to obtain the coordinate of seismic horizon and line of rent by structure elucidation from original earthquake data.As long as have geological data just passable during general gridding, but for the present invention, the coordinate of seismic horizon or line of rent has been arranged and be recorded in the gridding information, grid data has just had the information of seismic horizon and line of rent like this), then seismic section is divided into some grids by number of samples and road number, the ordinate of grid is number of samples, and horizontal ordinate is number, so that how much of late time data calculate.In general the sampling point with each seismic trace defines a grid data on a certain seismic section.For computerese, a grid data that is defined namely is an one-dimension array.Grid data is the internal data of storing with lattice structure in the computing machine, is the direct product of scan-type digitizer, is applicable to screen display and line printing output.In grid data, research range is divided into graticule mesh matrix of uniform size.The information of storage can be the point, line, surface entity, also can be the pointer that points to this relevant attribute in unit.Graticule mesh is less, and precision is higher, but storage capability is larger.Because graticule mesh is regular arrangement, the coordinate position of historical facts or anecdotes body can lie in the memory address of graticule mesh.Grid data is convenient to data processing, regional complex analysis and evaluation.Compare with vector data, its Software for Design is simpler, and shortcoming is that memory data output is large, and particularly sparse spatial data will be wasted many storage unit.Grid data is applicable to digital terrain model, the storage of the information such as remote sensing images." grid data that is defined " generally can be defined as one-dimension array, such as this one-dimension array 10 spaces arranged, and ordinate, the 2nd space preservation horizontal ordinate are preserved in the 1st space, and other space can save colors, the parameter of attribute.Use one-dimension array to be because it realizes that with computerized algorithm upper meeting is more convenient.
Fig. 1-1 is the line of rent in the original seismic section (tomography is geologic saying, the line segment of definition when line of rent refers to explain), seismic horizon and well synoptic diagram, and Fig. 1-2 is the seismic cross-section after the gridding.Well among Fig. 1 represents on this seismic section Liang Koujing is arranged, and is respectively well 1 and well 2, rather than refers to that seismic section is between Liang Koujing.The back all be to carry out on the basis of Fig. 1-2 in steps, for the work of each step below clear the demonstration, in Fig. 2, Fig. 3 and Fig. 4, removed those grids.
(2) define upper and lower control line:
For two parallel control lines in coboundary and lower boundary definition of zone of interest, when in the earthquake section a plurality of zone of interest being arranged, to set up successively the upper and lower control line of each zone of interest.Described zone of interest refers generally to the reservoir of oily, refer in the present invention the zone that formed by two seismic horizons, the zone that forms such as the seismic horizon h-1 among Fig. 1-1 and seismic horizon h-2 can be called zone of interest, wherein h-1 is called as the coboundary of this zone of interest, and h-2 is called as the lower boundary of this zone of interest.Fig. 2 is an example, and wherein, what horizontal curve h-1, h-2 and h-3 represented is 3 seismic horizons, and oblique curve f1, f3 and f2 are line of rent, represent 3 tomographies.
As shown in Figure 2, define this earthquake section with the computer screen coordinate system of standard, the computer screen coordinate system of described standard just refers to that the upper left corner is initial point, just is take the upper left corner of this picture in Fig. 2 on initial point (longitudinal and transverse coordinate all be on the occasion of).Upper control line L1 requires the minimum value less than the coboundary, and lower control line L2 requires the maximal value greater than lower boundary.The data recording of seismic horizon each sampling point on this seismic horizon in the position of seismic section (being on the corresponding grid one in length and breadth coordinate points), therefore seismic horizon is done traversal can obtain its seismic horizon with respect to maximal value and the minimum value of computer screen coordinate system, has namely obtained the minimum value of coboundary and the maximal value of lower boundary.When having a plurality of zone of interest, can set up successively section upper control line and the lower control line of each zone of interest.Upper control line L1 in the example shown in Figure 2 and lower control line L2 refer to the upper and lower control line of the zone of interest that is made of h-1 and h-3, and this zone of interest is actually all zone of interest that h-1, h-2, h-3 are controlled and has all comprised.If only need pay close attention to the zone of interest that h-1 and h-2 consist of, the maximal value place that descends so control line L2 to be arranged on greater than h-2 gets final product.
(3) ask for the intersection point of line of rent and control line:
As shown in Figure 3, find line of rent f1, f2, f3 to arrive the intersection point of the extended line of upper and lower control line L1 and L2, if line of rent and control line are without natural intersection point, the point that can utilize the so the line of rent least significant end therewith any point on the line of rent (better close to more from end) is done connecting line, and this connecting line is prolonged until intersect with upper and lower control line.
(4) do perpendicular line:
Cross line of rent and do the boost line vertical with control line with the intersection point of control line, as shown in Figure 4, be the line segment J1 vertical with L1 at f1 with the L1 intersection point, be the line segment J2 vertical with L2 at f1 with the L2 intersection point, J1 and J2 are boost line.Come the part line of rent of modeling to do perpendicular line to being necessary in the zone of interest according to line of rent fracture trend.The purpose that adds perpendicular line is to cut grid special processing zone smallerly, dealing with precision can be better, among Fig. 4, the zone that line of rent f2 and f3, h1, h2 control is the zone that is cut out by seismic horizon and line of rent, and the zone is also little, has therefore just no longer done perpendicular line.
(5) definition grid special processing zone:
As shown in Figure 4, utilize the up and down seismic horizon h-1 of boost line J1 and J2, line of rent f1 and zone of interest, regional C1, C2 that h-2 cuts out, C3, C4 and all be called as grid special processing zone by the zone that line of rent f2 and f3, h-1, h-2 cut out.Add that perpendicular line can cut smallerly to special area, dealing with precision can be better, and the zone itself that line of rent f2 and f3, h-1, h-2 cut out is also little, has therefore just no longer done perpendicular line.
The zone of only having zone that boost line, line of rent and seismic horizon cut out and seismic horizon and line of rent to cut out is grid special processing zone (being called processing region in Fig. 6 for short), and other zones are non-grid special processing zone (being called non-processing region in Fig. 4 for short).All dark, grayish zones all are grid special processing zones among Fig. 4, have only marked four of C1-C4 wherein among the figure, and the numbering in other special processing zone does not provide, and deep mixed sample is clear for Image Display among Fig. 4, does not have other particular meaning.
(6) grid property in definition special processing zone:
Cutting out grid special processing zone, be that (line of rent herein is the part of whole piece line of rent along line of rent in the zones such as C1, C2, C3, C4, be arranged in exactly the part in this zone, shown in the thick oblique line bar of rightmost of C1 piece among Fig. 5, its direction is shown in the line of rent trend of the mark of the four-headed arrow among Fig. 5) (width interval is take ready-portioned grid data as a unit to do the one fixed width interval, can between 1～10 unit, adjust, try concrete situation and decide) parallel lines be calibration line.Fig. 5 is take the C1 piece as example, and in general the width interval between all calibration lines is the same or different, if all be the same, the width interval of the last item may be smaller, because the zone not necessarily just in time comprises so much grid cell.
Described line of rent may be straight line or curve, if line of rent is straight line, then calibration line also all is straight line, and these straight lines are parallel to each other.If curve, then calibration line is exactly and the duplicate curve of described line of rent, and when the calibration line adjacent with line of rent and the Trendline of described tomography were etc., adjacent two calibration lines also were when waiting; The distance that refers in the horizontal direction the sampling point of sampling point to another line on the every line during described wait all is the same, i.e. described width interval.All be called parallel with two kinds among the present invention.
(7) interpolation method in special processing zone:
Utilize wave impedance curve in the log data to carrying out interpolation calculation in the special processing that the cuts out zone (being C1, C2, the zones such as C3, C4 among Fig. 4), the interpolation method in special processing zone can be with methods such as association's Kriging method, Based on Interpolating Spline, anti-distance weighted algorithm, nature interpolation algorithms at zero point.Describe as an example of anti-distance weighted algorithm example, anti-distance weighted algorithm (a kind of mathematical method) is inverse distance power gridding method, is a weighted mean method of interpolation.The degree parameter is being controlled weight coefficient and how to be descended along with the increase of leaving a graticule mesh node distance.For a larger degree, the given higher weight share of nearer data point, for a less degree, weight ratio is distributed to each data point more equably.It is proportional that the weights that give a particular data point when calculating a graticule mesh node and this node from the node to the observation station of specifying degree are endowed inverse distance.When calculating a graticule mesh node, the weight of dispensing is a mark, and the summation of all weights equals 1.0.When an observation station overlaps with a graticule mesh node, this observation station be given one actual be 1.0 weight, all other observation stations are given one and are almost 0.0 weight.In other words, this node is assigned to the value consistent with observation station, an accurate interpolation that Here it is.
Utilize the wave impedance curve in the log data specific as follows to the process of carrying out interpolation calculation in the special processing zone that cuts out: well medium wave impedance curve is positioned at a place, (is on the calibration line with the well middle impedance curve mapping of a direction to the c place parallel with line of rent.), well middle impedance curve such as a place has M1 sample value, c place curve contains M2 sample value, and M1 sample value with a place when calculating is evenly distributed on M2 the sample value at c place, namely obtained the impedance curve parallel with line of rent at c place by the well middle impedance curve at a place.Then M2 resistance value with the c place is evenly distributed on (being said calibration line) on its next bar line parallel with line of rent.Through after the mapping of above-mentioned well middle impedance curve values, obtained a, b, c, sample value on each bar straight line (being said calibration line) parallel with line of rent in the d enclosed region (as shown in Figure 5), then calculate each bar straight line and grid distance (straight line and grid distance refer to number of samples), the grid resistance value is by the acquisition (be exactly the impedance mean value that is arranged on the straight line of this grid with regard to the resistance value that refers to a lattice among Fig. 5) of averaging of the resistance value on the straight line that passes this grid.Thereby obtain the resistance value of grid place in the whole enclosed region (the grid place refers to a grid among Fig. 5).
It should be noted that using these methods in the special processing zone is not with the direction perpendicular to the computer screen coordinate, but with the calibration line direction in special processing zone (" correspondence " refers to all grid datas in the special processing zone) grid data classification of correspondence is calculated, namely utilize computerese that the data of correspondence are classified as the array of different length (because the line that special area is done is not necessarily isometric.), adopt above interpolation algorithm to carry out interpolation calculation to these arrays.The result of calculation in special processing zone (has only provided the result of calculation in part special processing zone) among Fig. 7-1 shown in Fig. 7-1, as we can see from the figure according to after the fault trend interpolation, model moves towards can be along the trend of seismic horizon.
(8) interpolation method in non-special processing zone:
Utilize the wave impedance curve in the log data to do vertical interpolation processing in non-special processing zone (being all white portions among Fig. 4), the interpolation method in non-special processing zone can be with assisting the methods such as Kriging method, Based on Interpolating Spline, anti-distance weighted algorithm, nature interpolation algorithm at zero point, that is to say that special area and non-special area all use identical wave impedance curve and identical interpolation method, just direction is different.In general, special area and non-special area must be selected the same interpolation algorithm, and to keep the consistance of whole numerical values recited, because different interpolation methods, the result after the interpolation has some differences.It should be noted that the direction that must be perpendicular to the computer screen horizontal ordinate when in non-special processing zone, using these methods, the namely what is often called directions X.The grid data classification of vertical direction is calculated, the result of calculation in non-special processing zone (has only provided the result of calculation in the non-special processing of part zone) among Fig. 7-1 shown in Fig. 7-2, not regional according to the non-special processing of line of rent trend interpolation as we can see from the figure, the data value of model is basically identical, utilize computerese that the data of correspondence are classified as the array of different length, adopt above interpolation algorithms to carry out interpolation calculation to these arrays.
These two-part data have all been calculated complete when the regional and non-special processing of special processing zone, namely finished modeling, utilize model that the inventive method generates shown in Fig. 8-2, light can find out that this model presses close to seismic horizon and fault strike trend from figure, and utilize model that prior art generates shown in Fig. 8-1, light can be found out this model level from figure, does not press close to seismic horizon and fault strike trend.
The inventive method is not only applicable to the two-dimension earthquake data and is applicable to too 3D seismic data, in fact can adopt the mode of two-dimension earthquake data to finish for 3-d seismic data set, because 3-d seismic data set can be divided into some line data, the line data just are equivalent to a two-dimension earthquake section, after finishing processing to these line data with method of the present invention successively, namely be equivalent to finish whole 3-D data volume.Account form is just the same with two dimension, just description of drawings no longer.
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 be not limited only to the described method of the above-mentioned embodiment of the present invention, therefore previously described mode is just preferred, and does not have restrictive meaning.
Claims (6)
1. inverting modeling method based on faults control is characterized in that: said method comprising the steps of:
(1) gridding seismic section: the coordinate that obtains seismic horizon and line of rent by structure elucidation from original earthquake data, then seismic section is counted seismic section after grid division obtains gridding by number of samples and road, the ordinate of grid is number of samples, and horizontal ordinate is the road number of seismic trace; On the seismic section after the described gridding, the corresponding grid data of a sampling point on each seismic trace;
(2) definition control line: according to the coboundary of the zone of interest in the seismic section after the described gridding and the control line of two levels of lower boundary definition, what be positioned at the top is upper control line, control line under being below being positioned at; Described zone of interest refers to the reservoir of oily, by two zones that seismic horizon forms in the seismic section after the gridding, the seismic horizon that is positioned at the top is called as the coboundary of this zone of interest, and the seismic horizon that is positioned at the below is called as the lower boundary of this zone of interest; Described upper control line is positioned at the top of described coboundary, and described lower control line is positioned at the below of described lower boundary;
(3) ask for the intersection point of line of rent and control line: from the seismic section after the described gridding, find each line of rent respectively with the intersection point of upper control line and lower control line;
(4) do perpendicular line: in zone of interest, cross each line of rent and do the boost line vertical with upper control line with the intersection point of upper control line, cross each line of rent and do the boost line vertical with lower control line with the intersection point of lower control line;
(5) definition grid special processing zone: in zone of interest, the zone that is surrounded by boost line, line of rent and seismic horizon and be exactly grid special processing zone by the zone that line of rent and seismic horizon surround; In zone of interest, the zone beyond the described grid special processing zone just is non-grid special processing zone;
(6) grid property in definition grid special processing zone: in the grid special processing zone that is obtained by step (5), from line of rent begin to do one group parallel with line of rent and the calibration line of width interval arranged each other; Described calibration line is the identical curve of straight line or shape with described line of rent;
(7) utilize the wave impedance curve in the log data that the interpolation calculation result that interpolation calculation obtains grid special processing zone is carried out in described grid special processing zone;
(8) utilize the wave impedance curve in the log data that the interpolation calculation result that vertical interpolation calculation obtains non-grid special processing zone is carried out in described non-grid special processing zone;
The interpolation calculation result in the interpolation calculation result in described grid special processing zone and non-grid special processing zone has consisted of inverse model.
2. the inverting modeling method based on faults control according to claim 1, it is characterized in that: described step (2) is achieved in that the seismic section after the computer screen coordinate system with standard defines described gridding, and namely the upper left corner of the seismic section after the gridding is as initial point; The value of the ordinate of upper control line is less than the minimum value of the value of the ordinate of described coboundary, and the value of the ordinate of lower control line is greater than the maximal value of the value of the ordinate of lower boundary; The seismic horizon data recording the position of each sampling point on the seismic section after the gridding on this seismic horizon, earthquake layer bit data is done traversal can obtain this seismic horizon with respect to maximal value and the minimum value of computer screen coordinate system; When the seismic section after the described gridding has a plurality of zone of interest, set up successively upper control line and the lower control of each zone of interest.
3. the inverting modeling method based on faults control according to claim 1, it is characterized in that: in the described step (4), if line of rent and control line are without natural intersection point, the point that then the utilizes the line of rent least significant end therewith any point on the line of rent is done connecting line, and this connecting line is prolonged until it is crossing with control line; Described any point is better close to more from the least significant end of line of rent.
4. the inverting modeling method based on faults control according to claim 1, it is characterized in that: described step (7) be achieved in that will be corresponding with calibration line with the calibration line direction in the grid special processing zone grid data sort out and to calculate, namely utilize computerese will the grid data corresponding with calibration line to classify as the array of different length, adopt interpolation algorithms to carry out the interpolation calculation result that interpolation calculation obtains grid special processing zone to these arrays.
5. the inverting modeling method based on faults control according to claim 1, it is characterized in that: described step (8) is achieved in that perpendicular to the direction of computer screen horizontal ordinate vertical grid data is sorted out and calculates, namely utilize computerese that vertical grid data is classified as different arrays, adopt interpolation algorithm to carry out the interpolation calculation result that interpolation calculation obtains non-grid special processing zone to these arrays.
6. it is characterized in that according to claim 4 or 5 described inverting modeling methods based on faults control: described interpolation algorithm adopts association's Kriging method or Based on Interpolating Spline or anti-distance weighted algorithm or nature interpolation algorithm at zero point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210110254.XA CN103376463B (en) | 2012-04-13 | 2012-04-13 | A kind of Inverse modeling method based on faults control |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210110254.XA CN103376463B (en) | 2012-04-13 | 2012-04-13 | A kind of Inverse modeling method based on faults control |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103376463A true CN103376463A (en) | 2013-10-30 |
| CN103376463B CN103376463B (en) | 2016-02-03 |
Family
ID=49461860
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210110254.XA Active CN103376463B (en) | 2012-04-13 | 2012-04-13 | A kind of Inverse modeling method based on faults control |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103376463B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105137485A (en) * | 2015-08-31 | 2015-12-09 | 中国石油天然气股份有限公司 | Automatic extraction and conversion method and device for interpretation information in seismic section bitmap |
| CN103914620B (en) * | 2014-03-31 | 2017-02-15 | 成都理工大学 | Method for computing distribution of opening spaces of fractures of fault rupture zones |
| CN107367762A (en) * | 2017-06-27 | 2017-11-21 | 中国石油天然气股份有限公司 | A kind of method and device for determining reservoir parameter |
| CN110879681A (en) * | 2018-09-06 | 2020-03-13 | 中国石油天然气股份有限公司 | Method and device for displaying fault fracture strength |
| CN110967737A (en) * | 2018-09-30 | 2020-04-07 | 中国石油化工股份有限公司 | Initial model construction method for construction constraint |
| CN111812729A (en) * | 2020-05-27 | 2020-10-23 | 中煤科工集团西安研究院有限公司 | Multi-dimensional geological background model construction method for optimizing electrical prospecting inversion |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030132934A1 (en) * | 2001-12-12 | 2003-07-17 | Technoguide As | Three dimensional geological model construction |
| US20030216897A1 (en) * | 2002-05-17 | 2003-11-20 | Schlumberger Technology Corporation | Modeling geologic objects in faulted formations |
| CN102385067A (en) * | 2011-08-03 | 2012-03-21 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Drawing method for isoline containing reverse fault |
-
2012
- 2012-04-13 CN CN201210110254.XA patent/CN103376463B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030132934A1 (en) * | 2001-12-12 | 2003-07-17 | Technoguide As | Three dimensional geological model construction |
| US20030216897A1 (en) * | 2002-05-17 | 2003-11-20 | Schlumberger Technology Corporation | Modeling geologic objects in faulted formations |
| CN102385067A (en) * | 2011-08-03 | 2012-03-21 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Drawing method for isoline containing reverse fault |
Non-Patent Citations (3)
| Title |
|---|
| 刘少兵: "《测井层析储层成像系统》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 魏嘉: "《地质建模技术》", 《勘探地球物理进展》 * |
| 魏嘉等: "《三维地质构造建模技术研究》", 《石油物探》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103914620B (en) * | 2014-03-31 | 2017-02-15 | 成都理工大学 | Method for computing distribution of opening spaces of fractures of fault rupture zones |
| CN105137485A (en) * | 2015-08-31 | 2015-12-09 | 中国石油天然气股份有限公司 | Automatic extraction and conversion method and device for interpretation information in seismic section bitmap |
| CN105137485B (en) * | 2015-08-31 | 2017-12-05 | 中国石油天然气股份有限公司 | The automatic Picking method for transformation and device of information are explained in a kind of seismic profile bitmap |
| CN107367762A (en) * | 2017-06-27 | 2017-11-21 | 中国石油天然气股份有限公司 | A kind of method and device for determining reservoir parameter |
| CN107367762B (en) * | 2017-06-27 | 2019-06-11 | 中国石油天然气股份有限公司 | A kind of method and device of determining reservoir parameter |
| CN110879681A (en) * | 2018-09-06 | 2020-03-13 | 中国石油天然气股份有限公司 | Method and device for displaying fault fracture strength |
| CN110879681B (en) * | 2018-09-06 | 2021-08-03 | 中国石油天然气股份有限公司 | Method and device for displaying fault fracture strength |
| CN110967737A (en) * | 2018-09-30 | 2020-04-07 | 中国石油化工股份有限公司 | Initial model construction method for construction constraint |
| CN110967737B (en) * | 2018-09-30 | 2021-12-24 | 中国石油化工股份有限公司 | Initial model construction method for construction constraint |
| CN111812729A (en) * | 2020-05-27 | 2020-10-23 | 中煤科工集团西安研究院有限公司 | Multi-dimensional geological background model construction method for optimizing electrical prospecting inversion |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103376463B (en) | 2016-02-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8655632B2 (en) | Gridless geological modeling | |
| US10838108B2 (en) | Method and system for interpolating discontinuous functions in a subsurface model | |
| Sech et al. | Three-dimensional modeling of a shoreface-shelf parasequence reservoir analog: Part 1. Surface-based modeling to capture high-resolution facies architecture | |
| EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
| US7340385B2 (en) | Method and program storage device for generating grids representing the architecture of fluvial reservoirs | |
| CN103376463B (en) | A kind of Inverse modeling method based on faults control | |
| Pyrcz et al. | A library of training images for fluvial and deepwater reservoirs and associated code | |
| CN103699751A (en) | Sand body reservoir architecture modeling method and system based on space vectors | |
| AU2014306493B2 (en) | Determining reserve estimates for a reservoir | |
| WO2018222331A1 (en) | Constructing structural models of the subsurface | |
| GB2474740A (en) | Gridless geological modeling of a structural framework | |
| CN103279986A (en) | Three-dimensional horizontal geological section map manufacturing method and applications thereof | |
| CN104459768A (en) | Three-dimensional space target geologic body tracking method based on visualization | |
| CN112150582B (en) | Multi-modal data-oriented geological profile approximate expression method | |
| CN102867332A (en) | Multi-level subdivided mesh surface fitting method based on complicated boundary constraint | |
| Keogh et al. | Data capture for multiscale modelling of the Lourinha Formation, Lusitanian Basin, Portugal: an outcrop analogue for the Statfjord Group, Norwegian North Sea | |
| Bigi et al. | Discrete fracture network of the Latemar carbonate platform | |
| CN106846481A (en) | Geological profile generation method | |
| CN115880455A (en) | Three-dimensional intelligent interpolation method based on deep learning | |
| Sahin et al. | Stochastic modeling of porosity distribution in a multi-zonal carbonate reservoir | |
| CN111562629A (en) | Saturation determination method and device based on equivalent pore section index | |
| Li et al. | Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China | |
| CN106646608A (en) | Layer model establishing method for observation system parameter argumentation | |
| CN115469361B (en) | Clastic rock stratum three-dimensional geological modeling method | |
| MAHJOUR | Using Data Management Approaches to Improve Decision Analysis in Petroleum Field Development under Uncertainty |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |