CN104330825A - Fault interpretation processing method by slicing along structural trend surface - Google Patents
Fault interpretation processing method by slicing along structural trend surface Download PDFInfo
- Publication number
- CN104330825A CN104330825A CN201410375253.7A CN201410375253A CN104330825A CN 104330825 A CN104330825 A CN 104330825A CN 201410375253 A CN201410375253 A CN 201410375253A CN 104330825 A CN104330825 A CN 104330825A
- Authority
- CN
- China
- Prior art keywords
- fault
- interpretation
- trend
- along
- seismic data
- 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.)
- Pending
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a fault interpretation processing method by slicing along structural trend surface and belongs to the technical field of seismic data processing. The method is characterized by, according to development characteristics of stratums in different areas, carrying out slice research based on stratum development rules; through carrying out horizon grid interpretation on a target stratum in the seismic data, determining time range of the target stratum to obtain structural trend of the stratum; and carrying out slicing along the structural trend surface of the stratum in the three-dimensional seismic data, and the obtained slice can truly reflect the fault formation shape and trend in the period, so that the purpose of fault interpretation is achieved. The method overcomes the defects that conventional slices cannot truly reflect the fault system distribution along the structure trend, so that the influence of man-made factors on horizon slice result is prevented; minor fault is effectively identified, especially the interlayer minor fault; seismic interpretation is allowed to be more reasonable; reliable basis is provided for fault system combination; minor fault block entrapment is realized precisely; and the method provides basis for the oil field exploration and development and well deployment.
Description
Technical field
The present invention relates to a kind of along structural trend surface section fault interpretation disposal route, belong to Seismic Data Processing Technique field.
Background technology
Problem maximum in seismic interpretation is exactly fault interpretation, the degree of reliability of fault interpretation and rationality directly affect the fulfilling degree of structure, because fault pattern is numerous, scale very different, Forming Mechanism and structural setting different, therefore, in seismic interpretation, the various technological means of integrated application identifies tomography.
At present, three dimensional seismic data is utilized to explain in tomography process, comprehensively fault interpretation is completed mainly through " point ", " line ", " face ", " body ", utilize individual well point drilling well layering to complete the explanation of single-point in conjunction with earthquake, utilize wall scroll seismic section to complete on line and explain, in conjunction with geological knowledge by multi-thread explanation, complete the explanation of fault surface, on this basis, make full use of 3-D seismics advantage, seismic data volume is carried out tomography integrated interpretation.At present, in the technical method of fault interpretation, mainly through utilizing leveled time (being concerned with) section or horizon slice to instruct the explanation of tomography, leveled time (being concerned with) microtomy is applicable for the area of stratum relative level, but just there is certain problem in this technology in the area for stratum relative tilt, because stratum is not level, when production technique section, inevitable effect when having occurred wearing, well can not reflect trend and the occurrence change of tomography; And horizon slice is excessive for explanation Personnel Dependence, and detailed level position interpretation work can horizon slice be made after must be completed, and make section with explanation personnel for a kind of subjective consciousness of fault interpretation, there is certain human factor, there is certain randomness.
Summary of the invention
The object of this invention is to provide a kind of along structural trend surface section fault interpretation disposal route, adopt leveled time (being concerned with) section or horizon slice to cause well reacting fault strike and occurrence change to solve existing fault interpretation.
The present invention provides a kind of along structural trend surface section fault interpretation disposal route for solving the problems of the technologies described above, and this explanation and processing method comprises the following steps:
1) gather three dimensional seismic data, and poststack pre-service is carried out to it;
2) carry out horizontal time slice to pretreated three dimensional seismic data to browse, determine spread scope and the trend of tomography;
3) carry out layer position grid to zone of interest to explain, determine objective interval approximate time scope, obtained the structure trend of zone of interest by 3-D seismics process;
4) structural trend surface along stratum in 3-d seismic data set is cut into slices, and with the section of the true spread of reaction tomography obtained, reaches the object to fault tectonic interpretation process.
Described step 1) in poststack pre-service comprise denoising and open up frequently process.
By by step 4) section that obtains shows in three dimensions, to be concerned with result along layer really to obtain being close to.
The invention has the beneficial effects as follows: the present invention is according to stratum, different regions development characteristics, the rule of growing along stratum carries out section research, explain by carrying out layer position grid to zone of interest in seismic data, determine the time range of zone of interest, obtain the structure trend of zone of interest; Structural trend surface along stratum is cut into slices to tomography, obtains fault parameter and trend, reaches the object explained fault tectonic.The method overcome tradition section not enough to minor fault recognition capability, fracture system is portrayed unclear, can not along the defect of structure trend true reflection fracture system spread, avoid artificial explanation factor to the impact of horizon slice result, effectively identify minor fault, particularly interlayer minor fault, make seismic interpretation more reasonable, for fracture system combination provides reliable basis, meticulously implement little Fault-block trap, for In Oil Field Exploration And Development well site deployment provides foundation.
Accompanying drawing explanation
Fig. 1 is that in the embodiment of the present invention, ancient city is regional along structural trend surface (T54 seismic reflection layer) time slice schematic diagram;
Fig. 2 is area, ancient city 880ms time slice in the embodiment of the present invention (T54 fundamental purpose layer corresponding time) schematic diagram;
Fig. 3 is T54 seismic reflection layer structural map schematic diagram in area, ancient city in the embodiment of the present invention;
Fig. 4 is area, ancient city T54 seismic reflection layer the old and new fracture system contrast schematic diagram in the embodiment of the present invention;
Fig. 5 is that in the embodiment of the present invention, king collects area along structural trend surface (T53 seismic reflection layer) time slice schematic diagram;
Fig. 6 is that in the embodiment of the present invention, king collects regional 852ms time slice (T53 fundamental purpose layer corresponding time) schematic diagram;
Fig. 7 is that in the embodiment of the present invention, king collects regional T53 seismic reflection layer organigram.
Embodiment
Below in conjunction with accompanying drawing, the specific embodiment of the present invention is further described.
The present invention is directed to Problems existing in background technology, utilize prior art, by three-dimensional visualization, according to stratum, different regions development characteristics, the rule of growing along stratum sees exhibition section research, the occurrence of meticulous depiction tomography and trend, make tomography cutting relation in order, find minor fault, put fracture system spread in order, reach the object of fine structural interpretation, its concrete implementation process comprises following step.
1. pair initial three-dimensional seismic data carries out denoising and opens up process frequently, to widen seismic data effective band, improves the signal to noise ratio (S/N ratio) of seismic data.
2., on the seismic data basis handled by step 1, determined roughly spread scope and the trend of tomography by browsing time section.
3. pair fundamental purpose layer carries out the grid explanation of layer position, determines the time range of this region zone of interest, is obtained the structure trend of zone of interest by 3-D seismics.
4. the achievement of step 3 is carried out 3-D display, the section along this interpretation horizon is made in 3-d seismic data set, this section does not rely on interpretation horizon, section is made according to the structural trend surface in step 3, the section obtained like this is close to along layer, the true Distribution Characteristics being equivalent to tomography in this obtained, the spread of tomography is consistent with the tomography of final structural map, and such section better can portray the faults distribution of study area in period.
5. on the basis of step 4, show coherent data volumes in three dimensions, to obtain being close to result relevant along layer really, and without the artificial subjective consciousness of explanation personnel, more effectively can instruct interpretation of small faults and fracture system combination.
Fig. 1-Fig. 4 is the applicable cases of the present invention in area, ancient city, oil field, ancient city, Fig. 1 is that area, ancient city is along structural trend surface (T54 seismic reflection layer) time slice, area, comparison diagram 2 ancient city 880ms time slice (T54 fundamental purpose layer corresponding time), by spread and the syntagmatic of the northwestward and southeast tomography (as position of drawing a circle in figure) can well be told shown in Fig. 1, and identify that tomography is comparatively difficult by Fig. 2 (circled positions), comparison diagram 3, the section (see Fig. 1) utilizing the present invention to obtain and the fracture system degree of agreement of structural map very good, and comparison diagram 2, the a lot of tomography in position of drawing a circle in figure does not all react, mainly because Fig. 2 is based on waiting T0 interface, the form of real construction can not be reflected, and pass through the present invention, the tomography f23 explained in figure 3, the clear enough that f24 and f25 tri-tomographies show on Fig. 1, Fig. 2 does not see the vestige of tomography, absolutely prove the section utilizing the method along structural trend surface to make, reflect the spread of fracture, give full play to the advantage of three dimensional seismic data, better reflect rupture this period spread and structure looks simultaneously.
Utilize the present invention to go to launch fine structural interpretation to this, as seen from Figure 4, before not utilizing the present invention, explain this district fracture totally 54, explain after application the present invention and implement tomography 121, illustrate and utilize the present invention can effective execution tomography, put fracture spread in order.
Fig. 5-Fig. 7 is the applicable cases of the present invention at Wangji oilfield king 9 block, and specific implementation process is carried out according to the method for description.
Fig. 5 is that king collects area along structural trend surface (T53 seismic reflection layer) time slice, Fig. 6 collects regional 852ms time slice (T53 fundamental purpose layer corresponding time) for king, contrast two figure, there is certain difference in the tomography identified, pass through tectonic interpretation, obtain Fig. 7, left figure is the structural map before tackling key problem, right figure is the structural map after application the present invention, comparative analysis can be found out, after application the present invention, fault interpretation is more meticulous, with major fault, explain a collection of minor fault, these occur changes mainly because utilize the present invention can construct looks by actual response, true reflection fracture spread.By the present invention, can the more approximate fracture system portraying a certain zone of interest, better carry out structure elucidation.
Claims (3)
1., along a structural trend surface section fault interpretation disposal route, it is characterized in that, this explanation and processing method comprises the following steps:
1) gather three dimensional seismic data, and poststack pre-service is carried out to it;
2) carry out horizontal time slice to pretreated three dimensional seismic data to browse, determine spread scope and the trend of tomography;
3) carry out layer position grid to zone of interest to explain, determine objective interval approximate time scope, obtained the structure trend of zone of interest by 3-D seismics process;
4) structural trend surface along stratum in 3-d seismic data set is cut into slices, and with the section of the true spread of reaction tomography obtained, reaches the object to fault tectonic interpretation process.
2. according to claim 1 along structural trend surface section fault interpretation disposal route, to it is characterized in that, described step 1) in poststack pre-service comprise denoising and open up and frequently process.
3. according to claim 1ly it is characterized in that along structural trend surface section fault interpretation disposal route, by by step 4) section that obtains shows in three dimensions, to be concerned with result along layer really to obtain being close to.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410375253.7A CN104330825A (en) | 2014-07-31 | 2014-07-31 | Fault interpretation processing method by slicing along structural trend surface |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410375253.7A CN104330825A (en) | 2014-07-31 | 2014-07-31 | Fault interpretation processing method by slicing along structural trend surface |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104330825A true CN104330825A (en) | 2015-02-04 |
Family
ID=52405584
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410375253.7A Pending CN104330825A (en) | 2014-07-31 | 2014-07-31 | Fault interpretation processing method by slicing along structural trend surface |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104330825A (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105717540A (en) * | 2016-03-14 | 2016-06-29 | 中国海洋石油总公司 | Precise prediction method for micro-amplitude structure |
| CN106338761A (en) * | 2015-07-07 | 2017-01-18 | 中国石油化工股份有限公司 | Three-dimensional embedded entity modeling method |
| CN106556869A (en) * | 2017-01-19 | 2017-04-05 | 中国石油大学(华东) | Method that tomography walk sliding displacement is portrayed under a kind of extention environment quantitatively |
| CN106814392A (en) * | 2015-11-27 | 2017-06-09 | 中国石油化工股份有限公司 | The velocity modeling method in three-dimensional secondary closely plast fusion face |
| CN107121700A (en) * | 2017-02-28 | 2017-09-01 | 中国石油化工股份有限公司 | The preparation method of structural map under a kind of thrust-nappe belt in front of the mountains |
| CN107870364A (en) * | 2016-09-23 | 2018-04-03 | 中国石油天然气股份有限公司 | Correction method and device for seismic data interpretation data |
| CN109425898A (en) * | 2017-08-21 | 2019-03-05 | 中国石油化工股份有限公司 | A method of utilizing slice, the identification tomography that vertically and horizontally develops of section |
| CN110633557A (en) * | 2019-10-30 | 2019-12-31 | 太原理工大学 | A method for identifying favorable areas of coalbed methane structure |
| CN111287726A (en) * | 2018-12-10 | 2020-06-16 | 中国石油天然气股份有限公司 | Well position deployment area determination method and device and storage medium |
| CN112394407A (en) * | 2020-10-28 | 2021-02-23 | 中国石油天然气集团有限公司 | Method and device for predicting steep slope fan body |
| CN113589367A (en) * | 2020-04-30 | 2021-11-02 | 中国石油化工股份有限公司 | Method for correcting structure trend near fracture based on area conservation |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0956533A1 (en) * | 1997-01-31 | 1999-11-17 | Phillips Petroleum Company | Method for determining distribution of reservoir permeability, porosity and pseudo relative permeability |
| WO2003052671A1 (en) * | 2001-12-13 | 2003-06-26 | Exxonmobil Upstream Research Company | Method for locally controlling spatial continuity in geologic models |
| CN102200588A (en) * | 2011-03-22 | 2011-09-28 | 成都理工大学 | Method for analyzing waveform similarity body curvature of seismic data |
| CN103901478A (en) * | 2012-12-28 | 2014-07-02 | 中国石油天然气集团公司 | Method for determining deposition characteristics and distribution of reservoirs by combining logging and seismic information |
-
2014
- 2014-07-31 CN CN201410375253.7A patent/CN104330825A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0956533A1 (en) * | 1997-01-31 | 1999-11-17 | Phillips Petroleum Company | Method for determining distribution of reservoir permeability, porosity and pseudo relative permeability |
| WO2003052671A1 (en) * | 2001-12-13 | 2003-06-26 | Exxonmobil Upstream Research Company | Method for locally controlling spatial continuity in geologic models |
| CN102200588A (en) * | 2011-03-22 | 2011-09-28 | 成都理工大学 | Method for analyzing waveform similarity body curvature of seismic data |
| CN103901478A (en) * | 2012-12-28 | 2014-07-02 | 中国石油天然气集团公司 | Method for determining deposition characteristics and distribution of reservoirs by combining logging and seismic information |
Non-Patent Citations (3)
| Title |
|---|
| 刘金平: "松辽盆地徐家围子断陷火山岩储层综合预测方法研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
| 杨立英等: "松辽盆地南部深层火山岩、火山机构和火山岩相地质-地震综合识别", 《吉林大学学报（地球科学版）》 * |
| 郎晓玲等: "XB地区火山岩岩相划分及储层精细刻画", 《石油地球物理勘探》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106338761A (en) * | 2015-07-07 | 2017-01-18 | 中国石油化工股份有限公司 | Three-dimensional embedded entity modeling method |
| CN106338761B (en) * | 2015-07-07 | 2018-03-02 | 中国石油化工股份有限公司 | Three-dimensional inserted body modeling method |
| CN106814392B (en) * | 2015-11-27 | 2018-10-30 | 中国石油化工股份有限公司 | The velocity modeling method in three-dimensional secondary closely plast fusion face |
| CN106814392A (en) * | 2015-11-27 | 2017-06-09 | 中国石油化工股份有限公司 | The velocity modeling method in three-dimensional secondary closely plast fusion face |
| CN105717540A (en) * | 2016-03-14 | 2016-06-29 | 中国海洋石油总公司 | Precise prediction method for micro-amplitude structure |
| CN105717540B (en) * | 2016-03-14 | 2017-12-15 | 中国海洋石油总公司 | A kind of microstructures detailed predicting method |
| CN107870364A (en) * | 2016-09-23 | 2018-04-03 | 中国石油天然气股份有限公司 | Correction method and device for seismic data interpretation data |
| CN106556869A (en) * | 2017-01-19 | 2017-04-05 | 中国石油大学(华东) | Method that tomography walk sliding displacement is portrayed under a kind of extention environment quantitatively |
| CN107121700A (en) * | 2017-02-28 | 2017-09-01 | 中国石油化工股份有限公司 | The preparation method of structural map under a kind of thrust-nappe belt in front of the mountains |
| CN107121700B (en) * | 2017-02-28 | 2019-04-26 | 中国石油化工股份有限公司 | The production method of structural map under a kind of thrust-nappe belt in front of the mountains |
| CN109425898A (en) * | 2017-08-21 | 2019-03-05 | 中国石油化工股份有限公司 | A method of utilizing slice, the identification tomography that vertically and horizontally develops of section |
| CN111287726A (en) * | 2018-12-10 | 2020-06-16 | 中国石油天然气股份有限公司 | Well position deployment area determination method and device and storage medium |
| CN110633557A (en) * | 2019-10-30 | 2019-12-31 | 太原理工大学 | A method for identifying favorable areas of coalbed methane structure |
| CN110633557B (en) * | 2019-10-30 | 2023-04-14 | 太原理工大学 | A Method for Identifying Favorable Areas of Coalbed Methane Structure |
| CN113589367A (en) * | 2020-04-30 | 2021-11-02 | 中国石油化工股份有限公司 | Method for correcting structure trend near fracture based on area conservation |
| CN112394407A (en) * | 2020-10-28 | 2021-02-23 | 中国石油天然气集团有限公司 | Method and device for predicting steep slope fan body |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104330825A (en) | Fault interpretation processing method by slicing along structural trend surface | |
| CN105717540B (en) | A kind of microstructures detailed predicting method | |
| Porter et al. | Pervasive lower-crustal seismic anisotropy in Southern California: Evidence for underplated schists and active tectonics | |
| CN104345345B (en) | A kind of shale reservoir total organic carbon TOC content prediction methods | |
| CN103336297B (en) | Micro rupture vector scan method | |
| WO2017020267A1 (en) | Interpretation method of hydraulic fracturing micro-seismic event | |
| CN105334534A (en) | Low order fault interpretation method based on construction mode guidance | |
| CN103454686A (en) | Method and system for reservoir prediction based on small-scale sedimentary facies of stratigraphic slice | |
| NO318302B1 (en) | Method of identifying large-scale / small-scale structures and migration pathways in hydrocarbon-containing areas / reservoirs | |
| CN105093297B (en) | A kind of method of quality control of microseism positioning precision | |
| CN104007468B (en) | The method of volcanics space exhibition cloth is portrayed based on amplitude-differential body seismic properties | |
| CN103941286A (en) | Rapid three-dimensional horizon interpretation method based on horizontal navigation | |
| CN106709987B (en) | Dynamic construction method of three-dimensional geological profile model | |
| CN106990436B (en) | The recognition methods of karst collapse col umn and device | |
| CN105866832A (en) | Method and device for identifying low-order fault of deep tight sandstone gas reservoir | |
| CN110361795A (en) | The exploitation method of the oil-gas reservoir controlled by Strike-slip faulted | |
| CN117830561A (en) | Gold ore prospecting method based on three-dimensional geological modeling | |
| Boncio et al. | Late Quaternary faulting in the southern Matese (Italy): implications for earthquake potential and slip rate variability in the southern Apennines | |
| CN103941287A (en) | Rapid three-dimensional fault interpretation method based on horizontal navigation | |
| Obermann et al. | Seismicity and 3-D body-wave velocity models across the Hengill geothermal area, SW Iceland | |
| Lavecchia et al. | Was the Mirandola thrust really involved in the Emilia 2012 seismic sequence (northern Italy)? Implications on the likelihood of triggered seismicity effects | |
| CN105572725B (en) | A kind of ground micro-seismic monitors station distribution design method | |
| Liu et al. | Study of roof water inrush forecasting based on EM-FAHP two-factor model | |
| CN105093298A (en) | Rapidly calculating method for minor tremor data four-dimensional search and reverse time lamination | |
| Yang et al. | Fractured-vuggy reservoir characterization of carbonate, Tarim Basin, Northwest 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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150204 |
|
| RJ01 | Rejection of invention patent application after publication |