CN108732623B - TTI medium tomography inversion imaging method and computer readable storage medium - Google Patents

TTI medium tomography inversion imaging method and computer readable storage medium Download PDF

Info

Publication number
CN108732623B
CN108732623B CN201710250240.0A CN201710250240A CN108732623B CN 108732623 B CN108732623 B CN 108732623B CN 201710250240 A CN201710250240 A CN 201710250240A CN 108732623 B CN108732623 B CN 108732623B
Authority
CN
China
Prior art keywords
parameters
anisotropy
inversion
delta
epsilon
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
Application number
CN201710250240.0A
Other languages
Chinese (zh)
Other versions
CN108732623A (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.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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 China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201710250240.0A priority Critical patent/CN108732623B/en
Publication of CN108732623A publication Critical patent/CN108732623A/en
Application granted granted Critical
Publication of CN108732623B publication Critical patent/CN108732623B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/67—Wave propagation modeling

Abstract

The invention discloses a chromatographic inversion imaging method of a TTI medium and a computer readable storage medium, comprising the following steps: inputting data; based on TTI mediaAnisotropy parameter (v)p0,ε,δ,θ0) Inverting the anisotropic parameters epsilon and delta to obtain accurate results of the parameters epsilon and delta; based on the results of the exact parameters ε and δ, for the anisotropy parameters ε, δ and vp0Carrying out inversion to obtain accurate parameters epsilon, delta and vp0The result of (1); based on the exact parameters epsilon, delta and vp0As a result of (3), for the anisotropy parameters ε, δ, vp0And theta0Performing inversion to obtain the anisotropic parameter (v)p0,ε,δ,θ0) The inversion result of (2). The advantages are that: by inverting different sensitivity degrees of different parameters in the TTI anisotropic medium during traveling, the accuracy of the chromatographic modeling of the TTI medium is greatly improved, and the accuracy of offset imaging can be improved.

Description

