CN103226210B - A kind of entrapment implementation method - Google Patents
A kind of entrapment implementation method Download PDFInfo
- Publication number
- CN103226210B CN103226210B CN201310110998.6A CN201310110998A CN103226210B CN 103226210 B CN103226210 B CN 103226210B CN 201310110998 A CN201310110998 A CN 201310110998A CN 103226210 B CN103226210 B CN 103226210B
- Authority
- CN
- China
- Prior art keywords
- data
- utilize
- depth migration
- velocity field
- stack
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention discloses a kind of entrapment implementation method, comprise step 1: utilize pre-stack depth migration velocity field V (x, y, z) that prestack depth migration data H (x, y, z) is converted to time domain data T (x, y, z) by Depth Domain; Step 2: make an explanation into figure to using in soil deeply data, obtains waiting T0 figure; Step 3: utilize areal geology to be familiar with erection rate field; Step 4: utilize revised pre-stack depth migration velocity field V
1(x, y, z) carries out variable speed mapping to T0 figure such as using in soil deeply data and obtains final structural map H
1(x, y, z).The present inventor utilizes using in soil deeply time domain data to carry out structure elucidation, again according to the non-earthquake in study area and drilling data to the control of speed, depth shift velocity field is carried out suitably adjusted variable speed mapping and obtain structural map, so just can not rely on seismic imaging speed and implement structure, improve trap research precision.
Description
Technical field
The present invention relates to technical field of petroleum exploration, particularly a kind of entrapment implementation method.
Background technology
In the existing complicated entrapment implementation work undertaken by seismic data, mainly contain two kinds of basic skills: time domain data variable speed mapping implements trap and pre-stack depth migration entrapment implementation method in front of the mountains.But above-mentioned two kinds of methods exist following not enough:
1. time domain data variable speed mapping entrapment implementation method
The method general time domain seismic data interpretation goes out three-dimensional structure model, then utilizes variable speed mapping technology to complete vertical time and depth transfer, obtains structural map.It is little at stratigraphic dip, the region entrapment implementation precision that lateral variation in velocity is less is higher, but the area comparatively large at inclination angle, velocity variations is violent, time domain seismic data can not by accurate for geologic body migration, only accurately cannot implement trap by vertical time and depth transfer, for proposing figure folk prescription method territory settling time data offset issue in this problem prior art, after time and depth transfer, recycle mathematical algorithm by geologic body migration.
2. pre-stack depth migration entrapment implementation method
The method, owing to principle existing obvious advantage, can implement trap more accurately, and in complex area in front of the mountains, prestack depth data migration is better than time migration after stack and pre-stack time migration, can implement structure more exactly.
But, but present inventor is in the process realizing invention technical scheme in the embodiment of the present application, finds that above-mentioned technology at least exists following technical matters:
1. time migration after stack and pre-stack time migration are owing to by geologic body migration, only accurately cannot cannot implement trap by vertical time and depth transfer by the restriction of self principle;
2. figure needs artificially to set side-play amount partially, and entrapment implementation precision is restricted;
3. although pre-stack depth migration migration is better than time domain data, affecting seriously by data quality and anisotropy, often there is trap in its image taking speed, causes structural feature to be forbidden, and comparatively large with real error of boring, entrapment implementation precision is on the low side;
4. the pre-stack depth migration cycle is longer, is difficult to meet more and more faster exploration rhythm;
5. pre-stack depth migration reduced price is high, if to carry out pre-stack depth migration process consumption ratio after terminating too low for every wheel well.
Summary of the invention
The object of the present invention is to provide a kind of entrapment implementation method to solve confinement problems in prior art in oil-gas exploration process, provide one more accurately, economic, entrapment implementation method rapidly.
In order to achieve the above object, the technical solution used in the present invention is:
A kind of entrapment implementation method, is characterized in that: step 1: utilize pre-stack depth migration velocity field V (x, y, z) that prestack depth migration data H (x, y, z) is converted to time domain data T (x, y, z) by Depth Domain; Step 2: make an explanation into figure to using in soil deeply data, obtains waiting T0 figure; Step 3: utilize areal geology to be familiar with erection rate field; Step 4: utilize revised pre-stack depth migration velocity field V
1(x, y, z) carries out variable speed mapping to T0 figure such as using in soil deeply data and obtains final structural map H
1(x, y, z).
Further, pre-stack depth migration velocity field V (x, y is utilized, z) prestack depth migration data H (x, y, z) is converted to time domain data T (x by Depth Domain, y, z) computing formula is: T (x, y, z)=H (x, y, z)/V (x, y, z).
Further, utilize areal geology to be familiar with erection rate field and be specially: utilize geological knowledge control rate horizontal change rule, revise pre-stack depth migration velocity field V (x, y, z) and obtain rational velocity field V
1(x, y, z).
Further, described geological knowledge control rate horizontal change rule is specially: drilling well, non-earthquake and the outfall fan regularity of distribution.
Further, revised pre-stack depth migration velocity field V is utilized
1(x, y, z) carries out variable speed mapping to T0 figure such as using in soil deeply data and obtains final structural map H
1(x, y, z): computing formula be: H
1(x, y, z)=T (x, y, z) * V
1(x, y, z).
Compared with prior art, the beneficial effect that technical solution of the present invention produces is: the present invention utilizes pre-stack depth migration imaging speed that prestack depth migration data is converted to time domain data, using in soil deeply time domain data is utilized to carry out structure elucidation, again according to the non-earthquake in study area and drilling data to the control of speed, depth shift velocity field is carried out suitably adjusted variable speed mapping and obtain structural map, so just can not rely on seismic imaging speed and implement structure, improve trap research precision.
Accompanying drawing explanation
Accompanying drawing 1 provides the process flow diagram of a kind of entrapment implementation method of the embodiment of the present invention.
Embodiment
A kind of entrapment implementation method in the embodiment of the present invention, the method utilizes pre-stack depth migration imaging speed that prestack depth migration data is converted to time domain data, using in soil deeply time domain data is utilized to carry out structure elucidation, again according to the non-earthquake in study area and drilling data to the control of speed, depth shift velocity field is carried out suitably adjusted variable speed mapping and obtain structural map, so just can not rely on seismic imaging speed and implement structure, improve trap research precision.
Embodiment one:
As shown in Figure 1, disclosed in this method embodiment, a kind of entrapment implementation method comprises:
Step 1: utilize pre-stack depth migration velocity field V (x, y, z) that prestack depth migration data H (x, y, z) is converted to time domain data T (x, y, z) by Depth Domain;
Step 2: make an explanation into figure to using in soil deeply data, obtains waiting T0 figure;
Step 3: utilize areal geology to be familiar with erection rate field;
Step 4: utilize revised pre-stack depth migration velocity field V
1(x, y, z) carries out variable speed mapping to T0 figure such as using in soil deeply data and obtains final structural map H
1(x, y, z).
Further, pre-stack depth migration velocity field V (x, y is utilized, z) prestack depth migration data H (x, y, z) is converted to time domain data T (x by Depth Domain, y, z) computing formula is: T (x, y, z)=H (x, y, z)/V (x, y, z).
Further, utilize areal geology to be familiar with erection rate field and be specially: utilize geological knowledge control rate horizontal change rule, revise pre-stack depth migration velocity field V (x, y, z) and obtain rational velocity field V
1(x, y, z).
Further, described geological knowledge control rate horizontal change rule is specially: drilling well, non-earthquake and the outfall fan regularity of distribution.
Further, revised pre-stack depth migration velocity field V is utilized
1(x, y, z) carries out variable speed mapping to T0 figure such as using in soil deeply data and obtains final structural map H
1(x, y, z): computing formula be: H
1(x, y, z)=T (x, y, z) * V
1(x, y, z).
Embodiment two:
The skew of the forward modeling result data in KS area shows, when known formation true velocity, time migration after stack is also difficult to geologic body migration, time migration after stack trap position offsets to stratum downdip direction relative to master pattern, pre-stack depth migration then can by geologic body migration, and corresponding using in soil deeply data also can substantially by geologic body migration.
Bore from real, using in soil deeply data significantly improves data migration precision.The embodiment of the present invention this area to realize effect as follows:
1, the calibration result of drilling well shows that time migration after stack can not make high steep stratum and the mild stratum of angle migration simultaneously simultaneously, if make high steep stratum migration, the stratum that then angle is less then can overmigration, if make angle layer migration smaller, angle significantly layer again can side-play amount not enough, using in soil deeply data then can make different angle stratum migration simultaneously, and error basic controlling is within tolerance interval.
2, real boring confirms that the skew of using in soil deeply zone of interest is more reasonable, and Gao Dian, position of fault are also more identical with drilling well.
3, using in soil deeply data affects less by pre-stack depth migration velocity field details, reliable and stable.
4, exploration rhythm is greatly accelerated.The cycle of a storehouse car wheel well is in front of the mountains very long, often will more than 1 year, therefore at once above must bore next round well location after every wheel well goes out result, if wait for that pre-stack depth migration processes again at least need 2 months, and using in soil deeply is implemented about trap only needed for 2 week again.
5, greatly saved cost, the reduced price of earthquake pre-stack depth migration is high in front of the mountains, and using in soil deeply to implement trap cheap.
2010-2011 utilizes this technological achievement to carry out primary study in multiple work area in Tarim Basin, implements multiple structural trapping, and upper brill many mouthfuls of well locations, shorten research cycle, have saved processing cost, advance the gas prospecting in this district.
Embodiment three:
Gram dark 1-2 3-D seismics work area can be divided into three major constituents layers from top to bottom: salt suprastructure, gypsum-salt rock synthem, pre-salt structure layer.Salt suprastructure is high steep, and exploration zone of interest is pre-salt structure layer, buries dark.Because ground floor height is steep, time migration after stack data migration is forbidden, and has caused KS1 well salt top to be designed and has bored 1000m shallow than real.Although pre-stack depth migration migration is better than time migration after stack data, because seismic data quality is poor, ground floor height is steep, anisotropy serious, causes data to differ greatly with real boring.This district's gas reservoir is by orecontrolling factor, and the real brill that bores of KS1 well takes off zone of interest 132m gas-bearing formation, non-water breakthrough, and on pre-stack depth migration structural map, Structural range only has 50m, is not obviously inconsistent with real boring.KS101 well then pre-stack depth migration has then bored 300m dark than real.Because well error difference is large, mainly concentrate on structure axis, be difficult to simulate a suitable correction error trend surface, do not correct by Error Trend face and obtain rational structural map, cause evaluating well site deployment difficulty, Gas reservoir eraluation difficulty.And if controlled the pre-stack depth migration process carrying out a new round by new well data, be then that cycle or expense all cannot satisfy the demands.Take using in soil deeply trap method for this district in this case to have re-started trap and implement.
Specifically, the embodiment of the present invention this area to realize effect as follows:
1, pre-stack depth migration velocity field is utilized to complete prestack depth migration data using in soil deeply
It is identical that using in soil deeply data and prestack depth migration data construct vertical position, by the demarcation of well, no matter be the high steep stratum of shallow-layer or zone of interest, using in soil deeply data deviation accuracy is all better than time migration after stack data, substantially can by geologic body migration.
2, the explanation of holostrome position is carried out to using in soil deeply data
Drilling data is utilized to carry out Fine calibration to earthquake, geologic horizon.Earthquake, well data, non-seismic facies combine and set up structure elucidation model, complete whole district's structure elucidation, and what obtain each layer of position waits T0 figure.
3, the T0 figure such as each layer of position are utilized to extract each layer position average velocity planimetric map.
From zone of interest average velocity and well velocity contrast, pre-stack depth migration velocity field is obviously unreasonable in western part, work area, and KSx4 well average velocity is lower than KSx1 well average velocity, and pre-stack depth migration speed trend then in contrast.
4, utilize well speed degree, gravel bed distribution understanding to carry out plane correction to each layer average velocity planimetric map, obtain final purpose layer average velocity planimetric map.
5, utilize each interval velocity planimetric map to carry out variable speed mapping and obtain final structural map.
The wells such as grams dark 203,204,205,206,207,208 are deployed according to structural map that using in soil deeply becomes, wherein grams dark 204 and grams dark 207 all bored and take off zone of interest, two mouthfuls of drilling wells all show that structural map that using in soil deeply becomes is that position of fault or structural feature substantially all match with real boring.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (2)
1. an entrapment implementation method, is characterized in that:
Step 1: utilize pre-stack depth migration velocity field V (x, y, z) that prestack depth migration data H (x, y, z) is converted to time domain data T (x, y, z) by Depth Domain;
Step 2: the explanation of holostrome position is carried out to using in soil deeply data;
Utilize drilling data to carry out Fine calibration to earthquake, geologic horizon, earthquake, well data, non-seismic facies combine and set up structure elucidation model, complete whole district's structure elucidation, and what obtain each layer of position waits T0 figure;
Step 3: utilize areal geology to be familiar with erection rate field;
Step 4: utilize revised pre-stack depth migration velocity field V
1(x, y, z) carries out variable speed mapping to T0 figure such as using in soil deeply data and obtains final structural map H
1(x, y, z);
Utilize pre-stack depth migration velocity field V (x, y, z) by prestack depth migration data H (x, y, z) computing formula being converted to time domain data T (x, y, z) by Depth Domain is: T (x, y, z)=H (x, y, z)/V (x, y, z);
Utilize areal geology to be familiar with erection rate field to be specially:
Utilize geological knowledge control rate horizontal change rule, revise pre-stack depth migration velocity field V (x, y, z) and obtain rational velocity field V
1(x, y, z);
Described geological knowledge control rate horizontal change rule is specially: drilling well, non-earthquake and the outfall fan regularity of distribution.
2. a kind of entrapment implementation method as claimed in claim 1, is characterized in that, utilizes revised pre-stack depth migration velocity field V
1(x, y, z) carries out variable speed mapping to T0 figure such as using in soil deeply data and obtains final structural map H
1the computing formula of (x, y, z) is: H
1(x, y, z)=T (x, y, z) * V
1(x, y, z).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310110998.6A CN103226210B (en) | 2013-04-01 | 2013-04-01 | A kind of entrapment implementation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310110998.6A CN103226210B (en) | 2013-04-01 | 2013-04-01 | A kind of entrapment implementation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103226210A CN103226210A (en) | 2013-07-31 |
| CN103226210B true CN103226210B (en) | 2016-04-06 |
Family
ID=48836717
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310110998.6A Active CN103226210B (en) | 2013-04-01 | 2013-04-01 | A kind of entrapment implementation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103226210B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103955007A (en) * | 2014-05-20 | 2014-07-30 | 中国石油化工股份有限公司胜利油田分公司西部新区研究院 | Comprehensive modeling method of complicated mountain front tectonic zone and established geologic structure model |
| CN104459783B (en) * | 2014-12-11 | 2017-05-10 | 中国石油天然气集团公司 | Method and device for recognizing structural trap |
| CN104678440A (en) * | 2015-02-15 | 2015-06-03 | 山东科技大学 | Well-constrained two-dimensional seismic variable velocity field nonlinear error correction method |
| CN106569260B (en) * | 2015-10-10 | 2019-01-08 | 中国石油化工股份有限公司 | For high-dip structure area two-dimension earthquake speed change sky school drawing methods |
| CN106918839B (en) * | 2015-12-25 | 2018-12-25 | 中国石油天然气股份有限公司 | The quantitative evaluation method and device of structural trap model |
| CN107229074A (en) * | 2016-03-24 | 2017-10-03 | 中国石油化工股份有限公司 | The variable speed mapping method and apparatus in low exploration activity field marine facies exploratory area |
| CN107688199B (en) * | 2017-08-28 | 2019-03-29 | 中国石油集团东方地球物理勘探有限责任公司 | A kind of stacking velocity field construction method suitable for two-dimentional low SNR data |
| CN108051860B (en) * | 2017-11-28 | 2021-08-03 | 北京金双狐油气技术有限公司 | Method for generating error correction surface by utilizing trend surface form |
| CN113933901B (en) * | 2020-06-29 | 2024-01-30 | 中国石油天然气股份有限公司 | Depth domain seismic data splicing method and device |
| CN113219533B (en) * | 2021-05-18 | 2023-10-20 | 中海石油(中国)有限公司深圳分公司 | Pre-stack time migration velocity modeling method and device, medium and electronic equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2148838C1 (en) * | 1999-06-16 | 2000-05-10 | Нижне-Волжский научно-исследовательский институт геологии и геофизики | Seismic data processing technique |
| CN102565853A (en) * | 2011-12-20 | 2012-07-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for modeling velocity model of geologic structure |
| CN102967881A (en) * | 2012-11-12 | 2013-03-13 | 中国石油天然气股份有限公司 | Geology single-layer data depth time conversion method and device for seismic data explanation |
| CN102967882A (en) * | 2012-11-16 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for building layer velocity model of stratum |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060047429A1 (en) * | 2004-08-24 | 2006-03-02 | Adams Steven L | Method of estimating geological formation depths by converting interpreted seismic horizons from the time domain to the depth domain |
-
2013
- 2013-04-01 CN CN201310110998.6A patent/CN103226210B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2148838C1 (en) * | 1999-06-16 | 2000-05-10 | Нижне-Волжский научно-исследовательский институт геологии и геофизики | Seismic data processing technique |
| CN102565853A (en) * | 2011-12-20 | 2012-07-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for modeling velocity model of geologic structure |
| CN102967881A (en) * | 2012-11-12 | 2013-03-13 | 中国石油天然气股份有限公司 | Geology single-layer data depth time conversion method and device for seismic data explanation |
| CN102967882A (en) * | 2012-11-16 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for building layer velocity model of stratum |
Non-Patent Citations (4)
| Title |
|---|
| BS盆地北部陡坡带复杂小断块地震资料解释方法与变速成图;徐春贵等;《科技咨询》;20080303(第7期);第65页第2-3栏 * |
| 变速成图方法及应用研究;王树华等;《中国海洋大学学报》;20040130;第34卷(第1期);第139-146页 * |
| 地震速度场建立与变速构造成图的一种方法;马海珍等;《石油地球物理勘探》;20020215;第37卷(第1期);第58页第1栏第13行-第41行 * |
| 山地复杂构造地震成图方法探讨;梁顺军等;《石油物探》;20031225;第42卷(第4期);第530-537页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103226210A (en) | 2013-07-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103226210B (en) | A kind of entrapment implementation method | |
| CN106597533B (en) | A kind of Depth Domain velocity modeling method for mountain front seism processing | |
| Guoqi et al. | The pre-Sinian rift in central Sichuan Basin and its control on hydrocarbon accumulation in the overlying strata | |
| Huang et al. | Genetic types and sequence stratigraphy models of Palaeogene slope break belts in Qikou Sag, Huanghua Depression, Bohai Bay Basin, eastern China | |
| Karp et al. | Structure and stratigraphy of the Lake Albert Rift, East Africa: Observations from seismic reflection and gravity data | |
| CN104459782B (en) | Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit | |
| CN104360385A (en) | Method for comprehensively establishing initial depth interval velocity model by combining seismogeology understanding | |
| CN107783187A (en) | It is a kind of that logging speed and seismic velocity are combined to the method for establishing 3D velocity field | |
| CN103643945A (en) | Technology for thin layer lithology reservoir recognition and horizontal well drilling tracking | |
| Longde et al. | Innovation and prospect of geophysical technology in the exploration of deep oil and gas | |
| CN107065009A (en) | A kind of buried hill oil gas field underground geological map preparation method based on seismic data | |
| Aadil et al. | 3D geological modeling of Punjab platform, Middle Indus Basin Pakistan through integration of Wireline logs and seismic data | |
| CN105911592B (en) | A kind of real three dimensional seismic data long wavelength's static correcting method for boring constraint | |
| CN108121009A (en) | Development late stage complex fault block velocity-variant depth mapping method | |
| CN103377196A (en) | Lithofacies paleogeography data processing method based on geographic information system | |
| CN104459798A (en) | Velocity modeling method based on RTM imaging | |
| Bishop | Saskatchewan Exploratory Progress and Problems: Oil and Gas Occurrences | |
| Tari et al. | Seismic expression of salt tectonics in the Sab’atayn Basin, onshore Yemen | |
| Ballard et al. | An integrated approach to de-risking exploration and appraisal in structurally complex fold-and-thrust belts: application to the Incahuasi field (Bolivia) | |
| Datta et al. | Capturing multiscale heterogeneity in paralic reservoir characterization: A study in Greater Burgan Field, Kuwait | |
| REN et al. | Geomodeling technology under multifactor control | |
| Jin et al. | Control of anticline crest zone on depositional system and its geological significance for petroleum in Changshaling, Yinger Sag, Eastern Jiuquan Basin | |
| Hou et al. | Tectonic-Sedimentary Interplay of a Multi-Sourced, Structurally-Confined Turbidite System in a Foreland Basin Setting: The Pennsylvanian Lower Atoka Formation, Ouachita Mountains, USA | |
| Hou et al. | Tectonic–sedimentary interplay of a confined deepwater system in a foreland basin setting: the Pennsylvanian lower Atoka Formation, Ouachita Mountains, USA | |
| Yang et al. | Distinguishing between faults and coal collapse columns based on sediment dating: a case study of the Huainan coal field, China |
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 |