CN108919344A - Frequency dividing suitable for layered medium is configured inversion method - Google Patents

Frequency dividing suitable for layered medium is configured inversion method Download PDF

Info

Publication number
CN108919344A
CN108919344A CN201810293307.3A CN201810293307A CN108919344A CN 108919344 A CN108919344 A CN 108919344A CN 201810293307 A CN201810293307 A CN 201810293307A CN 108919344 A CN108919344 A CN 108919344A
Authority
CN
China
Prior art keywords
frequency
well
pseudo
hyperfrequency
medium
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
CN201810293307.3A
Other languages
Chinese (zh)
Other versions
CN108919344B (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.)
BEIJING NUOKE SIDA PETROLEUM TECHNOLOGY CO LTD
Original Assignee
BEIJING NUOKE SIDA PETROLEUM TECHNOLOGY 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 BEIJING NUOKE SIDA PETROLEUM TECHNOLOGY CO LTD filed Critical BEIJING NUOKE SIDA PETROLEUM TECHNOLOGY CO LTD
Priority to CN201810293307.3A priority Critical patent/CN108919344B/en
Publication of CN108919344A publication Critical patent/CN108919344A/en
Application granted granted Critical
Publication of CN108919344B publication Critical patent/CN108919344B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis

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 kind of frequency dividings suitable for layered medium to be configured inversion method, including:Low-frequency range impedance body is established using Kriging regression; intermediate frequency impedance body is converted by seismic data using colored inversion; utilize the configuration modeling controlled based on seismic waveform similarity and Method of Stochastic; it is configured inversion method; high-frequency resistance body is established, low frequency, intermediate frequency, high-frequency resistance modeling or the method for inverting are isolated out, then by three parts impedance body integrator; high-resolution impedance body is obtained, realizes thin-inter bed reservoir prediction.According to the technical solution of the present invention, randomness and unstability in Interpolation Process can be caused less due to well by efficiently solving high band, improved the stability and resolution ratio of inversion result, improved efficiency of inverse process.

Description