TTI medium tomography inversion imaging method and computer readable storage medium
Technical Field
The invention relates to the field of seismic exploration, in particular to a Transmission Time Interval (TTI) medium tomography inversion imaging method and a computer readable storage medium.
Background
The traditional seismic migration imaging technology is mostly based on the assumption that the earth medium is isotropic, but a large number of observations and experiments show that the anisotropy of the medium in the earth is ubiquitous, and if the migration imaging method based on the isotropic assumption is used for processing seismic data of an anisotropic exploration area, the precision of velocity analysis is seriously influenced, and then the accuracy of seismic migration imaging is reduced.
Anisotropy can be classified into VTI perpendicular transverse anisotropy, TTI tilt transverse anisotropy, and HTI azimuthal anisotropy. Among them, VTI media are currently being studied more, while TTI and HTI are being studied less. In the field of tomographic inversion velocity modeling in seismic exploration, due to the problem that the magnitude of each parameter of a TTI medium is not uniform, the algorithm is always dedicated to research on the influence of different inversion strategies on TTI multi-parameter inversion, and therefore the research on the tomographic velocity modeling of the TTI medium is delayed to a certain extent.
Therefore, it is necessary to develop a TTI media tomography inversion imaging method based on the sensitivity of different TTI anisotropy parameters to travel.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention provides a TTI medium tomography inversion imaging method and a computer readable storage medium, which iteratively updates a velocity model by researching the sensitivity of different TTI anisotropic parameters to travel and performing tomography inversion by adopting a parameterization means and a strategy, provides an accurate TTI anisotropic velocity model and lays a good foundation for offset imaging.
According to one aspect of the present invention, a TTI media tomography method is provided, including:
inputting data;
anisotropy parameter (v) based on TTI Mediump0,ε,δ,θ0) Inverting the anisotropic parameters epsilon and delta to obtain accurate results of the parameters epsilon and delta;
based on the results of the exact parameters ε and δ, for the anisotropy parameters ε, δ and vp0Carrying out inversion to obtain accurate parameters epsilon, delta and vp0The result of (1);
based on the precise parameters ε, δ, and vp0For the anisotropy parameters ε, δ, vp0And theta0Performing inversion to obtain the anisotropic parameter (v)p0,ε,δ,θ0) The inversion result of (2);
wherein v isp0The qp wave phase velocity along the direction of the symmetry axis; ε and δ are Thomsen parameters; theta0Is the angle between the normal direction of the wave front and the symmetry axis.
Preferably, the anisotropy parameter v is inverted for the anisotropy parameters ε and δp0And theta0Set to 0, proceedAnd obtaining the inversion result of the anisotropy parameters epsilon and delta through multiple iterations.
Preferably, for the anisotropy parameters ε, δ and vp0The anisotropy parameter theta is obtained during inversion0Setting the parameters as 0, and carrying out multiple iterations to obtain the anisotropy parameters epsilon, delta and vp0The inversion result of (2).
Preferably, for the anisotropy parameters ε, δ, vp0And theta0Inverting, and performing multiple iterations to obtain the anisotropic parameters epsilon, delta and vp0And theta0The inversion result of (2).
According to another aspect of the invention, a computer-readable storage medium is proposed, on which a computer program is stored, wherein the program realizes the following steps when executed by a processor:
anisotropy parameter (v) based on TTI Mediump0,ε,δ,θ0) Inverting the anisotropic parameters epsilon and delta to obtain accurate results of the parameters epsilon and delta;
based on the results of the exact parameters ε and δ, for the anisotropy parameters ε, δ and vp0Carrying out inversion to obtain accurate parameters epsilon, delta and vp0The result of (1);
based on the precise parameters ε, δ, and vp0For the anisotropy parameters ε, δ, vp0And theta0Performing inversion to obtain the anisotropic parameter (v)p0,ε,δ,θ0) The inversion result of (2);
wherein v isp0The qp wave phase velocity along the direction of the symmetry axis; ε and δ are Thomsen parameters; theta0Is the angle between the normal direction of the wave front and the symmetry axis.
Preferably, the anisotropy parameter v is inverted for the anisotropy parameters ε and δp0And theta0And setting the parameters as 0, and carrying out multiple iterations to obtain the inversion results of the anisotropy parameters epsilon and delta.
Preferably, for the anisotropy parameters ε, δ and vp0The anisotropy parameter theta is obtained during inversion0Setting the value to be 0, and carrying out multiple iterations to obtain the anisotropic parametersε, δ and vp0The inversion result of (2).
Preferably, for the anisotropy parameters ε, δ, vp0And theta0Inverting, and performing multiple iterations to obtain the anisotropic parameters epsilon, delta and vp0And theta0The inversion result of (2).
According to the TTI medium tomography inversion imaging method and the computer readable storage medium, the advantages are that: by inverting different sensitivity degrees of different parameters in the TTI anisotropic medium during traveling, the accuracy of the chromatographic modeling of the TTI medium is greatly improved, and the accuracy of offset imaging can be improved.
The method and computer-readable storage medium of the present invention have other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the present invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts.
Figure 1 shows a graph of sensitivity of TTI medium anisotropy parameters to travel.
FIG. 2 shows a flow chart of the steps of a TTI media tomography inversion imaging method in accordance with the present invention.
FIG. 3 shows the parameter v of the anisotropy parameter of the TTI medium in the actual modelp0Schematic representation of imaging.
Fig. 4a and 4b show schematic diagrams of the results of an epsilon and delta two-parameter inversion, respectively, according to an exemplary embodiment of the present invention.
FIGS. 5a, 5b and 5c respectively show v according to an exemplary embodiment of the present inventionp0Schematic diagram of results of three-parameter inversion of epsilon and delta.
FIG. 6a and FIG. 6bb. FIGS. 6c and 6d show v, respectively, according to an exemplary embodiment of the present inventionp0ε, δ and θ0And (5) a four-parameter inversion result schematic diagram.
Detailed Description
The invention will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
Due to TTI anisotropy parameter (v)p0,ε,δ,θ0) The scale of the inversion update volume varies for different effects on travel, and therefore sensitivity analysis of these parameters is necessary.
The relationship between the travel time and the anisotropy parameter of the TTI medium is as follows:
Figure BDA0001271882090000041
wherein v isp0The qp wave phase velocity along the direction of the symmetry axis; ε and δ are Thomsen parameters; theta0Is the angle between the normal direction of the wave front and the symmetry axis.
Shifting the right derivative in equation (1) to the left, the warping yields:
Figure BDA0001271882090000051
figure 1 shows a graph of sensitivity of TTI medium anisotropy parameters to travel.
The sensitivity of the anisotropy parameter of the TTI medium to the travel time can be plotted by the equation (2), as shown in FIG. 1, and it can be seen that the sensitivity of the anisotropy parameter of the TTI medium to the travel time is different, where ε is the largest, δ is the second, vp0And theta0The method is least sensitive, so that the method determines to firstly carry out inversion strategy setting on two parameters of epsilon and delta which are most sensitive.
The TTI medium tomography inversion imaging method can comprise the following steps: inputting data; anisotropy parameter (v) based on TTI Mediump0,ε,δ,θ0) Inverting the anisotropic parameters epsilon and delta to obtain accurate results of the parameters epsilon and delta; based on the results of the exact parameters ε and δ, for the anisotropy parameters ε, δ and vp0Carrying out inversion to obtain accurate parameters epsilon, delta and vp0The result of (1); based on the exact parameters epsilon, delta and vp0As a result of (3), for the anisotropy parameters ε, δ, vp0And theta0Performing inversion to obtain the anisotropic parameter (v)p0,ε,δ,θ0) The inversion result of (2).
Wherein v isp0The qp wave phase velocity along the direction of the symmetry axis; ε and δ are Thomsen parameters; theta0Is the angle between the normal direction of the wave front and the symmetry axis.
According to different sensitivity degrees of anisotropic parameters of the TTI medium to traveling, a targeted inversion strategy is provided, and the accuracy of chromatographic modeling of the TTI medium can be greatly improved.
In the geophysical field, tomography inversion is that under the premise that known medium anisotropy parameters are possibly inaccurate, certain residual exists between the calculation travel time and the observation travel time (namely, actual travel time) of a ray passing through a medium, and the residual is caused by inaccuracy of the medium anisotropy parameters, so that the medium parameters can be inverted according to the theory that the integral of the medium anisotropy parameters along the ray path is the residual during travel. Therefore, according to the relationship between the anisotropy parameters and the travel time residual errors deduced above, the inversion of the anisotropy parameters can be realized through chromatographic inversion.
As a preferred scheme, when the anisotropy parameters epsilon and delta are inverted, the anisotropy parameter vp0And theta0And setting the value to be 0, and carrying out multiple iterations to obtain the inversion results of the anisotropy parameters epsilon and delta.
For epsilon and delta two-parameter inversion, according to a TTI medium chromatographic inversion equation shown in the formula (1), the left side of the equation is a residual error between travel time calculation and travel time observation under actual model parameters, the right side of the equation can be understood as an integral of anisotropic parameters of different TTI media along a ray, and according to an inversion means such as a conjugate gradient method and the like, equation solution can be realized to obtain the anisotropic parameter value of the inverted TTI medium.
Preferably, the anisotropy parameters ε, δ and vp0Anisotropy parameter θ in inversion0Setting the parameters as 0, and carrying out multiple iterations to obtain anisotropy parameters epsilon, delta and vp0The inversion result of (2).
Preferably, the anisotropy parameters ε, δ and vp0And theta0Performing inversion, and performing multiple iterations to obtain anisotropy parameters epsilon, delta and vp0And theta0The inversion result of (2).
Through testing of the model, the inversion strategy provided by the invention is verified to be capable of realizing accurate inversion of the anisotropic parameters of the TTI medium, and the result can be used as input data of subsequent processing steps, so that the offset imaging precision is improved.
Examples
To facilitate an understanding of the aspects of the embodiments of the present invention and their effects, a specific example is given below. It will be understood by those skilled in the art that this example is merely for the purpose of facilitating an understanding of the present invention and that any specific details thereof are not intended to limit the invention in any way.
FIG. 2 shows a flow chart of the steps of a TTI media tomography inversion imaging method in accordance with the present invention.
The invention relates to a TTI medium tomography inversion imaging method, which comprises the following steps:
inputting data;
anisotropy parameter (v) based on TTI Mediump0,ε,δ,θ0) Inverting the anisotropic parameters epsilon and delta to obtain accurate results of the parameters epsilon and delta;
based on the results of the exact parameters ε and δ, for the anisotropy parameters ε, δ and vp0Carrying out inversion to obtain accurate parameters epsilon, delta and vp0The result of (1);
based on the exact parameters epsilon, delta and vp0As a result of (3), for the anisotropy parameters ε, δ, vp0And theta0Performing inversion to obtain the anisotropic parameter (v)p0,ε,δ,θ0) The inversion result of (2);
wherein v isp0The qp wave phase velocity along the direction of the symmetry axis; ε and δ are Thomsen parameters; theta0Is the angle between the normal direction of the wave front and the symmetry axis.
FIG. 3 shows the parameter v of the anisotropy parameter of the TTI medium in the actual modelp0Schematic representation of imaging.
Since the imaging profiles of the anisotropy parameters of the TTI media are uniform, only the parameter v is used thereinp0For example, as shown in FIG. 3, the imaging of the anisotropy parameters of the TTI media is nearly circular.
In this embodiment, an abnormal body model is adopted, the model grid is 201 × 201, the size of the grid is 5m × 5m, and the measured anisotropy parameter of the actual model is vp0＝3000m·s-1,ε＝0.1,δ＝0,θ 00 DEG, and the anomalous anisotropy parameter is vp0＝3300m·s-1,ε＝0.2,δ＝0.1,θ030 ° is set. The observation mode of surface excitation and underground reception is adopted, 41 cannons are used, 41 cannon shots are used, and shot point receiving points are uniformly distributed.
Fig. 4a and 4b show schematic diagrams of the results of an epsilon and delta two-parameter inversion, respectively, according to an exemplary embodiment of the present invention.
The two parameters of epsilon and delta are inverted at the same time, the initial values of the other parameters are set to be 0, after 10 times of iteration, the inversion result is shown in figure 4a and figure 4b, the inverted epsilon and delta are close to the actual model, such as the parameter imaging shown in figure 3, and the inversion precision is high.
FIGS. 5a, 5b and 5c respectively show v according to an exemplary embodiment of the present inventionp0Schematic diagram of results of three-parameter inversion of epsilon and delta.
Taking the results of the above inversions as input data, and performing epsilon, delta and vp0Three parameters are inverted simultaneously, the initial values of the other parameters are set to be 0, after 10 iterations, the inversion results are shown in figure 5a, figure 5b and figure 5c, and the inverted epsilon, delta and vp0The inversion accuracy is higher as the parameter imaging is close to the actual model, such as shown in FIG. 3.
Fig. 6a, 6b, 6c and 6d show respectively v according to an exemplary embodiment of the present inventionp0ε, δ and θ0And (5) a four-parameter inversion result schematic diagram.
The epsilon, delta and v obtained by the inversionp0The results of (3) are input as data, and epsilon, delta and v are performedp0And theta0The 4 anisotropic parameters are inverted simultaneously, and after 10 iterations, the inversion results are shown in FIG. 6a, FIG. 6b, FIG. 6c and FIG. 6d, wherein ε, δ and vp0The inversion result of the method is close to the actual model, such as the parameter imaging shown in FIG. 3, and the inversion precision is high. Only theta0The inversion result of (2) has a gap from the parameters of the actual model, which is related to the fact that the parameters themselves are too weak to be sensitive to travel.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the disclosed embodiments. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the illustrated embodiments. The terminology used herein is chosen in order to best explain the principles of the embodiments, the practical application, or improvements made to the technology in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein.

