CN112711071A - Stratum inclination angle correction method and device - Google Patents

Stratum inclination angle correction method and device Download PDF

Info

Publication number
CN112711071A
CN112711071A CN201911021736.6A CN201911021736A CN112711071A CN 112711071 A CN112711071 A CN 112711071A CN 201911021736 A CN201911021736 A CN 201911021736A CN 112711071 A CN112711071 A CN 112711071A
Authority
CN
China
Prior art keywords
frequency
seismic
data
vertical
thickness
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911021736.6A
Other languages
Chinese (zh)
Other versions
CN112711071B (en
Inventor
林腾飞
王雪玲
田中元
刘杏芳
罗贝维
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201911021736.6A priority Critical patent/CN112711071B/en
Publication of CN112711071A publication Critical patent/CN112711071A/en
Application granted granted Critical
Publication of CN112711071B publication Critical patent/CN112711071B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • 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/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • 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
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/52—Move-out correction
    • 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/62—Physical property of subsurface
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

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 discloses a method and a device for correcting a stratum inclination angle, wherein the method comprises the following steps: carrying out frequency spectrum decomposition according to seismic data of a preset work area to obtain a frequency spectrum decomposition result; acquiring seismic dip angle data according to the seismic data; according to the relation among the real thickness of the stratum, the vertical thickness and the seismic dip data, the relation among the real frequency in the frequency domain, the vertical frequency and the seismic dip data is obtained by combining the seismic data; and correcting the frequency spectrum decomposition result by using the relation among the real frequency, the vertical frequency and the seismic dip angle data in the frequency domain. The invention can finish the correction of the frequency spectrum decomposition result and ensure the accuracy of the frequency spectrum decomposition result.

Description

Stratum inclination angle correction method and device
Technical Field
The invention relates to the technical field of oil and gas exploration, in particular to a method and a device for correcting a stratum dip angle.
Background
The frequency spectrum decomposition technology is an interpretation method of a frequency domain, is an important component in seismic attribute analysis, and is used for researching a thin layer structure and judging a deposition environment by using the change relation of frequency along with time.
Currently, the frequency in the spectral decomposition technique is calculated according to seismic traces, and the corresponding thickness is the vertical thickness of the formation. Under the condition that the horizontal or inclined angle of the stratum is ignored, the difference between the real thickness of the stratum and the vertical thickness is not large. However, when the formation is inclined, the true thickness of the formation has a proportional relationship with the vertical thickness and the formation dip angle. And its true thickness will be less than the vertical thickness. The greater the formation dip, the greater the difference between the true thickness of the formation and the vertical thickness, and the smaller the corresponding tuning frequency. The result of the spectral decomposition will be shifted overall towards lower frequencies due to the presence of the dip. Therefore, when the formation dip is present, it is important to correct the result of the spectral decomposition, and the prior art does not have such a correction method.
Disclosure of Invention
The embodiment of the invention provides a stratum dip angle correction method, which is used for completing the correction of a frequency spectrum decomposition result and ensuring the accuracy of the frequency spectrum decomposition result, and comprises the following steps:
carrying out frequency spectrum decomposition according to seismic data of a preset work area to obtain a frequency spectrum decomposition result;
acquiring seismic dip angle data according to the seismic data;
according to the relation among the real thickness of the stratum, the vertical thickness and the seismic dip data, combining the seismic data to obtain the relation among the real frequency, the vertical frequency and the seismic dip data in the frequency domain;
and correcting the frequency spectrum decomposition result by using the relation among the real frequency, the vertical frequency and the seismic dip angle data in the frequency domain.
Optionally, the relationship between the true thickness of the formation, the vertical thickness, and the seismic dip data is:
da＝dr/cosθ；
wherein d isaIs a vertical thickness, drIs true thickness and θ is seismic dip data.
Optionally, the relationship between the true frequency, the vertical frequency and the seismic dip data in the frequency domain is as follows:
Figure BDA0002247425690000021
where f is the true frequency, faIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the formation, λ is the true wavelength of the seismic waves, and drHas a preset proportional relation.
Optionally, the vertical frequency faThe calculation formula of (a) is as follows:
Figure BDA0002247425690000022
wherein f isaIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the stratum, lambdaaIs the vertical wavelength of the seismic wave, its and daHas a preset proportional relation.
The embodiment of the invention also provides a stratum inclination angle correction device, which is used for completing the correction of the frequency spectrum decomposition result and ensuring the accuracy of the frequency spectrum decomposition result, and the device comprises:
the frequency spectrum decomposition module is used for carrying out frequency spectrum decomposition according to the seismic data of the preset work area to obtain a frequency spectrum decomposition result;
the data acquisition module is used for acquiring seismic dip angle data according to the seismic data;
the relation derivation module is used for acquiring the relation among the real frequency in the frequency domain, the vertical frequency and the seismic dip angle data by combining the seismic data according to the relation among the real thickness of the stratum, the vertical thickness and the seismic dip angle data;
and the correction module is used for correcting the frequency spectrum decomposition result by utilizing the relation among the real frequency, the vertical frequency and the seismic dip angle data in the frequency domain.
Optionally, the relationship between the true thickness of the formation, the vertical thickness, and the seismic dip data is:
da＝dr/cosθ；
wherein d isaIs a vertical thickness, drIs true thickness and θ is seismic dip data.
Optionally, the relationship between the true frequency, the vertical frequency and the seismic dip data in the frequency domain is as follows:
Figure BDA0002247425690000023
where f is the true frequency, faIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the formation, λ is the true wavelength of the seismic waves, and drHas a preset proportional relation.
Optionally, the vertical frequency faThe calculation formula of (a) is as follows:
Figure BDA0002247425690000024
wherein f isaIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the stratum, lambdaaIs the vertical wavelength of the seismic wave, its and daHas a preset proportional relation.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the method when executing the computer program.
An embodiment of the present invention further provides a computer-readable storage medium, in which a computer program for executing the above method is stored.
In the embodiment of the invention, the relationship among the real frequency in the frequency domain, the vertical frequency and the seismic dip angle data can be obtained by obtaining the frequency spectrum decomposition result and the seismic dip angle data according to the relationship among the real thickness of the stratum, the vertical thickness and the seismic dip angle data and combining the seismic data, so that the correction of the frequency spectrum decomposition result is completed, the accuracy of the frequency spectrum decomposition result is ensured, and a foundation is laid for the follow-up study of a thin layer structure and the judgment of a deposition environment.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
FIG. 1 is a flow chart of a method of correcting formation dip in an embodiment of the invention;
FIG. 2 is a schematic diagram of a structure of a formation dip correction apparatus according to an embodiment of the present invention;
FIG. 3 is an illustration of a dip in a formation in an embodiment of the present invention;
FIG. 4 is a two-dimensional tilted stratigraphic model diagram according to an embodiment of the present invention;
FIG. 5 is a schematic time-migration diagram of a test data seismic section in an embodiment of the present invention;
FIG. 6 is a schematic depth migration diagram of a test data seismic section in an embodiment of the present invention;
FIG. 7 is a cross-sectional view of uncorrected peak frequencies for test data in a frequency domain according to an embodiment of the present invention;
FIG. 8 is a cross-sectional view of uncorrected peak frequencies for test data in the depth domain in an embodiment of the present invention;
FIG. 9 is a cross-sectional view of the peak frequency of the test data after calibration in the frequency domain according to an embodiment of the present invention;
FIG. 10 is a cross-sectional view of the peak frequency of the test data in the depth domain after being corrected according to the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
The frequency in the spectral decomposition technique is calculated according to the seismic traces, and the corresponding thickness is the vertical thickness of the stratum. Under the condition that the horizontal or inclined angle of the stratum is ignored, the difference between the real thickness of the stratum and the vertical thickness is not large. However, when the formation is inclined, the true thickness of the formation has a proportional relationship with the vertical thickness and the formation dip angle. And its true thickness will be less than the vertical thickness. The greater the formation dip, the greater the difference between the true thickness of the formation and the vertical thickness, and the smaller the corresponding tuning frequency. The result of the spectral decomposition will be shifted overall towards lower frequencies due to the presence of the dip.
Fig. 1 is a flowchart of a method for correcting a formation dip according to an embodiment of the present invention, as shown in fig. 1, the method includes:
step 101, performing spectrum decomposition according to seismic data of a preset work area to obtain a spectrum decomposition result.
In the present embodiment, the spectral decomposition result can be understood as the vertical frequency mentioned below.
And 102, acquiring seismic dip angle data according to the seismic data.
In this embodiment, in the depth domain, the inclination of the seismic inclination in the X-axis direction is represented as θxThe inclination angle of the seismic inclination angle in the X-axis direction is represented as thetay. In the time domain, the inclination angle of the earthquake inclination angle in the X-axis direction is recorded as p, the inclination angle of the earthquake inclination angle in the Y-axis direction is recorded as q, and then theta is calculatedx、θyP, q are as follows:
Figure BDA0002247425690000041
q＝2tanθy/v＝ky/ω；
Figure BDA0002247425690000042
wherein the content of the first and second substances,kxand kyThe instantaneous wave numbers of the seismic waves in the X-axis direction and the Y-axis direction are respectively represented.
And 103, acquiring the relation among the real frequency, the vertical frequency and the seismic dip angle data in the frequency domain by combining the seismic data according to the relation among the real thickness, the vertical thickness and the seismic dip angle data of the stratum.
In this embodiment, fig. 3 is an exemplary diagram of the formation inclination, and as can be seen from fig. 3, the relationship between the true thickness and the vertical thickness of the formation and the seismic dip data is:
da＝dr/cosθ；
wherein d isaIs a vertical thickness, drIs true thickness and θ is seismic dip data.
The relationship between the true frequency, vertical frequency and seismic dip data in the frequency domain is:
Figure BDA0002247425690000043
where f is the true frequency, faIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the formation, λ is the true wavelength of the seismic waves, and drHas a preset proportional relation.
The vertical frequency faThe calculation formula of (a) is as follows:
Figure BDA0002247425690000044
wherein f isaIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the stratum, lambdaaIs the vertical wavelength of the seismic wave, its and daHas a preset proportional relation.
It can be understood that the values of λ and d are due torHaving a predetermined proportional relationship, λaAnd daThe method has a preset proportional relation, so that the method can be used for measuring the relation between the real thickness and the vertical thickness of the stratum under the depth domain and the seismic dip angle dataAnd deducing the relation among the real frequency, the vertical frequency and the seismic dip angle data in the frequency domain.
FIG. 4 is a two-dimensional tilted formation model diagram, in FIG. 4, the vertical thickness of the formation is 100 feet, the tuning frequency is 25 Hz, and the corresponding velocity is 10000 feet/second. However, we can find that the dip angles of the sections 1, 2 and 3 in fig. 4 are different, the vertical thickness is the same (100 feet), and the true thickness should be the vertical thickness divided by the dip angle, and this method can also be applied in the frequency domain.
And step 104, correcting the frequency spectrum decomposition result by using the relation among the real frequency in the frequency domain, the vertical frequency and the seismic dip angle data.
According to the stratum inclination angle correction method provided by the embodiment of the invention, the frequency spectrum decomposition result and the seismic inclination angle data are obtained, and the relation among the real frequency in the frequency domain, the vertical frequency and the seismic inclination angle data can be obtained according to the relation among the real thickness of the stratum, the vertical thickness and the seismic inclination angle data by combining the seismic data, so that the correction of the frequency spectrum decomposition result is completed, the accuracy of the frequency spectrum decomposition result is ensured, and a foundation is laid for the follow-up study of a thin layer structure and the judgment of a deposition environment.
When theta is small, f is approximately equal to faAnd has little influence on the frequency decomposition result. In addition, because the river reservoir is always relatively flat (segment 1 in fig. 4) or has only a small inclination angle (segment 3 in fig. 4), we can obtain good results, and as shown in fig. 3, when the inclination angle is large enough that we cannot ignore, we must correct the video frequency obtained by the vertical thickness to be the real frequency.
The following test data are examples to illustrate some aspects of the invention:
FIG. 5 shows time migration data for a test data seismic section, while FIG. 6 shows depth migration data for a test data seismic section. The shaded coverage is the same location in both sections, and it can be seen that the stratigraphic layers in the depth offset data are much deeper than the stratigraphic layers in the time offset data. This is because the velocity increases with time (depth). The sampling increments for the time offset data and the depth offset data are 0.002 seconds and 0.01 kilometer, respectively.
By respectively calculating the peak frequency (wave number) and the dip angle of the time domain and the depth domain offset data, so as to obtain the corresponding real frequency (wavelength) through calculation, it can be obtained that the peak frequency (wave number) is almost not different before and after the dip angle correction in the region with smaller formation dip angle; however, in the dip angle increasing region (both sides of the section), the peak frequency (wave number) is obviously increased, which has great significance for predicting the thickness of the inclined thin stratum. And, the result of the inclination correction is more obvious because the inclination of the depth domain offset data is relatively larger.
FIG. 7 is a cross-sectional view of uncorrected peak frequencies of test data in the frequency domain, FIG. 8 is a cross-sectional view of uncorrected peak frequencies of test data in the depth domain, FIG. 9 is a cross-sectional view of corrected peak frequencies of test data in the frequency domain, and FIG. 10 is a cross-sectional view of corrected peak frequencies of test data in the depth domain. It is obvious that the corrected test data peak frequency is closer to the true frequency.
Based on the same inventive concept, the embodiments of the present invention also provide a formation inclination correction apparatus, as described in the following embodiments. Because the principle of the formation dip angle correction device for solving the problems is similar to that of the formation dip angle correction method, the implementation of the formation dip angle correction device can refer to the implementation of the formation dip angle correction method, and repeated parts are not described again. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.
An embodiment of the present invention provides a formation dip angle correction device, as shown in fig. 2, the device includes:
the spectrum decomposition module 201 is configured to perform spectrum decomposition according to seismic data of a preset work area to obtain a spectrum decomposition result;
a data obtaining module 202, configured to obtain seismic dip angle data according to the seismic data;
the relation derivation module 203 is configured to obtain a relation between the real frequency in the frequency domain, the vertical frequency and the seismic dip angle data according to a relation between the real thickness of the stratum, the vertical thickness and the seismic dip angle data, in combination with the seismic data;
and the correcting module 204 is configured to correct the spectrum decomposition result by using a relationship between the real frequency in the frequency domain, the vertical frequency, and the seismic dip data.
In the embodiment of the invention, the relation between the real thickness and the vertical thickness of the stratum and the seismic dip angle data is as follows:
da＝dr/cosθ；
wherein d isaIs a vertical thickness, drIs true thickness and θ is seismic dip data.
In the embodiment of the invention, the relationship among the real frequency, the vertical frequency and the seismic dip data in the frequency domain is as follows:
Figure BDA0002247425690000061
where f is the true frequency, faIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the formation, λ is the true wavelength of the seismic waves, and drHas a preset proportional relation.
In the embodiment of the invention, the vertical frequency faThe calculation formula of (a) is as follows:
Figure BDA0002247425690000062
wherein f isaIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the stratum, lambdaaIs the vertical wavelength of the seismic wave, its and daHas a preset proportional relation.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the above method when executing the computer program.
An embodiment of the present invention further provides a computer-readable storage medium, in which a computer program for executing the above method is stored.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (10)

1. A method of correcting formation dip, comprising:
carrying out frequency spectrum decomposition according to seismic data of a preset work area to obtain a frequency spectrum decomposition result;
acquiring seismic dip angle data according to the seismic data;
according to the relation among the real thickness of the stratum, the vertical thickness and the seismic dip data, combining the seismic data to obtain the relation among the real frequency, the vertical frequency and the seismic dip data in the frequency domain;
and correcting the frequency spectrum decomposition result by using the relation among the real frequency, the vertical frequency and the seismic dip angle data in the frequency domain.
2. The method of claim 1, wherein the relationship between true thickness of the earth formation, vertical thickness, and seismic dip data is:
da＝dr/cosθ；
wherein d isaIs a vertical thickness, drIs true thickness and θ is seismic dip data.
3. The method of claim 2, wherein the relationship between true frequency, vertical frequency and seismic dip data in the frequency domain is:
Figure FDA0002247425680000011
where f is the true frequency, faIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the formation, λ is the true wavelength of the seismic waves, and drHas a preset proportional relation.
4. Method according to claim 3, wherein said vertical frequency faThe calculation formula of (a) is as follows:
Figure FDA0002247425680000012
wherein f isaIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the formation, λ a is the vertical wavelength of the seismic waves, and daHas a preset proportional relation.
5. A formation dip angle correcting device, comprising:
the frequency spectrum decomposition module is used for carrying out frequency spectrum decomposition according to the seismic data of the preset work area to obtain a frequency spectrum decomposition result;
the data acquisition module is used for acquiring seismic dip angle data according to the seismic data;
the relation derivation module is used for acquiring the relation among the real frequency in the frequency domain, the vertical frequency and the seismic dip angle data by combining the seismic data according to the relation among the real thickness of the stratum, the vertical thickness and the seismic dip angle data;
and the correction module is used for correcting the frequency spectrum decomposition result by utilizing the relation among the real frequency, the vertical frequency and the seismic dip angle data in the frequency domain.
6. The apparatus of claim 5, wherein the relationship between the true thickness of the earth formation, the vertical thickness, and the seismic dip data is:
da＝dr/cosθ；
wherein d isaIs a vertical thickness, drIs true thickness and θ is seismic dip data.
7. The apparatus of claim 6, wherein the relationship between true frequency, vertical frequency and seismic dip data in the frequency domain is:
Figure FDA0002247425680000021
where f is the true frequency, faIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the formation, λ is the true wavelength of the seismic waves, and drHas a preset proportional relation.
8. The apparatus of claim 7, wherein the vertical frequency faThe calculation formula of (a) is as follows:
Figure FDA0002247425680000022
wherein f isaIs the vertical frequency, T is the period, v is the propagation velocity of seismic waves in the formation, λ a is the vertical wavelength of the seismic waves, and daHas a preset proportional relation.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 4 when executing the computer program.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 4.
CN201911021736.6A 2019-10-25 2019-10-25 Stratum inclination angle correction method and device Active CN112711071B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911021736.6A CN112711071B (en) 2019-10-25 2019-10-25 Stratum inclination angle correction method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911021736.6A CN112711071B (en) 2019-10-25 2019-10-25 Stratum inclination angle correction method and device

Publications (2)

Publication Number Publication Date
CN112711071A true CN112711071A (en) 2021-04-27
CN112711071B CN112711071B (en) 2023-04-07

Family

ID=75541433

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911021736.6A Active CN112711071B (en) 2019-10-25 2019-10-25 Stratum inclination angle correction method and device

Country Status (1)

Country Link
CN (1) CN112711071B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6735527B1 (en) * 2003-02-26 2004-05-11 Landmark Graphics Corporation 3-D prestack/poststack multiple prediction
CN102109611A (en) * 2009-12-23 2011-06-29 中国石油天然气集团公司 Fast and convenient method for predicting high-quality petroleum reservoir in virtue of seism attributes
CN102466819A (en) * 2010-11-03 2012-05-23 中国石油天然气集团公司 Spectrum analysis method of seismic signal and apparatus thereof
CN103454679A (en) * 2013-08-12 2013-12-18 中国石油天然气股份有限公司 Method and system for determining seismic reflection event isochronism
WO2015102640A1 (en) * 2014-01-03 2015-07-09 Halliburton Energy Services, Inc. Systems and methods for relative dip correction
CN104950334A (en) * 2015-06-16 2015-09-30 中国石油天然气集团公司 Method and device for predicting reservoir distribution
CN105425299A (en) * 2015-11-17 2016-03-23 中国石油天然气股份有限公司 Method and device for determining formation fracture distribution
CN106918838A (en) * 2017-03-06 2017-07-04 中国科学院地质与地球物理研究所 Gaussian beam offset imaging method and device under the conditions of relief surface

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6735527B1 (en) * 2003-02-26 2004-05-11 Landmark Graphics Corporation 3-D prestack/poststack multiple prediction
CN102109611A (en) * 2009-12-23 2011-06-29 中国石油天然气集团公司 Fast and convenient method for predicting high-quality petroleum reservoir in virtue of seism attributes
CN102466819A (en) * 2010-11-03 2012-05-23 中国石油天然气集团公司 Spectrum analysis method of seismic signal and apparatus thereof
CN103454679A (en) * 2013-08-12 2013-12-18 中国石油天然气股份有限公司 Method and system for determining seismic reflection event isochronism
WO2015102640A1 (en) * 2014-01-03 2015-07-09 Halliburton Energy Services, Inc. Systems and methods for relative dip correction
CN104950334A (en) * 2015-06-16 2015-09-30 中国石油天然气集团公司 Method and device for predicting reservoir distribution
CN105425299A (en) * 2015-11-17 2016-03-23 中国石油天然气股份有限公司 Method and device for determining formation fracture distribution
CN106918838A (en) * 2017-03-06 2017-07-04 中国科学院地质与地球物理研究所 Gaussian beam offset imaging method and device under the conditions of relief surface

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨林: "地震频谱分解技术应用中有关问题的讨论", 《石油物探》 *

Also Published As

Publication number Publication date
CN112711071B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
US8483964B2 (en) Estimating reservoir properties from 4D seismic data
CN105182408A (en) Manufacturing method and device for synthesizing earthquake record
US20120140593A1 (en) Time-lapse seismic comparisons using pre-stack imaging and complex wave field comparisons to improve accuracy and detail
US20170184748A1 (en) A method and a computing system for seismic imaging a geological formation
US7447114B2 (en) Non-linear seismic trace matching to well logs
CN105626056B (en) A kind of method and apparatus of determining target area three-dimensional formation pore pressure
US20160320507A1 (en) Time lapse seismic data processing
CN105425298B (en) The method and apparatus of numerical solidification during a kind of elimination finite difference forward modeling
CN112415591A (en) Diffracted wave imaging method and device, electronic equipment and storage medium
US20130158878A1 (en) Device and method for locating microseismic events using array of receivers
CN111337977B (en) Quality factor determination method and device
CN110824564A (en) Attenuation curve chromatographic stripping method for near-surface quality factor Q value inversion
US7382683B1 (en) Computing an absorption parameter for a mode-converted seismic wave
CN105572735B (en) A kind of method and device for improving pre-stack depth migration imaging precision
CN111435174B (en) Method and device for compensating amplitude of seismic data in strong reflection area
CN112711071B (en) Stratum inclination angle correction method and device
CN103454679B (en) A kind of determination method and system of seismic reflection lineups tautochronism
CN109459790B (en) Method and system for establishing seismic velocity field of coal measure strata
CN115598700A (en) Seismic profile imaging method and device, storage medium and electronic equipment
CN111435172A (en) Method and device for chromatographic static correction
CN111999767B (en) Offset imaging method and device for undulating surface
CN116125535B (en) Three-dimensional VSP imaging method and device
CN114397699A (en) Undulating surface Q-migration seismic imaging method and device and terminal equipment
CN116224430A (en) Method and device for calculating dip angle of longitudinal and transverse wave combined stratum
CN112305589A (en) Method and device for imaging depth domain of anisotropic medium

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