Frequency dividing suitable for layered medium is configured inversion method
Technical field
The present invention relates to seismology technical fields more particularly to a kind of frequency dividing suitable for layered medium to be configured inverting side Method.
Background technique
Geostatistical inversion principle, low frequency part obtain low frequency variogram by well logging low-frequency data statistical variation, The conventional determinated back analysis method (Sparse Pulse Inversion) of intermediate-frequency section application obtains wave resistance antibody, to understand substantially dividing for reservoir Cloth, and for seeking intermediate frequency variogram.High frequency section obtains high frequency from well point and is deteriorated, and is closed by three parts variogram At, by stochastic simulation generate well between wave impedance, by the wave impedance of Whole frequency band be converted into reflection coefficient and with determinated back analysis side The wavelet that method acquires carries out convolution and generates synthetic seismogram, by iterating until synthetic seismogram reaches with seismic traces A degree of matching, inverting obtain multiple equiprobable Wave Impedance Data Volumes.
Conventional geostatistics are there are problem, and first, counting reasonable variogram need to number and distribution to fixed well It is required that stringent, it is desirable that well number is more and is evenly distributed, fast for cross directional variations or variation of lithological complex area precision of prediction is inclined It is low.Second, for high frequency stochastic simulation part, conventional geological statistics need stochastic simulation to be up to several hundred frequency range ( High cut-off frequency is shaken less than 100Hz).100Hz frequencies above variation of information Function Fitting difficulty is big, and effect is undesirable, inversion result with Machine is strong, and computational efficiency is low.
Summary of the invention
At least one of regarding to the issue above, the present invention provides a kind of frequency dividings suitable for layered medium to be configured inverting Method, using the thought of frequency dividing building, by it is low, in, the method for middle high, four part impedance modelings of hyperfrequency or inverting isolates It opens, each part is calculated using different modelings and inversion method, improves the stabilization of inversion result, obtains high-resolution The impedance body of rate realizes thin-inter bed reservoir prediction.The random of variogram statistics is replaced using seismic waveform similarity degree interpolation Simulation enhances the stability of a part of result in high-frequency information, has more characteristics of seismic, continuity in inversion result plane It is better than conventional geostatistics efficiency of inverse process.
To achieve the above object, the present invention provides a kind of frequency dividings suitable for layered medium to be configured inversion method, including: Step S101 is filtered the well curve of true well to obtain low-frequency impedance, counts low frequency variogram, and utilize Kriging regression Method carries out interpolation modeling, obtains low-frequency impedance body;Step S102 obtains intermediate frequency impedance body using colored inversion method, and will be described Intermediate frequency impedance body and the low-frequency impedance body synthesize absolute wave impedance body;Step S103 extracts the well bypass road of the true well Well impedance curve simultaneously filters acquisition high-frequency resistance to form Jing Ji, and the well bypass road and serial number of pseudo- well are extracted in seismic data grid Generate the random walk that inverting is carried out to pseudo- well;Step S104, according to the random walk sequence and using the Jing Ji into Row interpolation modeling obtains the medium-high frequency information to be inserted for putting pseudo- well;Step S105, according to the random walk sequence to the pseudo- well The stochastic simulation of hyperfrequency is carried out, Bayesian Decision rule is established with the minimum constraint of composite traces error residual error, to use Multiple equiprobability results that MCMC algorithm carries out optimizing acquisition are averaged to obtain hyperfrequency information;Step S106, with described The Jing Ji will be added after the medium-high frequency information and the synthesis of hyperfrequency information of all pseudo- wells in the sequence of random walk;Step S107 carries out inverting to the remaining road in the seismic data grid with step S104 and step S105 method;Step S108, with The impedance information of the absolute wave impedance body and impedance information of high-frequency resistance matches respectively with true well corresponding frequency band and group is combined into Whole frequency band impedance carries out the seismic data grid to form high-resolution inversion body by road iteration;Wherein, the Frequency Range is 0≤f < 10Hz, and the intermediate frequency is 10≤f < 80Hz, and the high frequency includes the medium-high frequency and the hyperfrequency, institute Stating medium-high frequency is 80≤f < 200Hz, and the hyperfrequency is 200≤f < 500Hz.
In the above-mentioned technical solutions, it is preferable that in the step S104, utilize the well bypass road to be inserted for putting pseudo- well and described true The well bypass road similarity and distance of well are weight, are modeled using the interpolation that the Jing Ji carries out medium-high frequency, and acquisition is to be inserted to put pseudo- well Medium-high frequency information.
In the above-mentioned technical solutions, it is preferable that according to the random walk sequence to the pseudo- well in the step S105 The stochastic simulation for carrying out hyperfrequency specifically includes:It is successively simulated by the random walk sequence, is counted and logged well by probability analysis The mean value and variance of hyperfrequency instruct the stochastic simulation of hyperfrequency in data, while being referred to using earthquake plane variate difference analysis Lead the simulation that the plane of high frequency extends.
In the above-mentioned technical solutions, it is preferable that in the step S106, will first be generated with the sequence of the random walk The medium-high frequency information and hyperfrequency information of pseudo- well synthesize and are added to the well and concentrates, by path points all in the random walk Iterative inversion is completed, to complete the expansion of the Jing Ji.
Compared with prior art, beneficial effects of the present invention are:Using the thought of frequency dividing building, by it is low, in, it is middle high, super The method of four part impedance modelings of high frequency or inverting is isolated out, each part is carried out using different modelings and inversion method It calculates, improves the stabilization of inversion result, obtain high-resolution impedance body, realize thin-inter bed reservoir prediction.Utilize seismic wave Shape similarity degree interpolation replaces the stochastic simulation of variogram statistics, enhances the stability of a part of result in high-frequency information, Characteristics of seismic is had more in inversion result plane, continuity is better than conventional geostatistics efficiency of inverse process.
Detailed description of the invention
Fig. 1 is the flow chart of the disclosed frequency dividing configuration inversion method suitable for layered medium of an embodiment of the present invention;
Fig. 2 is the process signal of the disclosed frequency dividing configuration inversion method suitable for layered medium of an embodiment of the present invention Block diagram;
Fig. 3 is that the frequency spectrum of the disclosed frequency dividing configuration inverting of an embodiment of the present invention constructs form schematic diagram.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is A part of the embodiments of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people Member's every other embodiment obtained without making creative work, shall fall within the protection scope of the present invention.
The present invention is described in further detail with reference to the accompanying drawing:
As shown in Figure 1 to Figure 3, a kind of frequency dividing suitable for layered medium provided according to the present invention is configured inversion method, Including:Step S101 carries out bandpass filtering to the well curve of true well and obtains low-frequency impedance (0~10Hz), counts low frequency variation letter Number, and interpolation modeling is carried out using Kriging regression method, obtain low-frequency impedance body;Step S102, determinated back analysis part utilize Colored inversion method obtains intermediate frequency impedance body (10~80Hz), and intermediate frequency impedance body and low-frequency impedance body are synthesized absolute wave impedance Body (frequency range 0Hz~80Hz);Step S103, data preparation extract the well impedance curve of the well bypass road of true well and filter and obtain High-frequency resistance (80~500Hz of frequency range) is obtained to form Jing Ji, by the well that pseudo- well is extracted in seismic data grid (50*50) Road and serial number generate the random walk that inverting is carried out to pseudo- well;Step S104, pseudo- well are configured interpolation, i.e. the interpolation of medium-high frequency is built Mould, implementation method are:It is weight using the well bypass road similarity and distance of the well bypass road to be inserted for putting pseudo- well and true well, utilizes Jing Ji The interpolation modeling for carrying out medium-high frequency (80~200Hz of frequency range) obtains the medium-high frequency information to be inserted for putting pseudo- well;Step S105, The stochastic simulation of pseudo- well hyperfrequency (200~500Hz of frequency range):According to random walk sequence to pseudo- well carry out hyperfrequency with Machine simulation establishes Bayesian Decision rule, using MCMC (Markov Chain with the minimum constraint of composite traces error residual error Monte Carlo, Markov chain Monte-Carlo) algorithm progress optimizing, it obtains multiple equiprobability results and is averaged, surpassed High-frequency information;After the medium-high frequency information of pseudo- well and the synthesis of hyperfrequency information Jing Ji is added, with random walk in step S106 Sequentially, the pseudo- well first generated is added to the Jing Ji of true well, influences the interpolation result of next path point puppet well, all path points are anti- After the completion of drilling, Jing Ji expansion is finished;Step S107, with step S104 and step S105 method to the residue in seismic data grid Road carries out inverting;Step S108, with the resistance of the impedance information of 0~80Hz absolute wave impedance body and the high-frequency resistance of 80~500Hz Anti- information is matched with the variance of true well corresponding frequency band and mean value respectively, and group is combined into the Whole frequency band impedance of 0~500Hz, right Seismic data grid carries out forming high-resolution inversion body by road iteration;Wherein, Frequency range is 0≤f < 10Hz, in Frequency is 10≤f < 80Hz, and high frequency includes medium-high frequency and hyperfrequency, and medium-high frequency is 80≤f < 200Hz, and hyperfrequency is 200≤f < 500Hz。
In the above-described embodiments, the thought for dividing statistical modeling and weighted band_wise, is different from conventional geostatistics The thought for dividing statistical variation, integrator inverting, is embodied as:Low frequency is established using Kriging regression in step s101 Section impedance body converts intermediate frequency impedance body for seismic data using colored inversion in step S102, and step S103-S107 is utilized Earthquake configuration inversion method establishes high-frequency resistance body, and three parts impedance body integrator is obtained high-resolution by final step S108 The impedance body of rate realizes thin-inter bed reservoir prediction.
Wherein, in the inverting of high frequency configuration, step S103, the pseudo- well building in S104 and S105, it is known that well and point to be inserted away from From smaller, influence degree is bigger, and well distribution consistency degree also determines the stability of surrounding interpolation point, and the building of pseudo- well ensure that Primary data is all uniformly distributed in seismic volume whole network line, efficiently solve high band can be caused less due to well in Interpolation Process with Machine and unstability.High frequency be configured inverting in be configured interpolation thought, basis be think lithology, the variation of physical property will lead to The variation of seismic data frequency, phase, amplitude, the form of waveform and the combination of reservoir are related, and in other words waveform is similar, surveys Well curve is similar, and seismic waveform similarity degree can determine whether log extrapolation degree.Therefore using the phase of seismic waveform structure The extrapolation of well logging high-frequency data is controlled like degree.Described in step S104, using seismic waveform similarity and to anti- It drills and a little carries out interpolation method at a distance from well point to calculate weight, utilize intermediate frequency certainty anti-instead of in geostatistical inversion Variation and the high frequency well logging drilled, which are deteriorated, calculates weight progress interpolation method, avoids high frequency well logging variogram and calculates randomness By force, the problem of plane extrapolability difference.
Specifically, frequency dividing configuration inverting is according to the thought of frequency dividing building, by it is low, in, middle high, four part impedances of hyperfrequency The method of modeling or inverting is isolated out, each part is calculated using different modelings and inversion method, and single inverting Method is not wide in range to data relevance grade, be easy to cause the unstable of inversion result.Low-frequency stability is strong, utilizes Kriging regression solution Certainly;Intermediate frequency information can stablize acquisition in the seismic data, can effectively keep the intermediate frequency of earthquake to believe using colored inversion method Breath;Medium-high frequency is effectively extrapolated using seismic data configural characteristic and carries out interpolation modeling;Hyperfrequency information randomness is strong, nothing in earthquake Method obtains, and is obtained using hyperfrequency stochastic simulation inversion method in log data.Medium-high frequency is configured in modeling method, utilizes earthquake Waveform similarity degree interpolation replaces the stochastic simulation of variogram statistics, enhances the stabilization of a part of result in high-frequency information Property, characteristics of seismic is had more in inversion result plane, continuity is better than conventional geostatistics efficiency of inverse process.
The above is embodiments of the present invention, and the frequency dividing configuration suitable for layered medium proposed according to the present invention is anti- Drill method, using the thought of frequency dividing building, by it is low, in, the method for middle high, four part impedance modelings of hyperfrequency or inverting isolates It opens, each part is calculated using different modelings and inversion method, improves the stabilization of inversion result, obtains high-resolution The impedance body of rate realizes thin-inter bed reservoir prediction.The random of variogram statistics is replaced using seismic waveform similarity degree interpolation Simulation enhances the stability of a part of result in high-frequency information, has more characteristics of seismic, continuity in inversion result plane It is better than conventional geostatistics efficiency of inverse process.
These are only the preferred embodiment of the present invention, is not intended to restrict the invention, for those skilled in the art For member, the invention may be variously modified and varied.All within the spirits and principles of the present invention, it is made it is any modification, Equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.