Claims (8)

1. A TTI media tomography inversion imaging method, comprising:
inputting data;
anisotropy parameter (v) based on TTI Mediump0,ε,δ,θ0) Inverting the anisotropic parameters epsilon and delta to obtain accurate results of the parameters epsilon and delta;
based on the results of the exact parameters ε and δ, for the anisotropy parameters ε, δ and vp0Carrying out inversion to obtain accurate parameters epsilon, delta and vp0The result of (1);
based on the precise parameters ε, δ, and vp0For the anisotropy parameters ε, δ, vp0And theta0Performing inversion to obtain eachParameter of anisotropy (v)p0,ε,δ,θ0) The inversion result of (2);
wherein v isp0The qp wave phase velocity along the direction of the symmetry axis; ε and δ are Thomsen parameters; theta0Is the angle between the normal direction of the wave front and the symmetry axis.
2. The TTI media tomography method of claim 1, wherein the anisotropy parameter v is inverted when inverting the anisotropy parameters epsilon and deltap0And theta0And setting the parameters as 0, and carrying out multiple iterations to obtain the inversion results of the anisotropy parameters epsilon and delta.
3. The TTI media tomography method of claim 1, wherein the anisotropy parameters epsilon, delta, and v are measuredp0The anisotropy parameter theta is obtained during inversion0Setting the parameters as 0, and carrying out multiple iterations to obtain the anisotropy parameters epsilon, delta and vp0The inversion result of (2).
4. The TTI media tomography method of claim 1, wherein the anisotropy parameters epsilon, delta, v are scaledp0And theta0Inverting, and performing multiple iterations to obtain the anisotropic parameters epsilon, delta and vp0And theta0The inversion result of (2).
5. A computer-readable storage medium, on which a computer program is stored, wherein the program, when executed by a processor, performs the steps of:
anisotropy parameter (v) based on TTI Mediump0,ε,δ,θ0) Inverting the anisotropic parameters epsilon and delta to obtain accurate results of the parameters epsilon and delta;
based on the results of the exact parameters ε and δ, for the anisotropy parameters ε, δ and vp0Carrying out inversion to obtain accurate parameters epsilon, delta and vp0The result of (1);
based on the precise parameters ε, δ, and vp0As a result, with respect to the anisotropy parameter ε,δ、vp0And theta0Performing inversion to obtain the anisotropic parameter (v)p0,ε,δ,θ0) The inversion result of (2);
wherein v isp0The qp wave phase velocity along the direction of the symmetry axis; ε and δ are Thomsen parameters; theta0Is the angle between the normal direction of the wave front and the symmetry axis.
6. The computer-readable storage medium of claim 5, wherein the anisotropy parameter, ν when inverting the anisotropy parameters, ε and δp0And theta0And setting the parameters as 0, and carrying out multiple iterations to obtain the inversion results of the anisotropy parameters epsilon and delta.
7. The computer-readable storage medium of claim 5, wherein the anisotropy parameters ε, δ, and v are measuredp0The anisotropy parameter theta is obtained during inversion0Setting the parameters as 0, and carrying out multiple iterations to obtain the anisotropy parameters epsilon, delta and vp0The inversion result of (2).
8. The computer-readable storage medium of claim 5, wherein the anisotropy parameters ε, δ, v are measuredp0And theta0Inverting, and performing multiple iterations to obtain the anisotropic parameters epsilon, delta and vp0And theta0The inversion result of (2).
CN201710250240.0A 2017-04-17 2017-04-17 TTI medium tomography inversion imaging method and computer readable storage medium Active CN108732623B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710250240.0A CN108732623B (en) 2017-04-17 2017-04-17 TTI medium tomography inversion imaging method and computer readable storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710250240.0A CN108732623B (en) 2017-04-17 2017-04-17 TTI medium tomography inversion imaging method and computer readable storage medium

