CN1158539C - Optimizing design method of 3D seismic observation system based on geologic geophysical model - Google Patents
Optimizing design method of 3D seismic observation system based on geologic geophysical model Download PDFInfo
- Publication number
- CN1158539C CN1158539C CNB011204451A CN01120445A CN1158539C CN 1158539 C CN1158539 C CN 1158539C CN B011204451 A CNB011204451 A CN B011204451A CN 01120445 A CN01120445 A CN 01120445A CN 1158539 C CN1158539 C CN 1158539C
- Authority
- CN
- China
- Prior art keywords
- earthquake
- seismic
- imaging results
- geophysical model
- geologic
- 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.)
- Expired - Fee Related
Links
Images
Abstract
The present invention relates to a method of optimization design of a three-dimensional earthquake observing method based on a geologic geophysical model. The method comprises the steps that (A) a priori geologic geophysical model is established; (B) two or more types of observation systems of three-dimensional earthquake data acquisition are designed; (C) data acquisition is completed by the combined forward modeling of three-dimensional wave equation numerical simulation of a virtual spectral method and three-dimensional earthquake geophysical model simulation of a method; (D) depth migration before stacking is carried out on acquired data on the basis of a common reflective point, and some imaging results of earthquake treatment are obtained; (E) the imaging qualities of the imaging results of the earthquake treatment are comprehensively compared and evaluated, and one optimal earthquake imaging result is determined through comparison. The method has the advantages of wide range of application, and good applied effect, and the method has operation significance of guiding earthquake prospecting actual production of oil and natural gas, and enhancing oil-gas prospecting successful rate.
Description
Technical field
The present invention relates to a kind of technical design optimization method of 3-d seismic exploration in petroleum gas and other mineral resources and the Geological Engineering, particularly a kind of optimizing design method of 3 D seismic observation system based on geologic geophysical model.
Background technology
Seismic prospecting three big links are to gather, handle and explain.Collection is the basis; The receiving system that activating system that epicenter excitation point position is formed on the ground and wave detector acceptance point position are formed must be by correct geometric distributions, and promptly the seismic observation system of You Huaing just can collect effectively complete subsurface geological structure information.Therefore, recording geometry plays an important role in seismic prospecting, and its design quality is directly connected to the success or failure of exploration task.Because the 3-D seismics workload is big, investment is high, three-dimensional work area often is arranged in the complex structure area again, so the stereo observing system optimal design is important doubly just.
So far, both at home and abroad mostly the stereo observing system design software that widely adopts is based on horizontal layer geologic model and common midpoint (CMP) superposition principle, and they are not suitable for the 3-D seismics technical design on complex structure, high spud angle stratum.Under the complex structure situation, use this class design proposal, not only can not obtain high-quality seismic imaging, even the seismic imaging result that may make mistake.
Summary of the invention
Therefore, purpose of the present invention makes provides a kind of optimizing design method of 3 D seismic observation system based on geologic geophysical model.
This method is based on the three-dimensional geological and the geophysical model of exploratory area priori, adopt numerical simulation and physical simulation to unite positive artistic skills art, according to common reflection point (CRP) principle, be imaged as evaluation criterion with the earthquake of geologic objective high-quality, the image quality of three-dimensional seismic acquisition techniques design proposal and geological data processing is contacted directly, realized optimizing stereo observing system with this and design this acquisition technique target.
Optimizing design method of 3 D seismic observation system based on geologic geophysical model provided by the invention may further comprise the steps:
A) according to the existing exploration achievement in exploratory area, promptly complex geological structure and complicated geological volume morphing and geophysical parameters thereof are set up the geologic geophysical model of priori; B) be exploration targets with above-mentioned geologic geophysical model, design the recording geometry that two or more 3D seismic datas are gathered; C) with this model as the 3-D seismics object of observation, design desired technical parameter and construction figure in the open-air earthquake data acquisition process of laboratory simulation according to recording geometry, obtain the open-air original seismic data of simulation, the concrete practice is to adopt the simulation of the three-dimensional fluctuation of empty spectrometry equation number value just drilling with two kinds of method associatings of 3-D seismics physical model simulation to finish this data collection task; D) data of gathering are carried out 3D seismic data and handle,,, obtain earthquake and handle several imaging results the data preceding depth shift processing that change wherein based on common reflection point; E) with the geologic geophysical model of priori as estimating and contrast standard, with signal to noise ratio (S/N ratio), resolution and structural feature feature as the technical Analysis main contents, the comprehensive comparative evaluation that several imaging results are carried out image quality is handled in earthquake, through contrast, determine one of them best earthquake imaging results.
Step e wherein) also comprise:
E1) if imaging results does not also reach the target call of seismic prospecting, then revise original stereo observing system parameter, again carry out the work of 3D seismic data acquisition process subsequently, described stereo observing system parameter comprises 3-D seismics observation width degree, bin size, covering or illumination number of times, illuminance, seismic trace sum, maximum offset, wherein observes width degree N be defined as:
Wherein: L
c/ L
iThe horizontal vertical ratio of-arrangement sheet,
n
c/ n
iThe horizontal vertical ratio of-degree of covering,
C
1-with arrange the relevant coefficient of the horizontal vertical ratio of sheet,
C
2-with the relevant coefficient of the horizontal vertical ratio of degree of covering, and C
2=1-C
1, generally get C
1=C
2=0.5,
θ-center, subarea bin azimuthal distribution scope (radian);
E2), then export 3 D seismic observation system optimal design process flow diagram if imaging results has reached the target call of seismic prospecting.
The present invention is described in detail below in conjunction with accompanying drawing, so that the advantage of further understanding above-mentioned purpose of the present invention, characteristics and bringing thus.
Description of drawings
Fig. 1 is the process flow diagram of the optimizing design method of 3 D seismic observation system based on geologic geophysical model of the present invention;
Fig. 2 (a) is the Qianmiqiao buried hill structural physical model according to the priori data make;
Fig. 2 (b) is the two dimension 337 survey line log sheets that North and South direction is passed the top, buried hill;
Fig. 3 adopts the three-dimensional Wave equation forward modeling of empty spectrometry respectively and adopts the resulting seismologic record of three-dimensional seismic physical simulation on 337 surveys line;
Fig. 4 is the horizontal superposition imaging results of geological data (337 line) that three kinds of recording geometrys are gathered, and wherein Fig. 4 (a) is vertical 1 recording geometry of 24 horizontal strokes of 8 lines, 2 big guns; Fig. 4 (b) is vertical 2 recording geometrys of 24 horizontal strokes of 8 lines, 4 big guns; Fig. 4 (c) is vertical 2 recording geometrys of 24 horizontal strokes of 8 lines, 8 big guns;
Fig. 5 is the earthquake data offset imaging results (the dropping cut slice figure of cima top imaging) of different width degree three-dimensional acquisitions, and wherein, Fig. 5 (a) is the earthquake data offset of narrow three-dimensional acquisition, and Fig. 5 (b) is the earthquake data offset of wide three-dimensional acquisition.
Embodiment
It is target since 1996 at exploratory area concrete complex geological structure that physical prospecting key lab of China National Petroleum Corporation (CNPC) explored a kind of, based on the CRP seismic imaging, new three-dimensional acquisition technical design method, this method workflow as shown in Figure 1.
At first set up the geologic geophysical model of priori, set up the seismic velocity model of underground medium usually according to the existing exploration achievement in exploratory area (complex geological structure and complicated geological volume morphing and geophysical parameters thereof); Second step was an exploration targets with above-mentioned geologic geophysical model, designed the recording geometry that two or more 3D seismic datas are gathered; The 3rd goes on foot with this model as the 3-D seismics object of observation, design desired technical parameter and construction figure in the open-air earthquake data acquisition process of laboratory simulation according to recording geometry, this method characteristics are to adopt mathematical computations and two kinds of methods of Physical Experiment to simulate open-air earthquake construction, obtaining the open-air original seismic data (referring to Fig. 3) of simulation, specifically is to adopt the simulation of the three-dimensional fluctuation of empty spectrometry equation number value just drilling with two kinds of method associatings of 3-D seismics physical model simulation to finish this data collection task; The 4th step was the data of gathering to be carried out 3D seismic data handle, and handling principle is based on common reflection point (CRP), and core technology is repeatedly preceding depth shift processing, i.e. a kind of seismic imaging technology of current popular, thus obtaining earthquake handles several imaging results; The 5th step with the geologic geophysical model of priori as estimating and contrast standard, with signal to noise ratio (S/N ratio), resolution and structural feature feature as the technical Analysis main contents, the comprehensive comparative evaluation that several imaging results are carried out image quality is handled in earthquake, through contrast, determine one of them best earthquake imaging results; If this imaging results does not also reach the target call of seismic prospecting, then to revise original stereo observing system parameter (design proposal), carry out the work of 3D seismic data acquisition process subsequently again; If this imaging results has reached the target call of seismic prospecting, then its pairing stereo observing system is exactly the end product of this area's 3-d seismic exploration optimal design, exports 3 D seismic observation system optimal design workflow diagram subsequently.At this moment, whole optimal design work also declares to be finished simultaneously.
The 3 D seismic observation system main design contents is that three-dimensional observation template (several lines, a few process of preparing Chinese medicine) is selected and the stereo observing system parameter designing.The stereo observing system parameter mainly contains 3-D seismics observation width degree, bin size, covering or illumination number of times, illuminance, seismic trace sum, maximum offset etc.When definite 3-D seismics is observed width degree N, the following formula that utilizes the applicant to propose:
Wherein: L
c/ L
IThe horizontal vertical ratio of-arrangement sheet,
n
c/ n
IThe horizontal vertical ratio of-degree of covering,
C
1-with arrange the relevant coefficient of the horizontal vertical ratio of sheet,
C
2-with the relevant coefficient of the horizontal vertical ratio of degree of covering, and C
2=1-C
1, generally get C
1=C
2=0.5,
θ-center, subarea bin azimuthal distribution scope (radian).
When: N<0.5 is narrow orientation recording geometry;
N 〉=0.5 o'clock is the wide-azimuth recording geometry;
N 〉=0.85 o'clock is the omnibearing observation system;
Oil-gas exploration task in conjunction with China the western and eastern complex area, physical prospecting key lab of China National Petroleum Corporation (CNPC) (CNPC) repeatedly is applied to this method in the actual 3-D seismics technical design, huge port Qianmiqiao buried hill model for example, the contrary zone of fracture model of covering in car row area, Xinjiang, Xinjiang Yi Qikelike fold thrust belts model, the complex fault block models of Jurong, Jiangsu and Hongze etc. have all been obtained good effect.Be example with huge port Qianmiqiao buried hill structure 3 D seismic observation system optimal design below, this method is done simple the introduction.
Fig. 2 (a) is the Qianmiqiao buried hill structural physical model according to the priori data make, and Fig. 2 (b) passes the two dimension 337 survey line log sheets at top, buried hill for North and South direction.Fig. 3 adopts the three-dimensional Wave equation forward modeling of empty spectrometry respectively and adopts the resulting seismologic record of three-dimensional seismic physical simulation on 337 surveys line, relatively these two records can think that both are similar.Adopted several different recording geometrys to carry out the 3D seismic data collecting work to this area,, obtained different seismic imaging result (seeing Fig. 4 and Fig. 5) by handling.Subsequently with geologic objective as standard, various imaging results are estimated, optimize the recording geometry design.For example the signal noise ratio (snr) of image of Fig. 4 (c) should think then that than the height of Fig. 4 (b) and Fig. 4 (a) the corresponding recording geometry of Fig. 4 (c) is wherein best; For example the top imaging of Fig. 5 (b) cima is relatively good again, and focusing quality is than Fig. 5 (b) height, and explanation should select the wide-azimuth stereo observing system to carry out the earthquake construction.Some above-mentioned practical applications show: this method scope of application is extensive, and effect is good, has realistic meaning to instructing petroleum gas seismic prospecting actual production, raising oil-gas exploration success ratio.
Claims (1)
1, a kind of optimizing design method of 3 D seismic observation system based on geologic geophysical model may further comprise the steps:
A) according to the existing exploration achievement in exploratory area, i.e. complex geological structure and complicated geological volume morphing and geophysical parameters thereof are set up the geologic geophysical model of priori;
B) be exploration targets with above-mentioned geologic geophysical model, design the recording geometry that two or more 3D seismic datas are gathered;
C) with this model as the 3-D seismics object of observation, design desired technical parameter and construction figure in the open-air earthquake data acquisition process of laboratory simulation according to recording geometry, obtain the open-air original seismic data of simulation, the concrete practice is to adopt the simulation of the three-dimensional fluctuation of empty spectrometry equation number value just drilling with two kinds of method associatings of 3-D seismics physical model simulation to finish this data collection task;
D) data of above-mentioned collection are carried out 3D seismic data and handle,,, obtain earthquake and handle several imaging results the data of the gathering preceding depth shift processing that change wherein based on common reflection point;
E) with the geologic geophysical model of priori as estimating and contrast standard, with signal to noise ratio (S/N ratio), resolution and structural feature feature as the technical Analysis main contents, the comprehensive comparative evaluation that several imaging results are carried out image quality is handled in described earthquake, through contrast, determine one of them best earthquake imaging results; Step e wherein) also comprise:
E1) if imaging results does not also reach the target call of seismic prospecting, then revise original stereo observing system parameter, again carry out the work of 3D seismic data acquisition process subsequently, described stereo observing system parameter comprises 3-D seismics observation width degree, bin size, covering or illumination number of times, illuminance, seismic trace sum, maximum offset, wherein observes width degree N be defined as:
Wherein: Lc/Li-arranges the horizontal vertical ratio of sheet,
The horizontal vertical ratio of nc/ni-degree of covering,
C1-and the relevant coefficient of the arrangement horizontal vertical ratio of sheet,
C2-and the relevant coefficient of the horizontal vertical ratio of degree of covering, and C2=1-C1 generally get C1=C2=0.5,
θ-center, subarea bin azimuthal distribution scope (radian);
E2), then export 3 D seismic observation system optimal design process flow diagram if imaging results has reached the target call of seismic prospecting.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB011204451A CN1158539C (en) | 2001-07-13 | 2001-07-13 | Optimizing design method of 3D seismic observation system based on geologic geophysical model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB011204451A CN1158539C (en) | 2001-07-13 | 2001-07-13 | Optimizing design method of 3D seismic observation system based on geologic geophysical model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1337583A CN1337583A (en) | 2002-02-27 |
CN1158539C true CN1158539C (en) | 2004-07-21 |
Family
ID=4664139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB011204451A Expired - Fee Related CN1158539C (en) | 2001-07-13 | 2001-07-13 | Optimizing design method of 3D seismic observation system based on geologic geophysical model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1158539C (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102147480A (en) * | 2011-01-12 | 2011-08-10 | 中国石油集团川庆钻探工程有限公司 | Test method of three-dimensional observation system |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7286690B2 (en) * | 2003-10-23 | 2007-10-23 | Pgs Americas, Inc. | Method for imaging of pre-stack seismic data |
CN1299127C (en) * | 2004-03-26 | 2007-02-07 | 中国石油天然气集团公司 | Laminated dielectric double-focusing method for earthquake observation system optimized design |
CN100429527C (en) * | 2005-03-25 | 2008-10-29 | 中国石油天然气集团公司 | Three-dimensional seismic data processing quality monitoring technology |
CN101625417B (en) * | 2008-07-08 | 2011-09-07 | 中国石油集团东方地球物理勘探有限责任公司 | Method for optimizing design of vertical seismic profile observation system |
CN102023310B (en) * | 2009-09-19 | 2013-05-15 | 中国石油化工集团公司 | Uniformity analysis method of three-dimensional observation system |
CN102053271B (en) * | 2010-11-17 | 2012-09-26 | 中国石油天然气股份有限公司 | Earthquake collection quality analysis method and device |
CN102169187B (en) * | 2010-12-01 | 2013-01-23 | 中国科学院地质与地球物理研究所 | Quantitative analysis method for overall balance of three dimensional seismic observation system |
CN102495943B (en) * | 2011-10-31 | 2015-04-22 | 国家卫星海洋应用中心 | Modeling method for geophysical model |
CN103116183B (en) * | 2011-11-16 | 2015-05-27 | 中国石油天然气集团公司 | Method of oil earthquake collection surface element covering degree property body slicing mapping |
CN103777229A (en) * | 2012-10-23 | 2014-05-07 | 中国石油化工股份有限公司 | VSP observation system design method facing objective layer |
CN103149590B (en) * | 2013-02-26 | 2016-01-27 | 佟小龙 | geophysical imaging method and device |
CN104142518B (en) * | 2013-10-29 | 2017-05-03 | 中国石油化工股份有限公司 | Method for analyzing pre-stack time migration response of seismological observation system |
CN104635259B (en) * | 2013-11-06 | 2017-09-01 | 中国石油天然气集团公司 | It is a kind of to determine the method and device that observation system arranges piece scope |
CN103901465A (en) * | 2014-03-21 | 2014-07-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Design method of holographic three-dimensional seismic prospecting and observing system |
CN104050359B (en) * | 2014-05-30 | 2016-10-05 | 中国石油大学(华东) | A kind of the Forward Modeling based on stereo observing system arrangement sheet data segmentation |
CN104090297A (en) * | 2014-06-14 | 2014-10-08 | 吉林大学 | Reverse illumination method for optimizing earthquake collection observing system |
CN104536041B (en) * | 2014-12-17 | 2017-02-22 | 中国石油天然气集团公司 | Optimization method of seismological observation system parameters |
EP3413092B1 (en) * | 2017-06-08 | 2022-06-01 | Total Se | Method for evaluating a geophysical survey acquisition geometry over a region of interest, related process, system and computer program product |
CN111723517B (en) * | 2020-05-26 | 2022-09-02 | 河海大学 | Geophysical model optimization method integrating machine learning principle and gradual deformation method |
CN114722590B (en) * | 2022-03-23 | 2024-03-22 | 中石化石油工程技术服务有限公司 | Design optimization method of random acquisition observation system based on geophysical model |
-
2001
- 2001-07-13 CN CNB011204451A patent/CN1158539C/en not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102147480A (en) * | 2011-01-12 | 2011-08-10 | 中国石油集团川庆钻探工程有限公司 | Test method of three-dimensional observation system |
Also Published As
Publication number | Publication date |
---|---|
CN1337583A (en) | 2002-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1158539C (en) | Optimizing design method of 3D seismic observation system based on geologic geophysical model | |
CN100557464C (en) | Seismic prospecting horizon calibration method based on the prestack wave-field simulation | |
EP2250520B1 (en) | Method for indexing a subsurface volume for the purpose of inferring geologic information | |
CN1170170C (en) | Spectral decomposition for seismic interpretation | |
US8694261B1 (en) | 3D-well log invention | |
CN102937720B (en) | Well control improves the method for seismic data resolution | |
CN1189218A (en) | Method and apparatus for seismic signal processing and exploration | |
CN104678434A (en) | Method for predicting storage layer crack development parameters | |
CN106556869A (en) | Method that tomography walk sliding displacement is portrayed under a kind of extention environment quantitatively | |
CN108445538A (en) | The method and system of Depth Domain layer Q model is established based on reflected seismic information | |
CN104142518A (en) | Method for analyzing pre-stack time migration response of seismological observation system | |
CN111722284B (en) | Method for establishing speed depth model based on gather data | |
CN106777707B (en) | Method for quantitatively identifying well logging lithology by using improved spider web diagram | |
CN101359056A (en) | Method for generating longitudinal wave time domain high precision converted wave sectional plane | |
US20220221604A1 (en) | Method of analysing seismic data | |
Martuganova et al. | 3D deep geothermal reservoir imaging with wireline distributed acoustic sensing in two boreholes | |
CN107942379A (en) | A kind of method for improving complex fault block rate pattern precision | |
CN104459798A (en) | Velocity modeling method based on RTM imaging | |
Lafond et al. | Migration of wide‐aperture onshore‐offshore seismic data, central California: Seismic images of late stage subduction | |
CN113031070B (en) | Method for making depth domain synthetic seismic record | |
CN114609675A (en) | Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle | |
CN104237939B (en) | A kind of 3 D seismic observation system multi-frequency migration noise acquisition methods and device | |
AU2014394076B2 (en) | Methods and systems for identifying and plugging subterranean conduits | |
CN110954956B (en) | Method for evaluating acquisition trace of observation system and computer-readable storage medium | |
CN110531421B (en) | Sedimentary microfacies identification method |
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 | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |