CN117930359B - Method for correcting fault shadows in post-stack seismic data by establishing high-precision velocity fields - Google Patents
Method for correcting fault shadows in post-stack seismic data by establishing high-precision velocity fields Download PDFInfo
- Publication number
- CN117930359B CN117930359B CN202410329330.9A CN202410329330A CN117930359B CN 117930359 B CN117930359 B CN 117930359B CN 202410329330 A CN202410329330 A CN 202410329330A CN 117930359 B CN117930359 B CN 117930359B
- Authority
- CN
- China
- Prior art keywords
- well
- fault
- seismic
- virtual
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000004519 manufacturing process Methods 0.000 claims abstract description 10
- 238000011160 research Methods 0.000 claims abstract description 10
- 230000001502 supplementing effect Effects 0.000 claims abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims abstract description 9
- 230000002159 abnormal effect Effects 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 238000012937 correction Methods 0.000 description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 241000256602 Isoptera Species 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 239000003209 petroleum derivative Substances 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/32—Transforming one recording into another or one representation into another
- G01V1/325—Transforming one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a method for correcting fault shadows in post-stack seismic data by establishing a high-precision velocity field, which comprises the following steps: manufacturing and calibrating the synthetic seismic records of all actual wells in the research area, and performing horizon interpretation according to the calibration result; determining the range of a fault shadow zone according to the seismic identification characteristics of the fault shadow; supplementing a virtual well in an area with the two sides not exceeding 200 meters inside and outside the boundary of the fault shadow area; making and calibrating the synthetic seismic record of the virtual well by means of the adjacent actual well data; taking the time domain interpretation horizon as space constraint, and obtaining a high-precision speed field with virtual well participation by adopting a well interpolation method; and performing time-depth conversion on the time-domain seismic data by using the high-precision velocity field to obtain accurate depth-domain seismic data. The invention adopts a mode of supplementing a virtual well to obtain a high-precision velocity field, and utilizes the velocity field to perform time-depth conversion on post-stack seismic data, so that an accurate depth-domain seismic data body and interpretation results can be obtained.
Description
Technical Field
The invention relates to the field of geophysical prospecting of petroleum and natural gas, in particular to a method for correcting fault shadows in post-stack seismic data by establishing a high-precision velocity field, which is applied to guiding marine prospecting and development.
Background
The seismic data is a data body formed by the fact that seismic waves excited by the earth surface propagate downwards, are reflected to the earth surface after encountering an underground geologic body interface and are recorded by a receiver, and the measurement unit is double-pass reflection time. In general, the stratum speed is relatively uniform, the time domain of the geologic body is approximately equivalent to the structural form of the depth domain, and the interpretation result of the time domain seismic data can reflect the structural form of the depth domain. The seismic data can be used for revealing the spread information of geologic bodies such as underground stratum, faults and the like through interpretation work.
If a velocity anomaly layer (such as thick mudstone, magma rock and the like) exists in the stratum and a large break fault is generated by breaking the stratum, the thickness of the velocity anomaly layer is suddenly changed, and the average velocity of the stratum is suddenly changed. Due to the variation of the average velocity, the time for the surface-excited seismic waves to reach a certain depth of formation will also vary during seismic data acquisition. At this time, the geologic body morphology (time domain morphology) revealed by the seismic data will not be identical to the true morphology (depth domain morphology). This phenomenon is caused by the fault fracture velocity anomaly layer, and is characterized in the seismic data by the occurrence of "pull-up and pull-down" of the same phase axis of the seismic in the triangular area below the fault, and the phenomenon is called "fault shadow". Sun Weizhao et al in forward modeling, recognition and correction of tomosynthesis: the formation principle of 'fault shadow' is researched by taking Nile Termit basin as an example, and the 'fault shadow' phenomenon is pointed out to be commonly existed in various oil-gas-containing basins, and is easily interpreted as a main fault associated minor fault in the interpretation process or interpreted as a anticline structure to deploy drilling, so that the accuracy of trap evaluation and well position deployment is seriously influenced. In recent years, a large amount of oil and gas resources are found in the ocean field in China, but because of high ocean exploration and development risks, the accuracy requirement on interpretation results of seismic data is higher, and accurate identification is particularly important to correct the fault shadow phenomenon.
The former considers that prestack depth migration imaging techniques can fundamentally suppress the "tomosynthesis" phenomenon. Forward modeling, recognition and correction of fault shadows as disclosed by Sun Weizhao et al in geophysical progress 2022: taking Nile Termit basin as an example ", it is pointed out that in areas with high precision velocity modeling conditions, developing prestack depth migration is the fundamental solution to eliminate" fault shadows ". "application research of 3D velocity modeling method based on fault and horizon constraint in eliminating fault shadows" by Peng Hailong et al, published in geophysical progress 2017, points out that the velocity field precision can be improved by utilizing horizons and faults to constrain the chromatographic velocity modeling. "Daqing chlamydospore field fault shadow seismic forward modeling and correction method" of Jiang Yan et al was disclosed in 2019 of petroleum geophysical exploration, indicating that a high-precision three-dimensional air-speed field can be established to correct "fault shadows" by means of well speed correction.
Although the former performs a series of researches on the prestack depth migration imaging technology, the actual work area is constrained by the well pattern density, a high-precision speed field is difficult to obtain, and a fault shadow phenomenon still exists in post-stack seismic data. And Jiang Yan proposes that the well speed correction mode is only suitable for a dense well pattern area, the established speed field for the dense well pattern area is not accurate enough, and a fault shadow correction method with higher precision for post-stack seismic data is to be proposed.
Disclosure of Invention
The invention aims to provide a method for correcting fault shadows in post-stack seismic data by establishing a high-precision velocity field, which is used for obtaining the true structural form of a stratum under a fault.
The technical scheme adopted for solving the technical problems is as follows: the method for correcting fault shadows in post-stack seismic data by establishing a high-precision velocity field comprises the following steps:
Step one, manufacturing and calibrating synthetic seismic records of all actual wells in a research area, and performing horizon interpretation according to calibration results;
Step two, determining the range of a fault shadow zone according to the earthquake identification characteristics of the fault shadow; the seismic identification features of fault shadows include:
In the time domain seismic data, a speed abnormal layer exists on an overlying stratum and is broken by faults; a stratum earthquake event under the fault is locally pulled up and pulled down; vertical dislocation or distortion phenomenon occurs on stratum earthquake phase axis below the fault; on the plane of the coherence attribute, a false fault which is parallel to the main fault appears on the fault lower disc;
Supplementing a virtual well in an area with the two sides not exceeding 200 meters inside and outside the boundary of the fault shadow area;
step four, manufacturing and calibrating the synthetic seismic record of the virtual well by means of the adjacent actual well data;
Step five, taking the speed of each well as a benchmark, taking a time domain interpretation horizon as a space constraint, and obtaining a high-precision speed field with virtual well participation by adopting a well interpolation method;
and step six, performing time-depth conversion on the time-domain seismic data by using the high-precision velocity field to obtain accurate depth-domain seismic data.
The method for manufacturing and calibrating the synthetic seismic record in the step one of the scheme specifically comprises the following steps:
And (3) acquiring seismic data and well curve data of a research area, utilizing a sound wave time difference curve and a density curve to manufacture synthetic seismic records of all actual wells of the research area, matching and calibrating according to the similarity of the synthetic seismic records and the actual seismic records beside the well, ensuring that the well earthquake achieves the best matching effect, establishing a time-depth relation of each actual well, and determining the corresponding positions of geological layers on a time domain seismic section.
The method for supplementing the virtual well in the scheme step three comprises the following steps:
(1) Well position setting: establishing a virtual well at the positions of the inner side and the outer side of the boundary of the fault shadow zone and lacking well control;
(2) Curve data for virtual wells: and adopting the acoustic time difference curve and the density curve of the adjacent well as the curve data of the virtual well.
The scheme comprises the following steps:
When the virtual well synthetic seismic record is manufactured, the initial speed is consistent with that of the adjacent actual well, and the acoustic time difference curve and the density curve also adopt the adjacent actual well curve; during calibration, matching the synthetic seismic record of the virtual well with the actual seismic record beside the virtual well; for a virtual well in a fault shadow region, the single well speed obtained after matching is an abnormal speed in the fault shadow region; for the virtual wells outside the fault shadow area, the single well speed obtained after matching is the normal speed outside the fault shadow area; this results in a single well velocity at each virtual well that matches the actual velocity.
Advantageous effects
1. The invention adopts a mode of supplementing the virtual well to obtain a high-precision velocity field, and the velocity field is utilized to carry out time-depth conversion on the post-stack seismic data, so that an accurate depth domain seismic data body and interpretation results can be obtained.
2. The invention corrects the shadow of the post-stack seismic data interrupt layer by establishing the high-precision velocity field, improves the accuracy of seismic data interpretation, and is suitable for seismic data interpretation work in the field of marine and land petroleum and natural gas exploration.
Drawings
FIG. 1 is a well A synthetic seismic record calibration diagram;
FIG. 2 is an example study area "tomosynthesis" impact range;
FIG. 3 is a position diagram of a virtual well replenishment;
FIG. 4 is a high precision velocity field profile;
FIG. 5 is a depth domain profile after time-depth conversion using a velocity field without supplementing a dummy well;
fig. 6 is a depth domain profile after time-depth conversion using a high precision velocity field after replenishing a virtual well.
Detailed Description
The invention is further described below with reference to the accompanying drawings:
this method of correcting "fault shadows" in post-stack seismic data by creating high-precision velocity fields, the method comprising the steps of:
step one, as shown in fig. 1, making and calibrating synthetic seismic records of all wells in an example research area, wherein the making of the synthetic seismic records uses a sound wave time difference curve and a density curve, and matching operation is carried out through the synthetic seismic records and the real seismic records beside the wells to obtain the time-depth relation of each well, and the positions of well layering on a time domain seismic section are clarified, so that the paraxial horizon interpretation of seismic data is developed.
Step two, determining the range of a shadow zone according to the earthquake recognition characteristics of 'fault shadows': and extracting coherence properties from the seismic data to obtain a layer-along coherence slice. And determining the influence range of the fault shadow zone according to the position of the pseudo fault on the plane coherent slice, the 'pull-up and pull-down' range of the seismic event below the fault on the section and the position of the dislocation or distortion of the seismic event. As shown in fig. 2, the left graph is a slice along the layer, the distribution of faults is visible on the slice, the right graph is a seismic reflection section, the triangle area below the faults has the characteristic of "pull-up", and the right boundary of the triangle area has the characteristic of dislocation of the same phase axis, so that the range of "fault shadows" is identified.
And thirdly, supplementing the virtual well in an area with the two sides not exceeding 200 meters inside and outside the boundary of the fault shadow area. As shown in fig. 3, the actual well positions of the example study area are unevenly distributed, so that virtual wells need to be supplemented at positions which are on the inner side and the outer side of the boundary of the shadow area and lack of well control, and an acoustic time difference curve and a density curve of adjacent wells are adopted as curve data of the virtual wells.
Step four, manufacturing and calibrating the synthetic seismic record of the virtual well by means of the adjacent actual well data; and (3) making a synthetic seismic record of the virtual well, wherein the initial speed is consistent with that of an adjacent actual well, and an adjacent well curve is adopted by the acoustic time difference and density curve. And matching the synthetic seismic record of the virtual well with the actual seismic record beside the virtual well, wherein the speed of the virtual well is the speed of the real stratum. The adjacent well is an adjacent actual well.
And fifthly, taking the speeds of all the wells as references, taking a time domain interpretation horizon as a space constraint, adopting a well interpolation method to obtain a high-precision speed field with participation of the virtual wells, and as shown in fig. 4, uniformly changing the high-precision speed field along a same phase axis, generating speed mutation at a shadow boundary, and generating abnormal high speed at a shadow region, so as to meet the speed change rule of a fault shadow phenomenon. Each well includes an actual well and a dummy well, in which a well is a well in a hatched area and B well is a well in a hatched area.
And step six, performing time-depth conversion on the time-domain seismic data by using the high-precision velocity field to obtain accurate depth-domain seismic data. As shown, FIG. 5 is a depth domain profile after a time-depth transition using velocity fields without supplementing the pseudo well, with "fault shadowing" phenomena uncorrected, and with the seismic reflection axis indicated by the arrows still having a distortion. And FIG. 6 is a depth domain section after time-depth conversion by using a high-precision velocity field after supplementing a virtual well, the common phase axis and the horizon of the earthquake are smooth, the correction effect of the 'fault shadow' phenomenon is good, and the real geological situation is met.
Claims (2)
1. A method for correcting fault shadows in post-stack seismic data by establishing a high-precision velocity field, comprising the steps of:
Step one, manufacturing and calibrating synthetic seismic records of all actual wells in a research area, and performing horizon interpretation according to calibration results;
Step two, determining the range of a fault shadow zone according to the earthquake identification characteristics of the fault shadow; the seismic identification features of fault shadows include:
In the time domain seismic data, a speed abnormal layer exists on an overlying stratum and is broken by faults; a stratum earthquake event under the fault is locally pulled up and pulled down; vertical dislocation or distortion phenomenon occurs on stratum earthquake phase axis below the fault; on the plane of the coherence attribute, a false fault parallel to the main fault appears on the fault lower disc;
Supplementing a virtual well in an area with the two sides not exceeding 200 meters inside and outside the boundary of the fault shadow area:
(1) Well position setting: establishing a virtual well at the positions of the inner side and the outer side of the boundary of the fault shadow zone and lacking well control;
(2) Curve data for virtual wells: adopting an acoustic time difference curve and a density curve of an adjacent well as curve data of a virtual well;
step four, manufacturing and calibrating the synthetic seismic record of the virtual well by means of the adjacent actual well data:
When the virtual well synthetic seismic record is manufactured, the initial speed is consistent with that of the adjacent actual well, and the acoustic time difference curve and the density curve also adopt the adjacent actual well curve; during calibration, matching the synthetic seismic record of the virtual well with the actual seismic record beside the virtual well; for a virtual well in a fault shadow region, the single well speed obtained after matching is an abnormal speed in the fault shadow region; for the virtual wells outside the fault shadow area, the single well speed obtained after matching is the normal speed outside the fault shadow area; obtaining the single well speed at each virtual well which is consistent with the actual speed;
Step five, taking the speed of each well as a benchmark, taking a time domain interpretation horizon as a space constraint, and obtaining a high-precision speed field with virtual well participation by adopting a well interpolation method;
and step six, performing time-depth conversion on the time-domain seismic data by using the high-precision velocity field to obtain accurate depth-domain seismic data.
2. The method of correcting fault shadows in post-stack seismic data by establishing a high-precision velocity field as recited in claim 1, wherein: the method for manufacturing and calibrating the synthetic seismic record in the first step comprises the following steps:
And (3) acquiring seismic data and well curve data of a research area, utilizing a sound wave time difference curve and a density curve to manufacture synthetic seismic records of all actual wells of the research area, matching and calibrating according to the similarity of the synthetic seismic records and the actual seismic records beside the well, ensuring that the well earthquake achieves the best matching effect, establishing a time-depth relation of each actual well, and determining the corresponding positions of geological layers on a time domain seismic section.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410329330.9A CN117930359B (en) | 2024-03-21 | 2024-03-21 | Method for correcting fault shadows in post-stack seismic data by establishing high-precision velocity fields |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410329330.9A CN117930359B (en) | 2024-03-21 | 2024-03-21 | Method for correcting fault shadows in post-stack seismic data by establishing high-precision velocity fields |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117930359A CN117930359A (en) | 2024-04-26 |
CN117930359B true CN117930359B (en) | 2024-06-14 |
Family
ID=90766808
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410329330.9A Active CN117930359B (en) | 2024-03-21 | 2024-03-21 | Method for correcting fault shadows in post-stack seismic data by establishing high-precision velocity fields |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117930359B (en) |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109387867B (en) * | 2017-08-10 | 2020-05-22 | 中国石油化工股份有限公司 | Compact sandstone reservoir modeling method |
CN108828666B (en) * | 2018-08-07 | 2019-12-27 | 中国石油大学(华东) | Method for recovering buried history of stratum near normal fault |
CN109828305A (en) * | 2019-03-13 | 2019-05-31 | 西安恒泰艾普能源发展有限公司 | The prediction technique of deep reservoir in a kind of shorter situation of well logging sound wave curve |
CN110673209B (en) * | 2019-10-13 | 2021-06-04 | 东北石油大学 | Well-seismic calibration method |
CN112255686B (en) * | 2020-10-14 | 2023-08-11 | 东北石油大学 | Fault edge speed modeling method based on regression algorithm |
CN114578428A (en) * | 2020-12-01 | 2022-06-03 | 中国石油天然气集团有限公司 | Efficient construction method and device for velocity field |
CN112558180B (en) * | 2021-01-05 | 2023-09-15 | 东北石油大学 | Method for rapidly checking seismic horizon calibration accuracy by utilizing horizontal isochronal plane |
CN114966856B (en) * | 2022-08-02 | 2022-12-02 | 中国科学院地质与地球物理研究所 | Carbon sequestration site location optimization method, system and equipment based on multiband seismic data |
-
2024
- 2024-03-21 CN CN202410329330.9A patent/CN117930359B/en active Active
Non-Patent Citations (1)
Title |
---|
应用空间变速度体及虚拟井参与速度成图;张述;《内蒙古石油化工》;20171231(第6期);第54-55页 * |
Also Published As
Publication number | Publication date |
---|---|
CN117930359A (en) | 2024-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Landro et al. | The Gullfaks 4D seismic study | |
US8665667B2 (en) | Vertical seismic profiling velocity estimation method | |
US8705317B2 (en) | Method for imaging of targeted reflectors | |
US20100118654A1 (en) | Vertical seismic profiling migration method | |
CN102385066B (en) | Pre-stack earthquake quantitative imaging method | |
Ramachandran et al. | Imaging permafrost velocity structure using high resolution 3D seismic tomography | |
CN117930359B (en) | Method for correcting fault shadows in post-stack seismic data by establishing high-precision velocity fields | |
CN117930360B (en) | Method for correcting fault shadows in post-stack seismic data | |
Tilbury et al. | Pluto 4D—Australia’s first 4D over a gas field is an outstanding success | |
CN117950053B (en) | Method for identifying fault shadows in seismic data by using average speed | |
Giustiniani et al. | 3D seismic data for shallow aquifers characterisation | |
Buia et al. | Depth imaging Coil data: Multi azimuthal tomography earth model building and depth imaging the full azimuth Tulip coil project | |
CN108535779B (en) | A kind of optimization method for seismic exploration data | |
CN113589367B (en) | Method for correcting nearby structure trend based on area conservation and large fracture | |
Vigh et al. | An offshore Gulf of Mexico case study applying full-waveform inversion | |
Smit | Experiences with OBS node technology in the Greater Mars Basin | |
CN110941029B (en) | Speed modeling method related to geological capping | |
Zhang | Analysis of factors affecting residual moveout picking and solutions | |
Dalętka | Selected aspects of modern seismic imaging and near-surface velocity model building in the area of Carpathian fold and thrust belt | |
Alhajni | Comparison study between the elevation and datum statics in NC 210, western Libya | |
Roth et al. | Challenges related to oyster floatstones in Vaca Muerta development: Visualization on PSTM, PSDM and diffraction imaging | |
Ting et al. | The method and application of anisotropic velocity model building using trend velocity constraints | |
Zhao et al. | Foothills seismic imaging for deep exploration in Junggar Basin—A case study | |
Murty et al. | Delineation of the Trap and sub-trappean sediments in Kutch, Deccan syneclise and Bengal basins-An analysis | |
Gonzalez et al. | Refining 3-D velocity models for depth migration using tomography: Application to rapid permafrost variations in Alaska's North Slope transition zone province |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |