CN108957554B - Seismic inversion method in geophysical exploration - Google Patents
Seismic inversion method in geophysical exploration Download PDFInfo
- Publication number
- CN108957554B CN108957554B CN201810903194.4A CN201810903194A CN108957554B CN 108957554 B CN108957554 B CN 108957554B CN 201810903194 A CN201810903194 A CN 201810903194A CN 108957554 B CN108957554 B CN 108957554B
- Authority
- CN
- China
- Prior art keywords
- seismic
- data
- layer
- force
- fractal
- 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
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Abstract
The invention discloses a seismic inversion method in geophysical exploration, which comprises the steps of obtaining logging data in a seismic work area, and collecting and editing force and moment data of corresponding points of each seismic layer; layer-by-layer interpretation is carried out on each layer of seismic data, and interpretation result data are gridded to obtain a geological framework model based on depth; determining a logging curve according to the determined geological framework model, and performing fractal random interpolation on a fractal parameter body of the logging curve to obtain an initial logging curve; and determining an actual logging curve according to the seismic data, correcting the initial logging curve according to the comparison result, and determining a seismic inversion model according to the corrected curve. The method comprises the steps of acquiring logging data in a seismic work area, collecting and editing force and moment data of corresponding points of each seismic layer, introducing various parameters of force and moment, establishing an exploration model by combining a specific soil environment, and determining a final inversion model by correcting a final value.
Description
Technical Field
The invention relates to the technical field of seismic inversion, in particular to a seismic inversion method in geophysical exploration.
Background
In the oil and gas industry, seismic inversion is used as a core technology for reservoir prediction, the high-resolution advantage in the longitudinal direction of logging and the dense sampling advantage in the transverse direction of earthquake can be closely combined, characteristic information such as lithology, physical property and the like in the transverse direction of a stratum can be estimated, and the method is widely applied to various stages of oil exploration and development at present. However, the conventional time domain seismic inversion method has poor identification capability for thin reservoirs, and a great deal of valuable high-level information is lost in the process of converting the well logging curve from a depth domain to a time domain, so that the method has great limitations.
In recent years, the most common practice for seismic inversion reservoir prediction has been depth domain seismic data inversion based on convolution models. The method has the first problem of wavelet acquisition, namely, the time domain wavelet is converted into the depth domain wavelet through time-depth conversion, and then the convolution effect of the wavelet is eliminated according to the convolution and inverse wavelet of the seismic channel, so that an approximate stratum reflection coefficient series is obtained for inverting a high-resolution stratum type section.
The logging constraint seismic inversion technology plays an important role in lithology identification and reservoir prediction, the seismic data with intensive spatial sampling and the logging data with high vertical resolution are combined, the interface reflection type seismic data are inverted to reflect the wave impedance data volume of the stratum, the identification precision of the reservoir is improved, the thickness change condition of the reservoir can be more finely described, particularly, in the aspect of thin reservoir identification, the resolution is improved through logging constraint, and thin layers which cannot be identified by conventional earthquakes can be identified.
The existing inversion method is mainly based on an amplitude inversion method, namely the change of the seismic amplitude strength is considered to reflect the spatial change of a reservoir stratum, most of the inversion methods are based on a model, the building process of the model is mainly based on the spatial interpolation and extrapolation of logging attributes under the constraint of a seismic interpretation horizon, and the change condition of a geologic body in the transverse direction is not considered in the whole modeling process, so that although the model-based inversion method can improve the identification capability of a thin layer, errors exist in the identification of a spatial pinch-off point of the reservoir stratum.
Therefore, the seismic inversion method in the prior art carries out inversion according to the model and waveform identification, and does not consider the actual physical factors and stress factors.
Disclosure of Invention
The invention aims to provide a seismic inversion method in geophysical exploration, which aims to solve the technical problem that stress factors are not considered in a model.
In order to achieve the purpose, the invention provides the following technical scheme: the method comprises the following steps:
step a, acquiring logging data in a seismic work area, and acquiring and editing force and moment data of corresponding points of each seismic layer;
b, performing layer-by-layer interpretation on each layer of seismic data, and gridding interpretation result data to obtain a geological frame model based on depth;
step c, determining a logging curve according to the determined geological framework model, and performing fractal random interpolation on a fractal parameter body of the logging curve to obtain an initial logging curve;
d, determining an actual logging curve according to the seismic data, correcting the initial logging curve according to a comparison result, and determining a seismic inversion model according to the corrected curve;
in the step a, a plurality of sampling points of each seismic layer are established according to the following model: one of the sampling points is modeled as:
wherein the content of the first and second substances,the zero offset values of the detection force and the moment of the seismic layer sampling point,is the gravity and gravity moment of the sampling point,andexternal force and external moment applied to the sampling point are included, wherein the moment is the moment of the sampling point relative to the central line of the seismic layer cavity;
setting the vector from the origin of the central coordinate system of the seismic layer cavity to the centroid of the sampling point asThen
The detailed expression of equation (1) in the survey coordinate system based on the center of the seismic layer void is:
wherein the content of the first and second substances,is a transformation matrix from the world coordinate system to the central coordinate system of the seismic interval cavity, determined by the initial position of the bytes, xc,yc,zc;
When the sampling point is positioned on the upper surface layer of the seismic layer cavity:
detecting the existence relation between force and external force:
F1representing the detected value of force, m, at the position of the upper layerloadRepresenting the real-time gravity of the geology itself, e.g. the weight of the sand, gravel, sand itself, fextExternal force, F0x,F0y,F0zIndicating the resultant force, T, in each direction of the position0x,T0y,T0zIndicating the resultant moment, x, in each direction of the positionc,yc,zcRepresenting the coordinates at the location, enabling the location of the sampling point to be uniquely determined;
when the sampling point is positioned on the lower surface of the seismic layer cavity:
detecting the existence relation between force and external force:
when the sampling point is positioned in the middle layer of the seismic layer cavity:
detecting the existence relation between force and external force:
through collecting the force detection measurement data in the motion process, the force detection zero offset and the self load weight parameter to be identified can be synchronously obtained:
finally determine x by the above equation (9)c,yc,zcAnd the coordinates of the positions can uniquely determine the stress and the moment of the sampling points and the positions of the corresponding sampling points.
Further, in the step d, the correction is performed according to the following formula (10) in accordance with the correction of each selected point;
im=ρ×i0 (10)
wherein imRepresents the instantaneous value of the sample point after correction, rho represents the correction factor, i0Representing the instantaneous value of the sample point; the correction coefficient p is calculated by the following formula (11),
where ρ represents a correction coefficient, i01And i02The instantaneous values of two points at the position of deviation when the deviation occurs are shown, N represents the number of sampling times, and k represents the sampling sequence.
Furthermore, in the step b, a plurality of sets of reflection positions are interpreted and tracked for the seismic data, the quality of the interpreted positions is checked and modified, and the positions are subjected to full three-dimensional tracking interpretation, so that the position interpretation result is fine, and each earthquake has a fine position interpretation result.
Further, in the step c, calculating a fractal parameter of a logging curve model of a single well, and converting the fractal parameter of each seismic channel into a fractal parameter body of the logging curve; respectively calculating fractal parameters of a logging curve depth domain model of a single well, calibrating the seismic fractal parameters of each seismic channel into a fractal parameter body of the logging curve by methods of normalization, proportion weighting and the like, thus obtaining the fractal parameters of each seismic channel of the logging curve, and then calculating the random component funine of fractal random interpolation of each seismic channel according to the following formula:
k is a calibration coefficient, H is a fractal dimension parameter reflecting details and roughness of the hierarchical domain post-stack seismic data, sigma is the variance of the hierarchical domain post-stack seismic data normal distribution, G is a Gauss random variable obeying the standard normal distribution, Delta X is the distance between a known point and an interpolation point, and a is a one-dimensional coordinate of the interpolation point.
Compared with the prior art, the invention has the beneficial effects that: a seismic inversion method in geophysical exploration comprises the steps of acquiring logging data in a seismic work area, collecting and editing force and moment data of corresponding points of each seismic layer, introducing various parameters of force and moment, establishing an exploration model by combining a specific soil environment, and determining a final inversion model by correcting a final value.
Further, the calibration method of the embodiment of the invention comprises the following steps: the identification time is short, only force detection data need to be collected, and a simulated exploration-based earthquake stress model is determined by adopting detection force, moment, stress position and combining external force and gravity information of different soil textures.
The invention ensures the stability of the final value by adopting a correction mode and comparing and correcting the empirical value and the real-time value.
Drawings
FIG. 1 is a logic flow diagram of a method of seismic inversion in geophysical exploration according to the present invention.
Detailed Description
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 only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, an embodiment of the invention provides a seismic inversion method in geophysical exploration, which specifically includes:
step a, acquiring logging data in a seismic work area, and acquiring and editing force and moment data of corresponding points of each seismic layer;
b, performing layer-by-layer interpretation on each layer of seismic data, and gridding interpretation result data to obtain a geological frame model based on depth;
step c, determining a logging curve according to the determined geological framework model, and performing fractal random interpolation on a fractal parameter body of the logging curve to obtain an initial logging curve;
and d, determining an actual logging curve according to the seismic data, correcting the initial logging curve according to the comparison result, and determining a seismic inversion model according to the corrected curve.
In this embodiment, in step a, a plurality of sampling points of each seismic layer are established according to the following model: one of the sampling points is modeled as:
wherein the content of the first and second substances,the zero offset values of the detection force and the moment of the seismic layer sampling point,is the gravity and gravity moment of the sampling point,andthe external force and external moment applied to the sampling point are shown, wherein the moment is the moment of the sampling point relative to the central line of the seismic layer cavity.
Setting the vector from the origin of the central coordinate system of the seismic layer cavity to the centroid of the sampling point asThen
The detailed expression of equation (1) in the survey coordinate system based on the center of the seismic layer void is:
wherein the content of the first and second substances,is a transformation matrix from the world coordinate system to the central coordinate system of the seismic interval cavity, determined by the initial position of the bytes, xc,yc,zc。
When the sampling point is positioned on the upper surface layer of the seismic layer cavity:
detecting the existence relation between force and external force:
F1representing the detected value of force, m, at the position of the upper layerloadRepresenting the real-time gravity of the geology itself, e.g. the weight of the sand, gravel, sand itself, fextExternal force, F0x,F0y,F0zIndicating the resultant force, T, in each direction of the position0x,T0y,T0zIndicating the resultant moment, x, in each direction of the positionc,yc,zcRepresenting the coordinates at that location, the location of the sample point can be uniquely determined.
When the sampling point is positioned on the lower surface of the seismic layer cavity:
detecting the existence relation between force and external force:
when the sampling point is positioned in the middle layer of the seismic layer cavity:
detecting the existence relation between force and external force:
through collecting the force detection measurement data in the motion process, the force detection zero offset and the self load weight parameter to be identified can be synchronously obtained:
finally determine x by the above equation (9)c,yc,zcAnd the coordinates of the positions can uniquely determine the stress and the moment of the sampling points and the positions of the corresponding sampling points.
The calibration method of the embodiment of the invention comprises the following steps: the identification time is short, only force detection data need to be collected, and a simulated exploration-based earthquake stress model is determined by adopting detection force, moment, stress position and combining external force and gravity information of different soil textures.
B, performing layer-by-layer interpretation on each layer of seismic data, and gridding interpretation result data to obtain a geological frame model based on depth; the method comprises the steps of interpreting and tracking a plurality of sets of reflection positions on seismic data, checking and modifying the quality of the interpreted positions, and performing full three-dimensional tracking interpretation on the positions, so that the position interpretation result is fine, each earthquake has a fine position interpretation result, and particularly some unconformities and standard reflection layers need to be interpreted.
Specifically, on the basis of horizon interpretation, well logging curves of a target interval are analyzed in combination with well drilling layering, curve types such as sound waves, density, gamma, porosity, resistivity, lithology and the like are comprehensively analyzed, electrical characteristics, elasticity and radioactivity characteristics of a reservoir are analyzed, lithology and oil-gas content of the reservoir are analyzed, and geological deposition characteristics of the target interval are analyzed, which is not repeated.
And c, determining a logging curve according to the determined geological framework model, and performing fractal random interpolation on a fractal parameter body of the logging curve to obtain an initial logging curve.
Specifically, there may be various embodiments for calculating the fractal parameters of the logging curve model of a single well and converting the fractal parameters of each seismic trace into the fractal parameter body of the logging curve. Respectively calculating fractal parameters of a logging curve depth domain model of a single well, calibrating the seismic fractal parameters of each seismic channel into a fractal parameter body of the logging curve by methods of normalization, proportion weighting and the like, thus obtaining the fractal parameters of each seismic channel of the logging curve, and then calculating the random component funine of fractal random interpolation of each seismic channel according to the following formula:
k is a calibration coefficient, H is a fractal dimension parameter reflecting details and roughness of the hierarchical domain post-stack seismic data, sigma is the variance of the hierarchical domain post-stack seismic data normal distribution, G is a Gauss random variable obeying the standard normal distribution, Delta X is the distance between a known point and an interpolation point, and a is a one-dimensional coordinate of the interpolation point.
There are various ways to find the variance σ and the fractal parameter H. For example, in one embodiment: on the level domain post-stack seismic data volume or the derivative attribute volume thereof, converting one-dimensional seismic channel data into volume data in the longitudinal direction, sequentially establishing a seismic data cube, namely a subvolume, by taking a certain point in the longitudinal direction as a center and a certain distance in the transverse direction as a radius in a geological frame model constraint range, and normalizing the subvolume data.
And calculating the fractal parameters and the variance of each subvolume, so that the fractal parameters and the variance of the seismic data of each seismic channel can be calculated. In the same way, the fractal dimension parameter and variance of the logging curve to be inverted can be calculated, and further the fractal dimension parameter ratio and variance ratio of the logging curve and the seismic data or the derivative attribute body thereof on the well point in a certain interval range can be calculated.
And d, determining an actual logging curve according to the seismic data, correcting the initial logging curve according to the comparison result, and determining a seismic inversion model according to the corrected curve.
Specifically, in the present embodiment, each point selected is corrected according to the following formula (10);
im=ρ×i0 (10)
wherein imRepresents the instantaneous value of the sample point after correction, rho represents the correction factor, i0Representing the instantaneous value of the sample point; the correction coefficient p is calculated by the following formula (11),
where ρ represents a correction coefficient, i01And i02The instantaneous values of two points at the position of deviation when the deviation occurs are shown, N represents the number of sampling times, and k represents the sampling sequence.
The use of the correction mode ensures the stability of the final value by comparing and correcting the empirical value and the real-time value.
In the specific implementation process, the first-stage reactor,
and respectively extracting each well of all the inversion wells as a verification well, performing depth domain reservoir seismic inversion on the unextracted inversion wells, comparing the inversion result of the unextracted inversion wells with the logging curve of the verification well, and calculating the inversion error between the inversion result of the unextracted inversion wells and the logging curve of the verification well.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
Claims (4)
1. A method of seismic inversion in geophysical exploration, comprising:
step a, acquiring logging data in a seismic work area, and acquiring and editing force and moment data of corresponding points of each seismic layer;
b, performing horizon interpretation on each layer of seismic data, and gridding interpretation result data to obtain a geological frame model based on depth;
step c, determining a logging curve according to the determined geological framework model, and performing fractal random interpolation on a fractal parameter body of the logging curve to obtain an initial logging curve;
d, determining an actual logging curve according to the seismic data, correcting the initial logging curve according to a comparison result, and determining a seismic inversion model according to the corrected curve;
in the step a, a plurality of sampling points of each seismic layer are established according to the following model: the establishing model of any sampling point is as follows:
wherein the content of the first and second substances,the zero offset values of the detected force and moment of the sampling points,is the gravity and gravity moment of the sampling point,andexternal force and external moment applied to the sampling point; the moment is the moment of the sampling point relative to the central line of the seismic layer cavity;
the vector from the origin of the measurement coordinate system at the center of the seismic layer cavity to the centroid of the sampling point isWherein x isc,yc,zcRepresenting coordinates; then
The detailed expression of equation (1) in the survey coordinate system based on the center of the seismic layer void is:
wherein the content of the first and second substances,a transformation matrix which is a measurement coordinate system from the world coordinate system to the center of the seismic interval cavity;
when the sampling point is positioned on the upper surface layer of the seismic layer cavity:
detecting the existence relation between force and external force:
F1representing the detected value of force, m, at the position of the upper layerloadRepresenting the real-time gravity of the geological formation itself, fextDenotes an external force, F0x,F0y,F0zIndicating the resultant force, T, in each direction0x,T0y,T0zRepresenting resultant moments in various directions;
when the sampling point is located on the lower surface of the seismic layer cavity:
detecting the existence relation between force and external force:
when the sampling point is positioned in the middle layer of the seismic layer cavity:
detecting the existence relation between force and external force:
by collecting force detection measurement data in the movement process, force detection zero offset and self load weight parameters to be identified are synchronously obtained:
finally determine x by the above equation (9)c,yc,zcAnd acquiring the stress and the moment of the sampling point and the position of the corresponding sampling point by the coordinates at the position.
2. The seismic inversion method in geophysical exploration according to claim 1, wherein in said step d, each point selected is corrected according to the following formula (10);
im=ρ×i0 (10)
wherein imRepresents the instantaneous value of the sample point after correction, rho represents the correction factor, i0Representing the instantaneous value of the sample point; the correction coefficient p is calculated by the following formula (11),
where ρ represents a correction coefficient, i01And i02The instantaneous values of two points at the position of deviation when the deviation occurs are shown, N represents the number of sampling times, and k represents the sampling sequence.
3. The method of seismic inversion in geophysical exploration according to claim 1, wherein in step b, interpretation and tracking of sets of reflection horizons are performed on the seismic data; and checking and modifying the quality of the interpreted horizon, and performing full three-dimensional tracking interpretation on the horizon to ensure that the horizon interpretation result is fine, so that each seismic channel has a fine horizon interpretation result.
4. The seismic inversion method in geophysical exploration according to claim 1, wherein in step c, a single-well log model fractal parameter is calculated, and the fractal parameter of each seismic trace is converted into a log fractal parameter body; respectively calculating fractal parameters of a logging curve depth domain model of a single well, calibrating the seismic fractal parameters of each seismic channel into a fractal parameter body of the logging curve by a normalization and proportion weighting method, further obtaining the fractal parameters of each seismic channel of the logging curve, and then calculating the random component funine of fractal random interpolation of each seismic channel according to the following formula:
k is a calibration coefficient, H is a fractal dimension parameter reflecting details and roughness of the hierarchical domain post-stack seismic data, sigma is the variance of the hierarchical domain post-stack seismic data normal distribution, G is a Gauss random variable obeying the standard normal distribution, Delta X is the distance between a known point and an interpolation point, and a is a one-dimensional coordinate of the interpolation point.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810903194.4A CN108957554B (en) | 2018-08-09 | 2018-08-09 | Seismic inversion method in geophysical exploration |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810903194.4A CN108957554B (en) | 2018-08-09 | 2018-08-09 | Seismic inversion method in geophysical exploration |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108957554A CN108957554A (en) | 2018-12-07 |
CN108957554B true CN108957554B (en) | 2021-03-12 |
Family
ID=64468930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810903194.4A Active CN108957554B (en) | 2018-08-09 | 2018-08-09 | Seismic inversion method in geophysical exploration |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108957554B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111624652B (en) * | 2020-04-17 | 2023-04-25 | 中国石油天然气股份有限公司 | Method for determining geological parameters in a formation |
CN112147709B (en) * | 2020-08-03 | 2022-07-29 | 中国海洋石油集团有限公司 | Gravity gradient data three-dimensional inversion method based on partial smoothness constraint |
CN115434688B (en) * | 2022-08-16 | 2024-01-30 | 成都捷科思石油天然气技术发展有限公司 | Drilling curve control method for logging while drilling of horizontal well |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2369566A1 (en) * | 1999-04-02 | 2000-10-12 | Alan Royce Huffman | A method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
CN103490511B (en) * | 2013-09-13 | 2015-08-05 | 北京师范大学 | A kind of power distribution network communication terminal detection system and method |
CN104570067B (en) * | 2013-10-11 | 2017-06-20 | 中国石油化工股份有限公司 | Phased seismic inversion method in a kind of geophysical exploration |
CN104850467A (en) * | 2015-05-23 | 2015-08-19 | 许昌学院 | Computer self-protection system and computer self-protection method |
CN104960170A (en) * | 2015-05-29 | 2015-10-07 | 张东海 | Intelligent injection molding machine and injection molding method thereof |
CN107329171A (en) * | 2017-06-07 | 2017-11-07 | 中国石油天然气股份有限公司 | Depth Domain reservoir seismic inversion method and device |
-
2018
- 2018-08-09 CN CN201810903194.4A patent/CN108957554B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN108957554A (en) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7062072B2 (en) | Methods of producing images of underground formations surrounding a borehole | |
Gringarten et al. | Teacher's aide variogram interpretation and modeling | |
CN101086535B (en) | A process and program for characterising evolution of an oil reservoir over time | |
CN107329171A (en) | Depth Domain reservoir seismic inversion method and device | |
CN106405651B (en) | Full waveform inversion initial velocity model construction method based on logging matching | |
CN104516018B (en) | Porosity inversion method under lithological constraint in geophysical exploration | |
CN105277978B (en) | A kind of method and device for determining near-surface velocity model | |
Dost et al. | Probabilistic moment tensor inversion for hydrocarbon‐induced seismicity in the Groningen gas field, the Netherlands, part 2: Application | |
US10732310B2 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
US20220291418A1 (en) | An integrated geomechanics model for predicting hydrocarbon and migration pathways | |
CN108957554B (en) | Seismic inversion method in geophysical exploration | |
CN109884710B (en) | Micro-logging tomography method aiming at excitation well depth design | |
WO2016008105A1 (en) | Post-stack wave impedance inversion method based on cauchy distribution | |
WO2017035104A1 (en) | Velocity model seismic static correction | |
CN104122581B (en) | A kind of poststack sound impedance inversion method | |
CN109799540B (en) | Volcanic rock type uranium deposit magnetic susceptibility inversion method based on geological information constraint | |
CN113031068B (en) | Reflection coefficient accurate base tracking prestack seismic inversion method | |
EP3978961A1 (en) | System and method for quantitative seismic integration modeling workflow | |
CN111308558B (en) | Shale gas horizontal well longitudinal wave time difference correction method | |
US20120099396A1 (en) | System and method for characterization with non-unique solutions of anisotropic velocities | |
CN113608258B (en) | Self-consistent deep learning method for constructing high-resolution wave impedance inversion label | |
Rotimi et al. | Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis | |
CN109031423B (en) | Pre-stack seismic texture analysis method based on gradient co-occurrence matrix | |
CN114861515A (en) | Method, device, equipment and medium for calculating layer speed data volume | |
CN113806674A (en) | Method and device for quantifying longitudinal dimension of ancient river channel, electronic equipment and storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |