CN108919344B - Frequency division configuration inversion method suitable for layered medium - Google Patents
Frequency division configuration inversion method suitable for layered medium Download PDFInfo
- Publication number
- CN108919344B CN108919344B CN201810293307.3A CN201810293307A CN108919344B CN 108919344 B CN108919344 B CN 108919344B CN 201810293307 A CN201810293307 A CN 201810293307A CN 108919344 B CN108919344 B CN 108919344B
- Authority
- CN
- China
- Prior art keywords
- frequency
- well
- medium
- impedance
- inversion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000004088 simulation Methods 0.000 claims abstract description 16
- 238000001914 filtration Methods 0.000 claims description 6
- 230000002194 synthesizing effect Effects 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 2
- 230000008569 process Effects 0.000 abstract description 2
- 238000010276 construction Methods 0.000 description 6
- 208000035126 Facies Diseases 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 238000013213 extrapolation Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000002948 stochastic simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
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
-
- 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
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 frequency division configuration inversion method suitable for a laminated medium, which comprises the following steps: the method comprises the steps of establishing a low-frequency-band impedance body by using kriging interpolation, converting seismic data into a medium-frequency impedance body by using colored inversion, establishing a high-frequency impedance body by using a configuration modeling and random simulation method based on seismic waveform similarity control, namely a configuration inversion method, splitting the low-frequency impedance body, the medium-frequency impedance body and the high-frequency impedance modeling or inversion method, and then combining the three impedance bodies into a whole to obtain the high-resolution impedance body so as to realize thin interbed reservoir prediction. Through the technical scheme of the invention, the randomness and instability in the interpolation process caused by few wells in a high-frequency section are effectively solved, the stability and the resolution of the inversion result are improved, and the inversion effect is improved.
Description
Technical Field
The invention relates to the technical field of seismology, in particular to a frequency division configuration inversion method suitable for a laminated medium.
Background
According to the geostatistical inversion principle, a low-frequency part obtains a low-frequency variation function through the statistical variation of logging low-frequency data, and a medium-frequency part obtains a wave impedance body by applying a conventional deterministic inversion method (sparse pulse inversion) so as to know the approximate distribution of a reservoir and is used for solving the medium-frequency variation function. The high-frequency part obtains high-frequency variation from a well point, the high-frequency variation is synthesized through three-part variation function, inter-well wave impedance is generated through random simulation, wave impedance of a full frequency band is converted into a reflection coefficient and is subjected to convolution with wavelets obtained by a deterministic inversion method to generate a synthetic seismic channel, and a plurality of equal-probability wave impedance data volumes are obtained through repeated iteration until the synthetic seismic channel is matched with an original seismic channel to a certain degree.
The conventional geostatistics have problems, and firstly, the variation function with reasonable statistics has strict requirements on the number and distribution of known wells, requires more wells and uniform distribution, and has low prediction precision in areas with quick transverse change or complicated lithological change. Second, for the high frequency stochastic simulation section, conventional geostatistical requires stochastic simulation of frequency ranges up to several hundred hertz (seismic high cut-off frequency less than 100 Hz). The fitting difficulty of the variation function of the frequency information above 100Hz is high, the effect is not ideal, the randomness of the inversion result is high, and the calculation efficiency is low.
Disclosure of Invention
Aiming at least one of the problems, the invention provides a frequency division configuration inversion method suitable for a layered medium, which cuts the impedance modeling or inversion method of four parts of low, medium and high frequency and ultrahigh frequency by using the idea of frequency division construction, each part is calculated by using different modeling and inversion methods, the stability of the inversion result is improved, a high-resolution impedance body is obtained, and the prediction of a thin interbed reservoir is realized. The method has the advantages that random simulation of variation function statistics is replaced by seismic waveform similarity interpolation, stability of a part of results in high-frequency information is enhanced, seismic facies characteristics are better on an inversion result plane, and the continuity is better than that of conventional geostatistics.
In order to achieve the above object, the present invention provides a frequency division configuration inversion method suitable for a layered medium, comprising: s101, filtering a well curve of a true well to obtain low-frequency impedance, counting a low-frequency variation function, and performing interpolation modeling by using a Krigin interpolation method to obtain a low-frequency impedance body; step S102, obtaining an intermediate frequency impedance body by using a colored inversion method, and synthesizing the intermediate frequency impedance body and the low frequency impedance body into an absolute wave impedance body; s103, extracting a well impedance curve of a well side channel of the true well, filtering to obtain high-frequency impedance to form a well set, and extracting the well side channel and a serial number of a pseudo well in a seismic data grid to generate a random path for inverting the pseudo well; s104, performing interpolation modeling by using the well set according to the sequence of the random path to obtain medium-high frequency information of the pseudo well to be interpolated; step S105, performing ultrahigh frequency random simulation on the pseudo well according to the sequence of the random path, establishing a Bayesian discrimination rule by taking the minimum synthetic recorded error residual as a constraint, and averaging a plurality of equal probability results obtained by optimizing by adopting an MCMC algorithm to obtain ultrahigh frequency information; step S106, synthesizing the medium-high frequency information and the ultrahigh frequency information of all the pseudo wells in the order of the random path, and adding the synthesized medium-high frequency information and ultrahigh frequency information into the well set; step S107, carrying out inversion on the residual traces in the seismic data grid by the methods of step S104 and step S105; step S108, matching and combining the impedance information of the absolute wave impedance body and the impedance information of the high-frequency impedance with the corresponding frequency bands of a true well to form full-band impedance, and performing channel-by-channel iteration on the seismic data grids to form a high-resolution inversion body; the low-frequency range is more than 0 and less than 10Hz, the medium frequency is more than or equal to 10 and less than 80Hz, the high frequency comprises the medium-high frequency and the ultrahigh frequency, the medium-high frequency is more than or equal to 80 and less than 200Hz, and the ultrahigh frequency is more than or equal to 200 and less than 500 Hz.
In the foregoing technical solution, preferably, in the step S104, the well side channel of the pseudo well to be interpolated and the well side channel of the real well are used as weights, and the well set is used to perform medium-high frequency interpolation modeling, so as to obtain medium-high frequency information of the pseudo well to be interpolated.
In the foregoing technical solution, preferably, the performing, in step S105, an ultrahigh frequency random simulation on the pseudo well according to the sequence of the random path specifically includes: and sequentially simulating according to the sequence of the random path, guiding the random simulation of the ultrahigh frequency by carrying out probability analysis and statistics on the mean value and the variance of the ultrahigh frequency in the logging data, and simultaneously guiding the simulation of the plane extension of the high frequency by using seismic plane variation analysis.
In the foregoing technical solution, preferably, in step S106, the medium-high frequency information and the ultrahigh frequency information of the pseudo well generated first are synthesized and added to the well set in the order of the random path, and all path points in the random path are iteratively inverted to complete the expansion of the well set.
Compared with the prior art, the invention has the beneficial effects that: by utilizing the idea of frequency division construction, the impedance modeling or inversion method of four parts of low, medium and high frequency is cut apart, each part is calculated by utilizing different modeling and inversion methods, the stability of the inversion result is improved, the impedance body with high resolution is obtained, and the prediction of the thin interbed reservoir is realized. The method has the advantages that random simulation of variation function statistics is replaced by seismic waveform similarity interpolation, stability of a part of results in high-frequency information is enhanced, seismic facies characteristics are better on an inversion result plane, and the continuity is better than that of conventional geostatistics.
Drawings
FIG. 1 is a flow chart of a frequency division configuration inversion method for a layered medium according to an embodiment of the present invention;
FIG. 2 is a schematic block diagram of a flow chart of a frequency division configuration inversion method for a layered medium according to an embodiment of the present invention;
fig. 3 is a schematic diagram of a spectrum construction form of inversion of the frequency division configuration according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
The invention is described in further detail below with reference to the attached drawing figures:
as shown in fig. 1 to 3, the present invention provides a method for inverting a frequency division configuration suitable for a layered medium, including: step S101, performing band-pass filtering on a well curve of a true well to obtain low-frequency impedance (0-10 Hz), counting a low-frequency variation function, and performing interpolation modeling by using a kriging interpolation method to obtain a low-frequency impedance body; step S102, a deterministic inversion part obtains an intermediate frequency impedance body (10-80 Hz) by using a colored inversion method, and synthesizes the intermediate frequency impedance body and the low frequency impedance body into an absolute wave impedance body (the frequency band range is 0 Hz-80 Hz); s103, preparing data, extracting well impedance curves of well side channels of the true wells, filtering to obtain high-frequency impedance (frequency band range is 80-500 Hz) to form a well set, extracting well side channels and serial numbers of the pseudo wells in a seismic data grid (50X 50), and generating a random path for inverting the pseudo wells; step S104, pseudo-well configuration interpolation, namely intermediate-high frequency interpolation modeling, and the realization method comprises the following steps: the method comprises the steps of utilizing the similarity and the distance of a well side channel of a pseudo well to be interpolated and a well side channel of a real well as weights, utilizing a well set to carry out interpolation modeling of medium-high frequency (frequency band range is 80-200 Hz), and obtaining medium-high frequency information of the pseudo well to be interpolated; step S105, randomly simulating the ultrahigh frequency (frequency band range of 200-500 Hz) of the pseudo well: performing ultrahigh frequency random simulation on the pseudo well according to a random path sequence, establishing a Bayesian discrimination rule by taking the minimum of a synthetic recording error residual as a constraint, optimizing by adopting an MCMC (Markov chain Monte Carlo) algorithm, and averaging to obtain a plurality of equal probability results to obtain ultrahigh frequency information; step S106, synthesizing the medium-high frequency information and the ultrahigh frequency information of the pseudo well, adding the synthesized pseudo well into a well set, adding the pseudo well generated firstly into the well set of the real well in the order of random paths to influence the interpolation result of the pseudo well at the next path point, and finishing the expansion of the well set after the inversion of all the path points is finished; step S107, carrying out inversion on the residual channels in the seismic data grid by the methods of step S104 and step S105; s108, matching the impedance information of the 0-80 Hz absolute wave impedance body and the impedance information of the 80-500 Hz high-frequency impedance with the variance and the mean value of the corresponding frequency band of the true well respectively, and combining the impedance information and the mean value into full-frequency-band impedance of 0-500 Hz, and performing channel-by-channel iteration on the seismic data grid to form a high-resolution inversion body; wherein the low-frequency range is more than or equal to 0 and less than 10Hz, the medium-frequency range is more than or equal to 10 and less than 80Hz, the high-frequency range comprises medium-high frequency and ultrahigh frequency, the medium-high frequency range is more than or equal to 80 and less than 200Hz, and the ultrahigh frequency range is more than or equal to 200 and less than 500 Hz.
In the above embodiment, the ideas of frequency division statistical modeling and frequency division inversion are different from the ideas of frequency division statistical variation and synthesis integrated inversion in conventional geostatistics, and are specifically embodied as follows: in the step S101, a low-frequency-band impedance body is established by using kriging interpolation, in the step S102, seismic data are converted into an intermediate-frequency impedance body by using colored inversion, in the steps S103-S107, a high-frequency impedance body is established by using a seismic configuration inversion method, and finally, in the step S108, the three impedance bodies are integrated to obtain a high-resolution impedance body, so that thin interbed reservoir prediction is realized.
In the high-frequency configuration inversion, pseudo-wells in steps S103, S104 and S105 are constructed, the distance between a known well and a point to be interpolated is smaller, the influence degree is larger, the well distribution uniformity also determines the stability of the peripheral point to be interpolated, the construction of the pseudo-wells ensures that the known data are uniformly distributed in the whole network line of the seismic volume, and the randomness and instability in the interpolation process caused by few wells in a high-frequency section are effectively solved. The idea of configuration interpolation in high-frequency configuration inversion is based on the idea that the change of lithology and physical properties can cause the change of frequency, phase and amplitude of seismic data, the form of a waveform is related to the combination of reservoirs, in other words, logging curves are similar if the waveforms are similar, and the extrapolation degree of the logging curves can be judged according to the similarity degree of the seismic waveforms. The degree of similarity of the seismic waveform structure can be used to control extrapolation interpolation of the logging high frequency data. The interpolation method described in step S104, which uses the seismic waveform similarity and the distance between the point to be inverted and the well point as the calculation weight, replaces the interpolation method using the variation of medium-frequency deterministic inversion and the calculation weight of high-frequency logging variation in geostatistical inversion, and avoids the problems of high randomness in calculation of the high-frequency logging variation function and poor plane extrapolation capability.
Specifically, the frequency division configuration inversion divides four impedance modeling or inversion methods of low, medium and high frequency parts into sections according to the idea of frequency division construction, each part is calculated by using different modeling and inversion methods, and a single inversion method has low applicability to data and is easy to cause instability of an inversion result. The low-frequency stability is strong, and the solution is solved by using kriging interpolation; the medium-frequency information can be stably obtained in the seismic data, and the medium-frequency information of the earthquake can be effectively kept by utilizing a colored inversion method; effectively extrapolating the configuration characteristics of the seismic data at medium-high frequency to perform interpolation modeling; the ultrahigh frequency information has strong randomness and cannot be obtained in the earthquake, and is obtained by utilizing an ultrahigh frequency random simulation inversion method in the logging data. In the medium-high frequency configuration modeling method, random simulation of variation function statistics is replaced by seismic waveform similarity interpolation, stability of a part of results in high-frequency information is enhanced, seismic facies characteristics are better on an inversion result plane, and the continuity is better than that of conventional geostatistics.
The method for frequency division configuration inversion suitable for the layered medium, provided by the invention, is used for splitting the impedance modeling or inversion method of the four parts of low, medium and high frequency by utilizing the idea of frequency division construction, each part is calculated by utilizing different modeling and inversion methods, the stability of the inversion result is improved, the high-resolution impedance body is obtained, and the thin interbed reservoir prediction is realized. The method has the advantages that random simulation of variation function statistics is replaced by seismic waveform similarity interpolation, stability of a part of results in high-frequency information is enhanced, seismic facies characteristics are better on an inversion result plane, and the continuity is better than that of conventional geostatistics.
The above is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes will occur to those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (4)
1. A method for inversion of a frequency division configuration for a layered medium, comprising:
s101, filtering a well curve of a true well to obtain low-frequency impedance, counting a low-frequency variation function, and performing interpolation modeling by using a Krigin interpolation method to obtain a low-frequency impedance body;
step S102, obtaining an intermediate frequency impedance body by using a colored inversion method, and synthesizing the intermediate frequency impedance body and the low frequency impedance body into an absolute wave impedance body;
s103, extracting a well impedance curve of a well side channel of the true well, filtering to obtain high-frequency impedance to form a well set, and extracting the well side channel and a serial number of a pseudo well in a seismic data grid to generate a random path for inverting the pseudo well;
s104, performing interpolation modeling by using the well set according to the sequence of the random path to obtain medium-high frequency information of the pseudo well to be interpolated;
step S105, performing ultrahigh frequency random simulation on the pseudo well according to the sequence of the random path, establishing a Bayesian discrimination rule by taking the minimum synthetic recorded error residual as a constraint, and averaging a plurality of equal probability results obtained by optimizing by adopting an MCMC algorithm to obtain ultrahigh frequency information;
step S106, synthesizing the medium-high frequency information and the ultrahigh frequency information of all the pseudo wells in the order of the random path, and adding the synthesized medium-high frequency information and ultrahigh frequency information into the well set;
step S107, carrying out inversion on the residual traces in the seismic data grid by the methods of step S104 and step S105;
step S108, matching and combining the impedance information of the absolute wave impedance body and the impedance information of the high-frequency impedance with the corresponding frequency bands of a true well to form full-band impedance, and performing channel-by-channel iteration on the seismic data grids to form a high-resolution inversion body;
the low-frequency range is more than 0 and less than 10Hz, the medium frequency is more than or equal to 10 and less than 80Hz, the high frequency comprises the medium-high frequency and the ultrahigh frequency, the medium-high frequency is more than or equal to 80 and less than 200Hz, and the ultrahigh frequency is more than or equal to 200 and less than 500 Hz.
2. The inversion method of frequency division configuration suitable for layered media according to claim 1, wherein in step S104, the well side channel of the pseudo well to be interpolated and the well side channel of the real well are used as weights, and the well set is used to perform interpolation modeling of medium and high frequencies, so as to obtain medium and high frequency information of the pseudo well to be interpolated.
3. The inversion method of frequency division configuration suitable for layered media according to claim 1, wherein the step S105 of performing the random simulation of the uhf on the pseudo well according to the sequence of the random path specifically comprises:
and sequentially simulating according to the sequence of the random path, guiding the random simulation of the ultrahigh frequency by counting the mean value and the variance of the ultrahigh frequency in the logging data through probability analysis, and simultaneously guiding the simulation of the plane extension of the high-frequency component through plane variation analysis of the seismic data.
4. The inversion method of frequency division configuration suitable for layered media according to claim 1, wherein in step S106, the medium-high frequency information and the ultrahigh frequency information of the pseudo-well generated in advance are synthesized and added to the well set in the order of the random path, and all path points in the random path are iteratively inverted to complete the expansion of the well set.
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 CN108919344A (en) | 2018-11-30 |
| CN108919344B true 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) |
Families Citing this family (3)
| 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 |
| CN113917529B (en) * | 2020-07-07 | 2024-06-18 | 中国石油化工股份有限公司 | Wave impedance inversion method and device, storage medium and electronic equipment |
| CN112784379B (en) * | 2021-03-09 | 2023-05-23 | 重庆邮电大学 | Optimal design method of impeller rotor system |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104570066B (en) * | 2013-10-10 | 2017-02-08 | 中国石油天然气股份有限公司 | 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 |
| CN105425289B (en) * | 2015-10-29 | 2018-08-14 | 中国石油天然气集团公司 | The method and apparatus for determining low frequency wave impedance |
| CN105891888B (en) * | 2016-03-28 | 2017-03-08 | 吉林大学 | Multiple domain divides multiple dimensioned full waveform inversion method parallel |
| CN106707339B (en) * | 2016-11-22 | 2018-09-25 | 中国石油天然气集团有限公司 | A kind of phased inversion method of Depth Domain |
-
2018
- 2018-03-30 CN CN201810293307.3A patent/CN108919344B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108919344A (en) | 2018-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108919344B (en) | Frequency division configuration inversion method suitable for layered medium | |
| Miyake et al. | Source characterization for broadband ground-motion simulation: Kinematic heterogeneous source model and strong motion generation area | |
| CN108645994B (en) | Geological random inversion method and device based on multipoint geostatistics | |
| CN101349764B (en) | Earthquake gyration analysis method | |
| CN108267784A (en) | A kind of seismic signal random noise compression process method | |
| CN103439740B (en) | Method and device for predicting relative impedance based on dipole seismic wavelet multiple integral | |
| CN103728660A (en) | Multi-channel matching tracking method based on seismic data | |
| CN109738950B (en) | The noisy-type data primary wave inversion method of domain inverting is focused based on sparse 3 D | |
| KR101931092B1 (en) | Hybrid method for full waveform inversion using simultaneous and sequential source method | |
| CN104122581A (en) | Poststack acoustic wave impedance inversion method | |
| CN103364826A (en) | An earthquake blind source deconvolution method based on independent component analysis | |
| Yin et al. | Improving horizontal resolution of high-frequency surface-wave methods using travel-time tomography | |
| Emolo et al. | Kinematic source parameters for the 1989 Loma Prieta earthquake from the nonlinear inversion of accelerograms | |
| US11143774B2 (en) | Method and system for separating blended seismic data | |
| CN106249282A (en) | A kind of frequency domain seismic channel set creation method being applicable to AVAF inverting | |
| CN105785443B (en) | The method of high-precision trace integral calculating Relative Wave Impedance | |
| CN102830432A (en) | Method for identifying weak reflection reservoir under cover of coal series strong earthquake reflection characteristics | |
| Mueller et al. | Optimizing near-orthogonal air-gun firing sequences for marine simultaneous source separation | |
| CN111914609B (en) | Well-seismic combined prestack geostatistical elastic parameter inversion method and device | |
| Dittrich et al. | An efficient method for avoiding shadow fading maps in system level simulations | |
| Mironov et al. | Strong motion record processing of the Baikal rift zone | |
| Fu et al. | USING THE WAVELET TRANSFORM FOR SEISMIC WAVE IMPEDANCE INVERSION | |
| CN114966826A (en) | High-low frequency fusion near-surface modeling method based on field near-surface survey data | |
| Vyzhva et al. | The advanced algorithm of statistical simulation of seismic noise in the multidimensional area for determination the frequency characteristics of geological environment | |
| Wang et al. | Stochastic inversion based on spectral simulation of non-stationary statistical characteristic parameters |
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 |