Claims (4)

1. a kind of frequency dividing suitable for layered medium is configured inversion method, which is characterized in that including:
Step S101 is filtered the well curve of true well to obtain low-frequency impedance, counts low frequency variogram, and utilize Ke Lijin Interpolation method carries out interpolation modeling, obtains low-frequency impedance body;
Step S102 obtains intermediate frequency impedance body using colored inversion method, and by the intermediate frequency impedance body and the low-frequency impedance body Synthesize absolute wave impedance body;
Step S103, extracts the well impedance curve of the well bypass road of the true well and filters and obtain high-frequency resistance to form Jing Ji, Well bypass road and serial number gencration that pseudo- well is extracted in seismic data grid carry out the random walk of inverting to pseudo- well;
Step S104 according to the sequence of the random walk and carries out interpolation modeling using the Jing Ji, obtains the pseudo- well of point to be inserted Medium-high frequency information;
Step S105 is carried out the stochastic simulation of hyperfrequency to the pseudo- well according to the random walk sequence, is missed with composite traces Bayesian Decision rule is established in the poor minimum constraint of residual error, to the multiple equiprobability results for carrying out optimizing acquisition using MCMC algorithm It is averaged to obtain hyperfrequency information;
Step S106 is synthesized the medium-high frequency information of all pseudo- wells and hyperfrequency information with the sequence of the random walk After the Jing Ji is added;
Step S107 carries out inverting to the remaining road in the seismic data grid with step S104 and step S105 method;
Step S108, with the impedance information of the absolute wave impedance body and the impedance information of high-frequency resistance frequency corresponding with true well respectively Section matches and group is combined into Whole frequency band impedance, carries out forming high-resolution inversion body by road iteration to the seismic data grid;
Wherein, the Frequency range is 0≤f < 10Hz, and the intermediate frequency is 10≤f < 80Hz, and the high frequency includes described Medium-high frequency and the hyperfrequency, the medium-high frequency are 80≤f < 200Hz, and the hyperfrequency is 200≤f < 500Hz.
2. the frequency dividing according to claim 1 suitable for layered medium is configured inversion method, which is characterized in that the step In S104, it is weight using the well bypass road similarity and distance of the well bypass road to be inserted for putting pseudo- well and the true well, utilizes the well Collection carries out the interpolation modeling of medium-high frequency, obtains the medium-high frequency information to be inserted for putting pseudo- well.
3. the frequency dividing according to claim 1 suitable for layered medium is configured inversion method, which is characterized in that the step It is specifically included in S105 according to the stochastic simulation that the random walk sequence carries out hyperfrequency to the pseudo- well:
It is successively simulated by the random walk sequence, the mean value and variance pair of hyperfrequency in log data is counted by probability analysis The stochastic simulation of hyperfrequency is instructed, while passing through the plane variate difference analysis of seismic data, and radio-frequency component plane is instructed to extend Simulation.
4. the frequency dividing according to claim 1 suitable for layered medium is configured inversion method, which is characterized in that the step In S106, the medium-high frequency information of the pseudo- well first generated and hyperfrequency information are synthesized and is added to the sequence of the random walk The well is concentrated, and path point iterative inversions all in the random walk is completed, to complete the expansion of the Jing Ji.
CN201810293307.3A 2018-03-30 2018-03-30 Frequency division configuration inversion method suitable for layered medium Active CN108919344B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810293307.3A CN108919344B (en) 2018-03-30 2018-03-30 Frequency division configuration inversion method suitable for layered medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810293307.3A CN108919344B (en) 2018-03-30 2018-03-30 Frequency division configuration inversion method suitable for layered medium

Publications (2)

Publication Number Publication Date
CN108919344A true CN108919344A (en) 2018-11-30
CN108919344B CN108919344B (en) 2020-04-21

Family

ID=64403400

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810293307.3A Active CN108919344B (en) 2018-03-30 2018-03-30 Frequency division configuration inversion method suitable for layered medium

Country Status (1)

Country Link
CN (1) CN108919344B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112444864A (en) * 2019-08-28 2021-03-05 中国石油化工股份有限公司 Thin-layer high-resolution inversion method and storage medium
CN112784379A (en) * 2021-03-09 2021-05-11 重庆邮电大学 Optimization design method of impeller rotor system
CN113917529A (en) * 2020-07-07 2022-01-11 中国石油化工股份有限公司 Wave impedance inversion method, wave impedance inversion device, storage medium and electronic equipment

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104570066A (en) * 2013-10-10 2015-04-29 中国石油天然气股份有限公司 Construction Method of Seismic Inversion Low-Frequency Model
CN105425289A (en) * 2015-10-29 2016-03-23 中国石油天然气集团公司 Method and device of determining low frequency wave impedance
WO2016044538A1 (en) * 2014-09-19 2016-03-24 Conocophillips Company Bandwidth extension beyond the vibrator sweep signal via a constrained simultaneous multiple vibrator inversion
CN105891888A (en) * 2016-03-28 2016-08-24 吉林大学 Multi-domain frequency-division parallel multi-scale full-waveform inversion method
CN106707339A (en) * 2016-11-22 2017-05-24 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Depth domain phase-control inversion method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104570066A (en) * 2013-10-10 2015-04-29 中国石油天然气股份有限公司 Construction Method of Seismic Inversion Low-Frequency Model
WO2016044538A1 (en) * 2014-09-19 2016-03-24 Conocophillips Company Bandwidth extension beyond the vibrator sweep signal via a constrained simultaneous multiple vibrator inversion
CN105425289A (en) * 2015-10-29 2016-03-23 中国石油天然气集团公司 Method and device of determining low frequency wave impedance
CN105891888A (en) * 2016-03-28 2016-08-24 吉林大学 Multi-domain frequency-division parallel multi-scale full-waveform inversion method
CN106707339A (en) * 2016-11-22 2017-05-24 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Depth domain phase-control inversion method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
于建国,等: "分频反演方法及应用", 《石油地球物理勘探》 *
张繁昌,等: "叠前地震道集AVA多频信息同时反演", 《石油物探》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112444864A (en) * 2019-08-28 2021-03-05 中国石油化工股份有限公司 Thin-layer high-resolution inversion method and storage medium
CN113917529A (en) * 2020-07-07 2022-01-11 中国石油化工股份有限公司 Wave impedance inversion method, wave impedance inversion device, storage medium and electronic equipment
CN112784379A (en) * 2021-03-09 2021-05-11 重庆邮电大学 Optimization design method of impeller rotor system
CN112784379B (en) * 2021-03-09 2023-05-23 重庆邮电大学 Optimal design method of impeller rotor system

Also Published As

Publication number Publication date
CN108919344B (en) 2020-04-21

Similar Documents

Publication Publication Date Title
CN105182444B (en) A kind of Inversion of geostatistics of High-Resolution Sequence Stratigraphic Framework constraint
CN108919344A (en) Frequency dividing suitable for layered medium is configured inversion method
CN106597532A (en) Pre-stack seismic data frequency band expanding method of combining well information and horizon information
CN108267784A (en) A kind of seismic signal random noise compression process method
CN104820243B (en) Three-dimensional post-stack seismic data amplitude compensation method and device
CN107272062A (en) A kind of underground medium Q methods of estimation of data-driven
CN109425900B (en) Seismic reservoir prediction method
CN104122581B (en) A kind of poststack sound impedance inversion method
CN108254731B (en) The multiple dimensioned staged layer of Coherent Noise in GPR Record removes full waveform inversion method
CN107861149B (en) Based on the prestack P-S wave velocity ratio analogy method under drive waveform
CN105301647B (en) The method for distinguishing grey matter mud stone and sandstone
CN104991272A (en) Earthquake speed disturbance modeling method for well-free earthquake reversion
CN108931811A (en) Seismic Reservoir inversion method based on multiple spot geological statistics
CN106033125A (en) Suppression prestack large-angle gather interference based frequency increasing method
CN105425286A (en) Earthquake time-travelling acquisition method and crosshole earthquake time-travelling tomography method based on the earthquake time-travelling acquisition method
CN107179550A (en) A kind of seismic signal zero phase deconvolution method of data-driven
CN106873036A (en) Denoising method based on well-seismic combination
CN111751877A (en) Seismic data multiple-integration coherent fracture prediction method and device
CN104297792B (en) The phased inversion method of water channel reservoir is stacked on a kind of fan
CN105954803A (en) Post-stack seismic inversion method and device
CN110095811A (en) Method and device for constructing and processing velocity model of paste rock stratum
CN103217715B (en) Multiple dimensioned regular grid Static Correction of Tomographic Inversion method
CN105487115A (en) Wavelet transform-based high frequency continuation method
CN106842291B (en) Unconformity trapped reservoir lithology prediction method based on prestack seismic ray impedance inversion
CN109884699A (en) Beach dam method for predicting reservoir based on super complete wavelet base adjacent objects function method

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