Publications (2)

Publication Number Publication Date
CN108732623A CN108732623A (en) 2018-11-02
CN108732623B true CN108732623B (en) 2020-05-19

Family

ID=63924955

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710250240.0A Active CN108732623B (en) 2017-04-17 2017-04-17 TTI medium tomography inversion imaging method and computer readable storage medium

Country Status (1)

Country Link
CN (1) CN108732623B (en)

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100135115A1 (en) * 2008-12-03 2010-06-03 Chevron U.S.A. Inc. Multiple anisotropic parameter inversion for a tti earth model
CN104407378B (en) * 2014-11-25 2017-05-10 中国石油天然气股份有限公司 Method and device for inversing anisotropy parameters

Also Published As

Publication number Publication date
CN108732623A (en) 2018-11-02

Similar Documents

Publication Publication Date Title
Chu et al. Pure acoustic wave propagation in transversely isotropic media by the pseudospectral method
CN105005076B (en) Seismic wave full waveform inversion method based on least square gradient update speed model
US8571842B2 (en) Method of determining parameter from sparse measurement data
US8594939B2 (en) Method for calibrating seismic imaging velocities
CN104775811B (en) A kind of extraction of formation anisotropy information and bearing calibration and system
US20130258810A1 (en) Method and System for Tomographic Inversion
CN106483559B (en) A kind of construction method of subsurface velocity model
CN107817516B (en) Near-surface modeling method and system based on first-motion wave information
CN107229071B (en) A kind of subsurface structure inversion imaging method
Zhou et al. Multiparameter joint tomography for TTI model building
Zhang et al. Moveout-based wave-equation migration velocity analysis
CN109946741A (en) Pure qP wave least square reverse-time migration imaging method in a kind of TTI medium
Zhou et al. Nonlinear traveltime inversion scheme for crosshole seismic tomography in tilted transversely isotropic media
Pšenčík et al. Reflection moveout approximations for P-waves in a moderately anisotropic homogeneous tilted transverse isotropy layer
US10302801B2 (en) Temperature correction of magnetic resonance data
CN110989017A (en) Ground nuclear magnetic resonance inversion method containing variable frequency offset
CN107229066B (en) VSP data full waveform inversion modeling method based on surface seismic structure constraint
CN112773396B (en) Medical imaging method based on full waveform inversion, computer equipment and storage medium
CN108732623B (en) TTI medium tomography inversion imaging method and computer readable storage medium
Gong et al. Combined migration velocity model-building and its application in tunnel seismic prediction
CN104316961A (en) Method for obtaining geological parameters of weathered layer
Mothi et al. Detecting and estimating anisotropy errors using full waveform inversion and ray-based tomography
CN115236755B (en) Tensor eigenvalue-based aeromagnetic anomaly boundary detection method and device
CN112731520A (en) Full waveform inversion method and system based on structure tensor diffusion filtering
Zhang et al. Wave-equation migration velocity analysis using partial stack-power maximization

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