CN116338786A - Earthquake inversion method and system for eliminating strong reflection shielding effect - Google Patents

Earthquake inversion method and system for eliminating strong reflection shielding effect Download PDF

Info

Publication number
CN116338786A
CN116338786A CN202111599576.0A CN202111599576A CN116338786A CN 116338786 A CN116338786 A CN 116338786A CN 202111599576 A CN202111599576 A CN 202111599576A CN 116338786 A CN116338786 A CN 116338786A
Authority
CN
China
Prior art keywords
wave impedance
pseudo
seismic
inversion
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
Application number
CN202111599576.0A
Other languages
Chinese (zh)
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 CN202111599576.0A priority Critical patent/CN116338786A/en
Publication of CN116338786A publication Critical patent/CN116338786A/en
Pending legal-status Critical Current

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/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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • 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/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616—Data from specific type of measurement
    • G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
    • 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/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616—Data from specific type of measurement
    • G01V2210/6169—Data from specific type of measurement using well-logging
    • 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
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6226—Impedance
    • 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
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (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 seismic inversion method for eliminating a strong reflection shielding effect, which comprises the following steps: acquiring post-stack seismic data and a logging wave impedance curve, and performing well shock calibration; obtaining unshielded seismic data based on the seismic data; obtaining a pseudo-wave impedance curve based on the log-wave impedance curve, and matching the synthetic seismic record of the pseudo-wave impedance curve with the unshielded seismic data; acquiring pseudo-wave impedance inversion data based on the unshielded seismic data and the pseudo-wave impedance curve; and developing reservoir comprehensive interpretation based on the pseudo-wave impedance inversion data. The method specifically eliminates and weakens the influence of strong reflection energy shielding and Gibbs effect, and can improve the resolution of seismic inversion and the accuracy of reservoir prediction.

Description

Earthquake inversion method and system for eliminating strong reflection shielding effect
Technical Field
The invention belongs to the technical field of interpretation of seismic data of petroleum geophysical exploration, and particularly relates to a seismic inversion method and system for eliminating a strong reflection shielding effect.
Background
Formations such as coal seams, source rocks, volcanic rocks, etc. that have a large difference in impedance from surrounding rocks often develop in the subsurface formation to form strong reflections on the seismic profile. The stronger the seismic reflection energy, the more the wave impedance between the formations varies. Discontinuities are exhibited when the wave impedance of the formation varies too strongly across the strongly reflective interface.
The vicinity of such discontinuous boundaries is prone to Gibbs effect (Gibbs phenomenon) limited by the seismic traces and wavelet lengths, the seismic sampling rate, and the seismic inversion frequency band. The Gibbs effect causes severe oscillations in the seismic inversion results near the strongly reflecting interface. Meanwhile, the wave impedance difference between the rock formations is very small in the upper and lower strata of the strong reflection interface, the wave impedance difference is far smaller than oscillation generated by Gibbs phenomenon, a serious shielding effect is generated, the inversion result is distorted, the real wave impedance characteristics of the strata are difficult to accurately reflect, and the accuracy of reservoir prediction is affected.
In the published literature on deterministic seismic inversion methods, 2 methods mainly applied in practice include model-based inversion, constraint sparse pulse inversion and the like. The method is limited by the realization thought of the conventional inversion method and the limited bandwidth of the seismic data, the geological features of discontinuous wave impedance are not considered, and accurate prediction of reservoirs near a strong reflection interface is difficult to achieve. Thus, no effective solution has been currently focused on and proposed for seismic inversion in such geological situations.
Therefore, a need exists for a seismic inversion method that eliminates the strong reflection shielding effect to overcome the shortcomings of current seismic inversion methods.
Disclosure of Invention
In order to further reduce the shielding effect of the strong reflection interface on the two side wave impedances in the seismic inversion, the seismic inversion method needs to be further optimized, and errors generated by the Gibbs effect are eliminated and weakened, so that the wave impedance characteristics of the reservoirs at the two sides of the strong reflection interface are better kept, and the accuracy of the seismic inversion prediction of the reservoirs is improved.
In view of the above, the present invention provides a seismic inversion method that eliminates the strong reflection shielding effect,
the method comprises the following steps:
acquiring post-stack seismic data and a logging wave impedance curve, and performing well shock calibration;
obtaining unshielded seismic data based on the seismic data;
obtaining a pseudo-wave impedance curve based on the log-wave impedance curve, and matching the synthetic seismic record of the pseudo-wave impedance curve with the unshielded seismic data;
acquiring pseudo-wave impedance inversion data based on the unshielded seismic data and the pseudo-wave impedance curve;
and developing reservoir comprehensive interpretation based on the pseudo-wave impedance inversion data.
Further, the obtaining the unshielded seismic data is specifically:
performing post-stack unshielded processing on the seismic data to obtain unshielded seismic data;
the post-stack de-masking treatment of the seismic data is to adopt a reflection coefficient fitting method to realize the de-strong reflection masking treatment of the seismic data.
Further, the obtaining a pseudo-wave impedance curve specifically includes:
and performing step-removing treatment on the logging wave impedance curve to obtain a pseudo wave impedance curve, so that the synthetic seismic record of the pseudo wave impedance curve is matched with the unshielded seismic data.
Further, the acquiring pseudo-wave impedance inversion data includes:
establishing a wave impedance initial model by utilizing interpolation of a pseudo-wave impedance curve;
based on the unshielded seismic data and the pseudo-wave impedance curve, extracting wavelets, performing post-stack seismic inversion, and obtaining pseudo-wave impedance inversion data.
Further, developing the reservoir comprehensive interpretation based on the pseudo-wave impedance inversion data includes:
determining a value range of the pseudo-wave impedance of the reservoir based on histogram analysis of the pseudo-wave impedance curve;
and carrying out reservoir identification on the pseudo-wave impedance inversion data according to the range of the pseudo-wave impedance of the reservoir.
The present invention also provides a seismic inversion system for eliminating strong reflection shielding effects, the system comprising:
the acquisition unit is used for acquiring post-stack seismic data and a logging wave impedance curve and carrying out well seismic calibration;
obtaining unshielded seismic data based on the seismic data;
obtaining a pseudo-wave impedance curve based on the log-wave impedance curve, and matching the synthetic seismic record of the pseudo-wave impedance curve with the unshielded seismic data;
acquiring pseudo-wave impedance inversion data based on the unshielded seismic data and the pseudo-wave impedance curve;
and the analysis and interpretation unit is used for analyzing inversion data based on the pseudo wave impedance and developing reservoir comprehensive interpretation.
Further, the obtaining the unshielded seismic data is specifically:
performing post-stack unshielded processing on the seismic data to obtain unshielded seismic data;
the post-stack de-masking treatment of the seismic data is to adopt a reflection coefficient fitting method to realize the de-strong reflection masking treatment of the seismic data.
Further, the obtaining a pseudo-wave impedance curve specifically includes:
and performing step-removing treatment on the logging wave impedance curve to obtain a pseudo wave impedance curve, so that the synthetic seismic record of the pseudo wave impedance curve is matched with the unshielded seismic data.
Further, the acquiring pseudo-wave impedance inversion data includes:
establishing a wave impedance initial model by utilizing interpolation of a pseudo-wave impedance curve;
based on the unshielded seismic data and the pseudo-wave impedance curve, extracting wavelets, performing post-stack seismic inversion, and obtaining pseudo-wave impedance inversion data.
Further, developing the reservoir comprehensive interpretation based on the pseudo-wave impedance inversion data includes:
determining a value range of the pseudo-wave impedance of the reservoir based on histogram analysis of the pseudo-wave impedance curve;
and carrying out reservoir identification on the pseudo-wave impedance inversion data according to the range of the pseudo-wave impedance of the reservoir.
The invention has the following technical effects:
(1) Aiming at the problem of seismic inversion numerical distortion near a strong reflection interface, different conventional seismic inversion technical ideas are adopted. Eliminating the influence of strong reflection energy shielding through seismic processing, and recovering and highlighting thin-layer weak reflection seismic signals; and performing step removal processing on the well logging data to eliminate the Gibbs effect, so that the seismic inversion result can better reflect the geophysical characteristics of the thin layer, and the reservoir prediction precision is improved.
(2) The technology is relatively independent, the processed seismic data and pseudo-wave impedance curves are utilized, the seismic inversion is realized through the existing inversion software, and the seismic data well seismic calibration, interpolation modeling technology, wavelet extraction technology and seismic inversion technology are adopted in the whole scheme, so that the technology feasibility is good.
(3) Compared with the current seismic attribute prediction method after the unmasking treatment, the method provided by the invention can realize quantitative reservoir prediction results.
(4) Compared with the traditional seismic inversion method and flow, the method provided by the invention can purposefully eliminate and weaken the influence of strong reflection energy shielding and Gibbs effect, and can improve the resolution of seismic inversion and the accuracy of reservoir prediction.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention may be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 illustrates a flow chart of a seismic inversion method of eliminating strong reflection shielding effects in accordance with an embodiment of the invention;
FIG. 2 shows the results of well A and measured wave impedance curve borehole shock calibration in accordance with an embodiment of the present invention;
FIG. 3 illustrates a pre-reflection-shielding stacked seismic section according to an embodiment of the invention;
FIG. 4 illustrates a post-de-reflection-shielding stacked seismic section in accordance with an embodiment of the invention;
FIG. 5 illustrates the well A pseudowave impedance curve and the results of unshielded seismic borehole seismic calibration in accordance with an embodiment of the invention;
FIG. 6 shows a pseudo-wave impedance inversion profile of well A according to an embodiment of the invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention provides a seismic inversion method for eliminating a strong reflection shielding effect, as shown in fig. 1, and fig. 1 shows a flow chart of the seismic inversion method for eliminating the strong reflection shielding effect according to an embodiment of the invention. The method comprises the following steps: acquiring post-stack seismic data and a logging wave impedance curve, and performing well shock calibration; performing post-stack unshielded processing on the seismic data to obtain unshielded seismic data; performing step-removing treatment on the logging curve to obtain a pseudo-wave impedance curve, so that the synthetic seismic record of the pseudo-wave impedance curve is matched with the unshielded seismic data; establishing a wave impedance initial model by utilizing interpolation of a pseudo-wave impedance curve; extracting wavelets based on the unshielded seismic data and the pseudo-wave impedance curve, and performing post-stack seismic inversion to obtain pseudo-wave impedance inversion data; determining a value range of the pseudo-wave impedance of the reservoir based on histogram analysis of the pseudo-wave impedance curve; and carrying out reservoir identification on the pseudo-wave impedance inversion data according to the range of the pseudo-wave impedance of the reservoir.
Specifically, acquiring post-stack seismic data and a logging wave impedance curve, and performing well seismic calibration; FIG. 2 shows the results of well shock calibration performed on post-stack seismic data and well A log impedance curves for an oilfield. The location shown by the strong reflection interface in FIG. 2a is a strong energy seismic reflection event in the region, the strong reflection energy being much higher than the lower formation reflection; FIG. 2b shows a sudden change in the nearby wave impedance curve corresponding to the strongly reflecting interface, with the reservoir being within 20ms of the strongly reflecting interface; FIG. 2c is a synthetic seismic log, where the corresponding weak peak seismic reflection of the reservoir is masked by wavelet sidelobes of the strong reflection interface (the position indicated by the arrow in the figure) due to weak energy, and it is difficult to truly reflect the transverse variation characteristics of the reservoir.
Performing post-stack unshielded processing on the seismic data to obtain unshielded seismic data; in specific implementation, the reflection coefficient fitting method is adopted to realize the anti-strong reflection shielding treatment of the seismic data. To meet the inversion requirement, 2 conditions are required to be satisfied simultaneously: the seismic response characteristics of other formations are maintained while removing the strong reflected energy generated by the discontinuous interface.
In specific implementation, according to the condition of the stratum nearby the strong reflection interface, the intensity of the strong reflection axis in the seismic data and the quality of the seismic data, under the condition that the seismic response characteristics of other stratum are maintained while the strong reflection energy is removed, the post-stack strong reflection removal processing method is not limited to reflection coefficient fitting, and methods such as matching tracking and lithology removal processing can be adopted.
Fig. 3 and fig. 4 are respectively a seismic section before and after the reflection coefficient fitting method is adopted to perform the strong reflection shielding treatment, the seismic reflection amplitude of the reservoir section at the lower part of the strong reflection interface is strengthened and the continuity is improved before the comparison treatment, and the seismic reflection energy of the reservoir is effectively recovered, which is indicated by an arrow in fig. 4.
Step-removing the logging curve to obtain a pseudo wave impedance curve matched with the unshielded post-stack seismic data;
the seismic data is subjected to strong reflection removal shielding treatment, which is equivalent to eliminating the wave impedance difference of stratum at two sides of a strong reflection interface, and the energy is greatly changed. In order to match the de-screened post-stack seismic data to perform seismic inversion, the wave impedance difference corresponding to the strong reflection needs to be eliminated from the log data to obtain a pseudo wave impedance curve.
In practice, to obtain log data that matches the unshielded post-stack seismic data, a pseudo-wave impedance curve is obtained by the following steps.
(1) Determining the top surface depth t0 and the bottom surface depth t1 of a discontinuous boundary corresponding to the strong reflection interface;
(2) Respectively obtaining average wave impedance of stratum at two sides of the strong reflection interface, and obtaining ladder wave impedance of the corresponding strong reflection interface according to the following formula;
Figure BDA0003432737980000061
wherein,,
Figure BDA0003432737980000062
representing the step wave impedance; PI (proportional integral) u Representing the average wave impedance of the formation above the strongly reflective interface; PI (proportional integral) d Representing the average wave impedance of the formation below the strongly reflective interface; t is t i Representing an i-th depth sampling point;
(3) To maintain the measured wave impedance characteristics of the target interval, calculating the pseudo wave impedance according to the following formula;
Figure BDA0003432737980000063
wherein,,
Figure BDA0003432737980000064
representing the pseudo-wave impedance; />
Figure BDA0003432737980000065
Representing the step wave impedance; PI (proportional integral) i Representing an initial wave impedance; PI (proportional integral) u Representing the average wave impedance of the formation above the strongly reflective interface;
fig. 5b is a pseudo-wave impedance curve of well a obtained using the above procedure. FIG. 5c is a synthetic seismic record corresponding to a pseudo-wave impedance curve, which is very similar to the unshielded seismic data of FIG. 5a, illustrating that the pseudo-wave impedance curve matches well with the actual seismic data while maintaining the reservoir wave impedance characteristics.
Establishing a wave impedance initial model by utilizing interpolation of a pseudo-wave impedance curve; based on the unshielded seismic data and the pseudo-wave impedance curve, extracting wavelets, performing post-stack seismic inversion, and obtaining pseudo-wave impedance inversion data. FIG. 6 is a pseudo-wave impedance inversion profile with a gray scale color depth indicating a relatively low wave impedance value indicative of reservoir development. In the figure, the inversion result of the pseudo-wave impedance accords with a low-value area of the wave impedance curve on the well point A well, namely the development section of the reservoir is facilitated. And further analyzing the pseudo-wave impedance distribution characteristics of the reservoir, and carrying out reservoir comprehensive interpretation based on the pseudo-wave impedance inversion data.
The invention also provides a seismic inversion system for eliminating the strong reflection shielding effect, the system comprises: the acquisition unit is used for acquiring post-stack seismic data and a logging wave impedance curve and carrying out well seismic calibration; obtaining unshielded seismic data based on the seismic data; obtaining a pseudo-wave impedance curve based on the log-wave impedance curve, and matching the synthetic seismic record of the pseudo-wave impedance curve with the unshielded seismic data; based on the unshielded seismic data and the pseudo-wave impedance curve, pseudo-wave impedance inversion data is obtained.
The system also comprises an analysis and interpretation unit which is used for analyzing inversion data based on the pseudo-wave impedance and developing reservoir comprehensive interpretation.
The method for obtaining the unshielded seismic data specifically comprises the following steps: and performing post-stack unshielded processing on the seismic data to obtain unshielded seismic data.
The pseudo-wave impedance curve is obtained specifically as follows: and performing step-removing treatment on the logging wave impedance curve to obtain a pseudo wave impedance curve, so that the synthetic seismic record of the pseudo wave impedance curve is matched with the unshielded seismic data.
Acquiring pseudo-wave impedance inversion data includes: establishing a wave impedance initial model by utilizing interpolation of a pseudo-wave impedance curve; based on the unshielded seismic data and the pseudo-wave impedance curve, extracting wavelets, performing post-stack seismic inversion, and obtaining pseudo-wave impedance inversion data. Based on the pseudo-wave impedance inversion data, developing reservoir comprehensive interpretation includes: determining a value range of the pseudo-wave impedance of the reservoir based on histogram analysis of the pseudo-wave impedance curve; and carrying out reservoir identification on the pseudo-wave impedance inversion data according to the range of the pseudo-wave impedance of the reservoir.
Although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (10)

1. A seismic inversion method of eliminating strong reflection shielding effects, the method comprising:
acquiring post-stack seismic data and a logging wave impedance curve, and performing well shock calibration;
obtaining unshielded seismic data based on the seismic data;
obtaining a pseudo-wave impedance curve based on the log-wave impedance curve, and matching the synthetic seismic record of the pseudo-wave impedance curve with the unshielded seismic data;
acquiring pseudo-wave impedance inversion data based on the unshielded seismic data and the pseudo-wave impedance curve;
and developing reservoir comprehensive interpretation based on the pseudo-wave impedance inversion data.
2. The method of seismic inversion with strong reflection shielding effect cancellation according to claim 1, wherein the obtaining of unshielded seismic data is specifically:
performing post-stack unshielded processing on the seismic data to obtain unshielded seismic data;
the post-stack de-masking treatment of the seismic data is to adopt a reflection coefficient fitting method to realize the de-strong reflection masking treatment of the seismic data.
3. The method of seismic inversion with strong reflection shielding effect cancellation according to claim 1, wherein the obtaining a pseudo-wave impedance curve is specifically:
and performing step-removing treatment on the logging wave impedance curve to obtain a pseudo wave impedance curve, so that the synthetic seismic record of the pseudo wave impedance curve is matched with the unshielded seismic data.
4. A method of seismic inversion to cancel a strong reflection shielding effect according to any of claims 1-3, wherein the obtaining pseudo-wave impedance inversion data comprises:
establishing a wave impedance initial model by utilizing interpolation of a pseudo-wave impedance curve;
based on the unshielded seismic data and the pseudo-wave impedance curve, extracting wavelets, performing post-stack seismic inversion, and obtaining pseudo-wave impedance inversion data.
5. The seismic inversion method of claim 4, wherein developing reservoir integrated interpretation based on the pseudo-wave impedance inversion data comprises:
determining a value range of the pseudo-wave impedance of the reservoir based on histogram analysis of the pseudo-wave impedance curve;
and carrying out reservoir identification on the pseudo-wave impedance inversion data according to the range of the pseudo-wave impedance of the reservoir.
6. A seismic inversion system that eliminates strong reflection shielding effects, the system comprising:
the acquisition unit is used for acquiring post-stack seismic data and a logging wave impedance curve and carrying out well seismic calibration;
obtaining unshielded seismic data based on the seismic data;
obtaining a pseudo-wave impedance curve based on the log-wave impedance curve, and matching the synthetic seismic record of the pseudo-wave impedance curve with the unshielded seismic data;
acquiring pseudo-wave impedance inversion data based on the unshielded seismic data and the pseudo-wave impedance curve;
and the analysis and interpretation unit is used for analyzing inversion data based on the pseudo wave impedance and developing reservoir comprehensive interpretation.
7. The seismic inversion system of claim 6 wherein said obtaining unshielded seismic data is specifically:
performing post-stack unshielded processing on the seismic data to obtain unshielded seismic data;
the post-stack de-masking treatment of the seismic data is to adopt a reflection coefficient fitting method to realize the de-strong reflection masking treatment of the seismic data.
8. The seismic inversion system of claim 6 wherein said obtaining a pseudo-wave impedance curve is specifically:
and performing step-removing treatment on the logging wave impedance curve to obtain a pseudo wave impedance curve, so that the synthetic seismic record of the pseudo wave impedance curve is matched with the unshielded seismic data.
9. The seismic inversion system of claim 6 wherein said obtaining spurious impedance inversion data comprises:
establishing a wave impedance initial model by utilizing interpolation of a pseudo-wave impedance curve;
based on the unshielded seismic data and the pseudo-wave impedance curve, extracting wavelets, performing post-stack seismic inversion, and obtaining pseudo-wave impedance inversion data.
10. The seismic inversion system of any of claims 6-9 wherein developing reservoir integrated interpretation based on the pseudo-wave impedance inversion data comprises:
determining a value range of the pseudo-wave impedance of the reservoir based on histogram analysis of the pseudo-wave impedance curve;
and carrying out reservoir identification on the pseudo-wave impedance inversion data according to the range of the pseudo-wave impedance of the reservoir.
CN202111599576.0A 2021-12-24 2021-12-24 Earthquake inversion method and system for eliminating strong reflection shielding effect Pending CN116338786A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111599576.0A CN116338786A (en) 2021-12-24 2021-12-24 Earthquake inversion method and system for eliminating strong reflection shielding effect

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111599576.0A CN116338786A (en) 2021-12-24 2021-12-24 Earthquake inversion method and system for eliminating strong reflection shielding effect

Publications (1)

Publication Number Publication Date
CN116338786A true CN116338786A (en) 2023-06-27

Family

ID=86879449

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111599576.0A Pending CN116338786A (en) 2021-12-24 2021-12-24 Earthquake inversion method and system for eliminating strong reflection shielding effect

Country Status (1)

Country Link
CN (1) CN116338786A (en)

Similar Documents

Publication Publication Date Title
US8923093B2 (en) Determining the quality of a seismic inversion
CN111722282B (en) Method for predicting natural gas hydrate reservoir top hydrate saturation by AVO
CN114994758B (en) Wave impedance extraction and structure characterization method and system for carbonate rock fracture control reservoir
EA001509B1 (en) A method of processing seismic data signals
CN105093306A (en) Method for automatic interpretation and thickness calculation of reservoir in geophysical exploration
CN108508489B (en) Seismic inversion method based on waveform micro-variation matching
CN116520419B (en) Hot fluid crack channel identification method
CN113109875B (en) Inversion method of carbonate rock reservoir under full waveform velocity field constraint
CN112946782B (en) Earthquake fine depicting method for dense oil-gas storage seepage body
CN108121008A (en) Method for predicting spatial distribution of river sand by seismic attributes
CN109283577B (en) Seismic horizon calibration method
CN114488293A (en) High-resolution inversion method based on sensitive logging curve
CN110244383B (en) Geological lithology comprehensive model establishing method based on near-surface data
CN116338786A (en) Earthquake inversion method and system for eliminating strong reflection shielding effect
CN111880223B (en) Attribute analysis method for well-side sound wave remote detection imaging fracture-cave carving
CN110632660B (en) Thin sand body characterization method and device based on seismic data body
CN113848593A (en) Method for quantitatively predicting rock slurry erosion area in coal-bearing stratum
CN107831537B (en) Method for describing high sand-to-land ratio nepheloid rock reservoir
CN112147676A (en) Method for predicting thickness of coal bed and gangue
CN113806674A (en) Method and device for quantifying longitudinal dimension of ancient river channel, electronic equipment and storage medium
CN112394398A (en) Stratum pinch-out line identification method based on lithofacies and seismic facies
CN116338784A (en) Method and system for inversion of earthquake before lamination of weathering crust karst reservoir
CN113589365A (en) Reservoir pinch-out line description method based on time-frequency domain information
CN113640873B (en) Sandstone uranium ore pre-stack earthquake prediction method and device
US20240069228A1 (en) Learning hydrocarbon distribution from seismic image

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