CN103913773A - Method for achieving encrypting microlog by means of residual error abnormity during three-dimensional exploration - Google Patents
Method for achieving encrypting microlog by means of residual error abnormity during three-dimensional exploration Download PDFInfo
- Publication number
- CN103913773A CN103913773A CN201310351921.8A CN201310351921A CN103913773A CN 103913773 A CN103913773 A CN 103913773A CN 201310351921 A CN201310351921 A CN 201310351921A CN 103913773 A CN103913773 A CN 103913773A
- Authority
- CN
- China
- Prior art keywords
- velocity
- point
- microlog
- velocity bed
- elevation
- 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 relates to a method for achieving encrypting microlog by means of a residual error abnormity during three-dimensional exploration. The method comprises the steps that microlog points are evenly distributed through diluting, the low-velocity layer, the decelerating layer, the velocity of the high-speed layer, and the thickness of the high-speed layer of each microlog point are calculated, the measured thickness of each triggered high-speed layer and the measured elevation of the top interface of each triggered high-speed layer are obtained, the thickness of each triggered high-speed layer is obtained, an abnormal area is found out, the microlog points in the abnormal area are encrypted till the abnormal area is eliminated, and reasonable arrangement of microlog points is achieved. By the adoption of the method for achieving encrypting microlog by means of the residual error abnormity during three-dimensional exploration, encrypting microlog can be effectively conducted in an abnormal area, caused by the undulating surface or underground geological conditions, on the surface layer, microlog points can be effectively arranged in the normal area in a diluting mode, as a result, the workload is lightened, the construction cost is saved, and the working efficiency is improved.
Description
Technical field
The invention belongs to petroleum gas seismic prospecting technique for investigation field, relate to one and be applied in the three dimensional seismic data collection of oil seismic exploration field, utilize residual error extremely to encrypt the method for micro logging.
Background technology
In seismic prospecting, the precision of surficial geology energy structure survey is most important, can it had both been related to provide reliable foundation for the selection of earthquake stimulation parameter, can be also related to low for top layer, the static correction of reduction of speed band provides accurate data, do the investigation of surface infrastructure model well, can determine best mode of excitation and receive mode, optimize Hydrographic General Line, lay the foundation for indoor and field static correcting value calculate.
In recent years, in high precision seismic exploration application, the investigation precision of surface structure is related to the success or failure of exploration effects especially, the method that the control points layout method adopting at present adopts Large-Area-Uniform to lay, be unfavorable for finding the catastrophe point of surface structure, surface structure model and the actual error therefore investigated out are larger.
Summary of the invention
The object of this invention is to provide in a kind of three-dimensional exploration and utilize residual error method extremely to encrypt the method for micro logging, can not be in the problem of the three-dimensional micro logging of relief surface Reasonable Arrangement to solve prior art.
For achieving the above object, in three-dimensional exploration of the present invention, utilize residual error method extremely to encrypt the step of method of micro logging as follows:
A, in three-dimensional work area to be onstructed, measure each shot point coordinate and earth's surface elevation H according to construction requirement
1in order to obtain the low velocity layer data of each shot point, first measuring the accuracy requirement of former work area, evenly to lay the density of micrometering well point be m*m, adopt mn*mn vacuate evenly to lay micrometering well point, n is greater than 1 number, and field carries out data acquisition, obtain the seismologic record of different depth of every mouthful of micrometering well point and coordinate, earth's surface elevation H
2data;
B, undertaken after primary wave contrast by the seismologic record to obtained, make the vertical T-X curve of t-h, according to the difference of t-h rate of curve, mark off different velocity layerings, obtain low velocity layer (LVL), reduction of speed layer, high-velocity bed speed and the high-velocity bed thickness D of each micrometering well point
1, and calculate high-velocity bed elevation of top surface H
3=earth's surface elevation H
2-high-velocity bed thickness D
1, wherein t refers to the primary wave time of the each acceptance point of micro logging, h refers to the degree of depth of each acceptance point;
C, use kriging interpolation method are to micrometering well point high-velocity bed thickness D
1, high-velocity bed elevation of top surface H
3data are carried out gridding, then carry out residual error with Yi Ce whole district shot point, obtain the each high-velocity bed thickness D exciting having measured in three-dimensional work area
2with high-velocity bed elevation of top surface H
4;
D, with the earth's surface elevation H of each shot point
1high-velocity bed elevation of top surface H with step c calculating
4carry out difference calculating, obtain each shot point high-velocity bed thickness D
3;
E, each shot point high-velocity bed thickness D that step c is obtained
2each shot point high-velocity bed thickness D with steps d acquisition
3carry out difference, and the threshold value of this difference and setting is compared, if it is abnormal area that difference exceedes the region of threshold value, in abnormal area, micrometering well point is encrypted at the position of difference maximum, quantity is determined as the case may be, then repeating step b-e, until eliminate abnormal area, realizes the Reasonable Arrangement of micro logging.
In described step a, the density of shot point is greater than the accuracy requirement of former work area and evenly lays the density of micrometering well point.
The threshold value of setting in described step e is 5m.
In described step e, difference becomes figure to show by isoline.
In three-dimensional exploration of the present invention, utilize residual error method extremely to encrypt the method for micro logging, having solved prior art can not be in the problem of the three-dimensional micro logging of relief surface Reasonable Arrangement, the method that in the past just strengthened the even layout density in micrometering well point for relief surface has increased construction cost and workload, micro logging is effectively encrypted in the variation abnormality region, top layer that this method not only can be intuitively causes relief surface or subsurface geology condition, and to without the abnormal area effectively laying of vacuate micrometering well point again, to reduce workload, save construction cost, increase work efficiency.
Accompanying drawing explanation
Fig. 1 is the method flow diagram of the embodiment of the present invention;
Fig. 2 is the shot point of the embodiment of the present invention, the micrometering well point location drawing;
Fig. 3 is the data recording figure of the single port micro logging of the embodiment of the present invention;
Fig. 4 is the vertical T-X curve figure of single port micro logging t-h of the embodiment of the present invention;
Fig. 5 be the embodiment of the present invention to the figure before the abnormal encryption of high-velocity bed;
Fig. 6 be the embodiment of the present invention to the figure after the abnormal encryption of high-velocity bed;
Fig. 7 is high-velocity bed thickness situation map before the micro logging of the embodiment of the present invention is encrypted;
Fig. 8 is high-velocity bed thickness situation map after the micro logging of the embodiment of the present invention is encrypted.
Embodiment
Utilize residual error method extremely encrypt micro logging method process flow diagram as shown in Figure 1, concrete steps comprise:
A, certain high-precision three-dimensional work area construction area 320km
2, measure each shot point coordinate and earth's surface elevation H according to construction requirement
1as table 1, in order to obtain the low velocity layer data of each shot point, first measuring the accuracy requirement of former work area, evenly to lay the density of micrometering well point be 1km*1km, lay altogether 285 of micrometering well points, density is far smaller than the density of shot point, as shown in Figure 2, adopts 2km*2km vacuate evenly to lay 80 mouthfuls of micrometering well points micro logging, and field carries out data acquisition, obtain every mouthful of micrometering well point different depth seismologic record as shown in Figure 3 with coordinate, earth's surface elevation H
2data are as shown in table 2;
Table 1 part shot point coordinate, earth's surface elevation H
1
Table 2 part micro logging coordinate, earth's surface elevation H
2
B, undertaken after primary wave contrast by the record to obtained, can make the vertical T-X curve of t-h as Fig. 4, according to the difference of t-h rate of curve, mark off different velocity layerings, can obtain low velocity layer (LVL), reduction of speed layer, high-velocity bed speed and the high-velocity bed thickness D of each micrometering well point
1as table 3, and calculate high-velocity bed elevation of top surface H
3=earth's surface elevation H
2-high-velocity bed thickness D
1, wherein t refers to the primary wave time of the each acceptance point of micro logging, h refers to the degree of depth of each acceptance point; (particular content of micro logging data acquisition and Interpretation On The Results elaborates on 120 pages of textbook " PRINCIPLE OF SEISMIC PROSPECTING WITH " first volumes);
Table 3 part micro logging interpretation results
C, use surfer software, this software has gridding and residual error function, selects kriging interpolation method, respectively to micrometering well point high-velocity bed thickness D
1, high-velocity bed elevation of top surface H
3data are carried out gridding, then carry out residual error with Yi Ce whole district shot point, obtain respectively the each high-velocity bed thickness D exciting having measured in three-dimensional work area
2with high-velocity bed elevation of top surface H
4as table 4;
Table 4 part shot point high-velocity bed thickness and high-velocity bed elevation of top surface data
| Shot point pile No. | High-velocity bed thickness D 2 | High-velocity bed elevation of top surface H 4 |
| 12611171 | 5.3 | 157.7 |
| 12611172 | 5.6 | 154.0 |
| 12611173 | 5.9 | 150.3 |
| 12611174 | 6.3 | 148.1 |
| 12611175 | 6.8 | 147.6 |
| 12611176 | 7.2 | 150.9 |
| 12611177 | 7.6 | 155.8 |
| 12611178 | 8.0 | 154.8 |
| 12611179 | 8.1 | 158.6 |
| 12611180 | 8.3 | 165.0 |
| 12611181 | 8.3 | 152.1 |
| 12611182 | 8.3 | 147.7 |
| 12611183 | 8.2 | 148.4 |
| 12611184 | 8.1 | 150.4 |
| 12611185 | 7.8 | 153.4 |
| 12611186 | 7.5 | 155.7 |
| 12611187 | 7.3 | 154.7 |
| 12611188 | 7.1 | 155.5 |
| 12611189 | 7.0 | 143.9 |
| 12611190 | 7.0 | 139.6 |
| 12611191 | 7.2 | 140.4 |
| 12611192 | 7.3 | 137.1 |
| 12611193 | 7.5 | 129.3 |
| 12611194 | 7.6 | 129.4 |
D, with the earth's surface elevation H of each shot point
1high-velocity bed elevation of top surface H with step c calculating
4carry out difference calculating, obtain each shot point high-velocity bed thickness D
3as table 5;
The high-velocity bed thickness data that table 5 part shot point another kind method calculates
Each shot point high-velocity bed thickness D that e, step c obtain
2each shot point high-velocity bed thickness D with steps d acquisition
3carry out difference as table 6, become figure to show as Fig. 5 by isoline, the threshold value 5m of difference and setting is compared, the setting of threshold value is set according to work area actual conditions and experience, what difference exceeded 5m is abnormal area, in abnormal area, the general 1-2 mouth in micrometering well point is encrypted at the position of difference maximum, if abnormal ranges is larger, also can add close several mouthfuls, quantity is determined as the case may be, then repeating step b-e, until eliminate abnormal area as Fig. 6, encrypts 125 mouthfuls of micrometering well points altogether.
The high-velocity bed thickness difference data that table 6 part shot point distinct methods obtains
By above step, reduce micro logging and lay 80 mouthfuls, saved construction cost, the surface-level model of acquisition can meticulous depiction subsurface picture, especially surface structure change large area had for pass point, investigation more accurately as Fig. 7, Fig. 8.
Claims (4)
1. in three-dimensional exploration, utilize residual error method extremely to encrypt a method for micro logging, it is characterized in that, the step of the method is as follows:
A, in three-dimensional work area to be onstructed, measure each shot point coordinate and earth's surface elevation H according to construction requirement
1in order to obtain the low velocity layer data of each shot point, first measuring the accuracy requirement of former work area, evenly to lay the density of micrometering well point be m*m, adopt mn*mn vacuate evenly to lay micrometering well point, n is greater than 1 number, and field carries out data acquisition, obtain the seismologic record of different depth of every mouthful of micrometering well point and coordinate, earth's surface elevation H
2data;
B, undertaken after primary wave contrast by the seismologic record to obtained, make the vertical T-X curve of t-h, according to the difference of t-h rate of curve, mark off different velocity layerings, obtain low velocity layer (LVL), reduction of speed layer, high-velocity bed speed and the high-velocity bed thickness D of each micrometering well point
1, and calculate high-velocity bed elevation of top surface H
3=earth's surface elevation H
2-high-velocity bed thickness D
1, wherein t refers to the primary wave time of the each acceptance point of micro logging, h refers to the degree of depth of each acceptance point;
C, use kriging interpolation method are to micrometering well point high-velocity bed thickness D
1, high-velocity bed elevation of top surface H
3data are carried out gridding, then carry out residual error with Yi Ce whole district shot point, obtain the each high-velocity bed thickness D exciting having measured in three-dimensional work area
2with high-velocity bed elevation of top surface H
4;
D, with the earth's surface elevation H of each shot point
1high-velocity bed elevation of top surface H with step c calculating
4carry out difference calculating, obtain each shot point high-velocity bed thickness D
3;
E, each shot point high-velocity bed thickness D that step c is obtained
2each shot point high-velocity bed thickness D with steps d acquisition
3carry out difference, and the threshold value of this difference and setting is compared, if it is abnormal area that difference exceedes the region of threshold value, in abnormal area, micrometering well point is encrypted at the position of difference maximum, quantity is determined as the case may be, then repeating step b-e, until eliminate abnormal area, realizes the Reasonable Arrangement of micro logging.
2. method according to claim 1, is characterized in that: in described step a, the density of shot point is greater than the accuracy requirement of former work area and evenly lays the density of micrometering well point.
3. method according to claim 1, is characterized in that: the threshold value of setting in described step e is 5m.
4. according to the method described in any one in claim 1-3, it is characterized in that: in described step e, difference becomes figure to show by isoline.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310351921.8A CN103913773B (en) | 2013-08-13 | 2013-08-13 | Three-dimensional exploration utilizes the abnormal method encrypting micro logging of residual error method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310351921.8A CN103913773B (en) | 2013-08-13 | 2013-08-13 | Three-dimensional exploration utilizes the abnormal method encrypting micro logging of residual error method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103913773A true CN103913773A (en) | 2014-07-09 |
| CN103913773B CN103913773B (en) | 2016-09-28 |
Family
ID=51039597
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310351921.8A Active CN103913773B (en) | 2013-08-13 | 2013-08-13 | Three-dimensional exploration utilizes the abnormal method encrypting micro logging of residual error method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103913773B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0801364A2 (en) * | 1996-04-12 | 1997-10-15 | Institut Francais Du Petrole | Method for generating a 3D meshing |
| US6106561A (en) * | 1997-06-23 | 2000-08-22 | Schlumberger Technology Corporation | Simulation gridding method and apparatus including a structured areal gridder adapted for use by a reservoir simulator |
| CN101046515A (en) * | 2007-04-22 | 2007-10-03 | 罗仁泽 | Method of raising seismic resolution with micro measuring well perpendicular to seismic profile and double well |
| CN101178439A (en) * | 2007-03-22 | 2008-05-14 | 孙素梅 | Accurate ground surface investigation method |
| CN102590864A (en) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | Near-surface modeling method using tomography inversion of two-step method |
-
2013
- 2013-08-13 CN CN201310351921.8A patent/CN103913773B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0801364A2 (en) * | 1996-04-12 | 1997-10-15 | Institut Francais Du Petrole | Method for generating a 3D meshing |
| US6106561A (en) * | 1997-06-23 | 2000-08-22 | Schlumberger Technology Corporation | Simulation gridding method and apparatus including a structured areal gridder adapted for use by a reservoir simulator |
| CN101178439A (en) * | 2007-03-22 | 2008-05-14 | 孙素梅 | Accurate ground surface investigation method |
| CN101046515A (en) * | 2007-04-22 | 2007-10-03 | 罗仁泽 | Method of raising seismic resolution with micro measuring well perpendicular to seismic profile and double well |
| CN102590864A (en) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | Near-surface modeling method using tomography inversion of two-step method |
Non-Patent Citations (1)
| Title |
|---|
| 朱艳保 等: "泌阳凹陷北部高精度三维地震采集技术", 《石油与化工设备》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103913773B (en) | 2016-09-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102967882B (en) | The modeling method of the interval velocity model on stratum | |
| CN105651676B (en) | A kind of reservoir heterogeneity characterizing method under horizontal well rule well pattern | |
| CN100429531C (en) | 3D electromagnetic fast inversion method of minimized target | |
| CN103472483B (en) | Velocity modeling method based on true earth's surface or floating datum | |
| CN105386756B (en) | A method of brittle formation porosity is calculated using dependent variable | |
| CN103576212B (en) | A kind of labyrinth borehole restraint three-dimensional density sequence inversion method | |
| CN104453834A (en) | Injection-production relation optimizing and adjusting method for well group | |
| CN104570102B (en) | Method for combining near-surface velocity model with middle-deep stratum velocity model | |
| CN108089227B (en) | Novel stratum pore pressure prediction method based on three-dimensional seismic data | |
| CN103064124B (en) | A kind of ratio approach correcting electromagnetic survey Considering Terrain Effect | |
| CN105631753A (en) | Attitude of stratum-based horizontal well oil reservoir profile modeling method | |
| CN106405664B (en) | A kind of magnetic anomaly normalizing pole method | |
| CN103293553B (en) | Upper and lower cable seismic acquisition data boundary element continuation bearing calibration at the bottom of a kind of Complex Sea | |
| CN108121009B (en) | Development late stage complex fault block velocity-variant depth mapping method | |
| Küperkoch et al. | Long‐term monitoring of induced seismicity at the Insheim geothermal site, Germany | |
| CN106324682B (en) | A kind of surface structure survey method applied to frozen soil layer area | |
| CN106768187A (en) | A kind of many tidal level station seaway topographic survey tidal level control methods and water level self recording instrument fixing device | |
| CN103954996A (en) | Device and method for determining fracture strike of strata fracture based on travel time method | |
| CN105626056A (en) | Method and device for determining three-dimensional stratum pore pressure of target area | |
| CN103487830B (en) | A kind of method compensating differential compaction in seismic amplitude generic attribute exploration | |
| CN105116447A (en) | A geological river direction discrimination method based on curvature-abnormal stripes | |
| CN104749630A (en) | Method for constructing microseism monitoring velocity model | |
| CN104297796A (en) | Short-refraction investigation method for undulating surface of terrain high-variability area in ultra-thick loess highland | |
| CN104237946A (en) | Single-layer reflected P-wave and reflection converted shear wave amplitude matching method based on well control | |
| CN104570070A (en) | Method and equipment for building two-dimensional near-surface geological model |
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 |