CN112305589A - Method and device for imaging depth domain of anisotropic medium - Google Patents
Method and device for imaging depth domain of anisotropic medium Download PDFInfo
- Publication number
- CN112305589A CN112305589A CN202011000068.1A CN202011000068A CN112305589A CN 112305589 A CN112305589 A CN 112305589A CN 202011000068 A CN202011000068 A CN 202011000068A CN 112305589 A CN112305589 A CN 112305589A
- Authority
- CN
- China
- Prior art keywords
- work area
- parameter
- anisotropic medium
- thickness
- stratum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 73
- 230000005012 migration Effects 0.000 claims abstract description 32
- 238000013508 migration Methods 0.000 claims abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 230000015572 biosynthetic process Effects 0.000 claims description 33
- 238000004590 computer program Methods 0.000 claims description 16
- 239000000126 substance Substances 0.000 claims description 12
- 238000004587 chromatography analysis Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 10
- 230000006870 function Effects 0.000 claims description 7
- 238000004441 surface measurement Methods 0.000 claims description 4
- 239000010410 layer Substances 0.000 description 108
- 238000010586 diagram Methods 0.000 description 16
- 230000008859 change Effects 0.000 description 12
- 239000003550 marker Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- VYMDGNCVAMGZFE-UHFFFAOYSA-N phenylbutazonum Chemical compound O=C1C(CCCC)C(=O)N(C=2C=CC=CC=2)N1C1=CC=CC=C1 VYMDGNCVAMGZFE-UHFFFAOYSA-N 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002344 surface layer Substances 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. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Abstract
The invention discloses a method and a device for imaging depth domains of anisotropic media, wherein the method comprises the following steps: explaining horizon data according to the isotropic depth domain of the work area, and calculating the thickness of the offset stratum of each horizon of the work area; performing interpolation calculation on each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position, and determining a well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum; according to the well seismic error curved surface after each layer of the work area is combined with the stratum thickness, modeling is carried out to obtain an anisotropic medium parameter model of the work area; and carrying out anisotropic prestack depth migration processing on the anisotropic medium parameter model of the work area, and determining anisotropic medium depth domain imaging data of the work area. The method can accurately perform anisotropic medium parameter model modeling by combining the thickness of the offset stratum, thereby improving the imaging precision of the anisotropic medium depth domain.
Description
Technical Field
The invention relates to the technical field of geophysical, in particular to a method and a device for imaging depth domains of anisotropic media.
Background
This section is intended to provide a background or context to the embodiments of the invention that are recited in the claims. The description herein is not admitted to be prior art by inclusion in this section.
Due to the common existence of anisotropy in the underground medium, the depth of a reflecting interface obtained by isotropic prestack depth migration treatment is not consistent with the depth of a real underground interface. To eliminate such borehole seismic errors, processing may be performed using an anisotropic prestack depth migration technique. While the basis of the anisotropic prestack depth migration technique is an anisotropic media parametric model. In a traditional anisotropic medium depth domain imaging method, TTI (oblique transverse Isotropic medium) anisotropic parameter modeling and migration processing are required.
In the traditional anisotropic medium depth domain imaging method, when modeling of the TTI anisotropic parameters is carried out, well seismic errors are mainly taken as the basis, so that the interpolation of the well seismic errors at a well position to the whole work area range is a very important step. However, the method only considers well seismic errors at well positions, and mass production practices show that the anisotropic medium parameter model obtained by inversion of the method is distorted, so that the result of anisotropic medium depth domain imaging is inaccurate.
Disclosure of Invention
The embodiment of the invention provides an anisotropic medium depth domain imaging method, which is used for improving the anisotropic medium depth domain imaging precision and comprises the following steps:
explaining horizon data according to the isotropic depth domain of the work area, and calculating the thickness of the offset stratum of each horizon of the work area;
performing interpolation calculation on each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position, and determining a well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum;
according to the well seismic error curved surface after each layer of the work area is combined with the stratum thickness, modeling is carried out to obtain an anisotropic medium parameter model of the work area;
and carrying out anisotropic prestack depth migration processing on the anisotropic medium parameter model of the work area, and determining anisotropic medium depth domain imaging data of the work area.
The embodiment of the invention also provides a device for imaging the depth domain of the anisotropic medium, which is used for improving the imaging precision of the depth domain of the anisotropic medium and comprises the following components:
the stratum thickness acquisition module is used for explaining horizon data according to the isotropic depth domain of the work area and calculating the thickness of the offset stratum of each horizon of the work area;
the well-seismic error curved surface determining module is used for carrying out interpolation calculation on each layer of the work area according to the offset stratum thickness of each layer of the work area and the well-seismic error at the well position, and determining a well-seismic error curved surface after each layer of the work area is combined with the stratum thickness;
the anisotropic medium parameter model modeling module is used for modeling according to the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum to obtain an anisotropic medium parameter model of the work area;
and the depth migration processing module is used for performing anisotropic prestack depth migration processing on the anisotropic medium parameter model of the work area and determining the anisotropic medium depth domain imaging data of the work area.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the processor implements the method for imaging the depth domain of the anisotropic medium.
Embodiments of the present invention further provide a computer-readable storage medium, which stores a computer program for executing the above method for anisotropic media depth domain imaging.
In the embodiment of the invention, the horizon data are explained according to the isotropic depth domain of the work area, and the offset stratum thickness of each horizon of the work area is calculated; performing interpolation calculation on each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position, and determining a well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum; according to the well seismic error curved surface after each layer of the work area is combined with the stratum thickness, modeling is carried out to obtain an anisotropic medium parameter model of the work area; and performing anisotropic prestack depth migration processing on the anisotropic medium parameter model of the work area, and determining the anisotropic medium depth domain imaging data of the work area, so that the anisotropic medium parameter model can be accurately modeled in combination with the thickness of a migration stratum during the anisotropic medium parameter model modeling, and further the anisotropic medium depth domain imaging precision is improved. Compared with the prior art, the method not only considers the well seismic error at the well position, but also considers the influence of the thickness of the offset stratum on modeling, and avoids the condition that the obtained anisotropic medium parameter model is distorted to cause inaccurate imaging result of the anisotropic medium depth domain.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
FIG. 1 is a schematic diagram illustrating superimposed display of isotropic depth domain interval velocity, isotropic depth domain interpreted horizon data, and well logging zonal data for a work zone crossing a well trajectory in an embodiment of the present disclosure;
FIG. 2 is a diagram illustrating a parameter v extracted by applying a conventional method for depth domain imaging of an anisotropic medium according to an embodiment of the present inventionP0A cross-sectional view of;
FIG. 3 is a flow chart illustrating a method for depth domain imaging of anisotropic media in an embodiment of the invention;
FIG. 4 is a schematic flow chart of an example of a method for depth domain imaging of anisotropic media in an embodiment of the invention;
FIG. 5 is a schematic flow chart of an example of a method for depth domain imaging of anisotropic media in an embodiment of the invention;
FIG. 6 is a schematic diagram of an example of a method of depth domain imaging of an anisotropic medium in an embodiment of the invention;
FIG. 7 is a schematic diagram of an example of a method of depth domain imaging of an anisotropic medium in an embodiment of the invention;
FIG. 8 is a schematic diagram of an example of a method for depth domain imaging of an anisotropic medium in an embodiment of the invention
Fig. 9 is a schematic structural diagram of an anisotropic medium depth domain imaging device in an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
The inventor finds that: in the traditional anisotropic medium depth domain imaging method, when modeling of the TTI anisotropic parameters is carried out, well seismic errors are mainly taken as the basis, so that the interpolation of the well seismic errors at a well position to the whole work area range is a very important step. However, the method only considers well seismic errors at well positions, and mass production practices show that the anisotropic medium parameter model obtained by inversion of the method is distorted, so that the result of anisotropic medium depth domain imaging is inaccurate. As shown in FIG. 1, FIG. 1 shows the interval velocity of the isotropic depth domain of a work zone crossing a well lineAnd (4) superimposed display of the degree, isotropic depth domain interpretation horizon data and the logging hierarchical data. The left ordinate D/(m) in FIG. 1 is the subsurface depth; the abscissa Xline is a tie line number indicating a change in the lateral position; color code v on the rightP0/(m·s-1) Is the longitudinal wave vertical velocity; t is3x1、T1j22、T1f1、P2l and P1l is the name of geological horizon, which is the bottom of the xu family river group, the bottom of the two sections of the Jialing river group, the bottom boundary of the Feixian group, the bottom boundary of the Longtan group and the bottom boundary of the Liangshan group. The first layer of the well log hierarchical data (labeled T) can be seen3x1Dotted lines in the word) and the floating reference plane (solid lines marked with the floating reference plane), but due to the formation tilt up, the thinnest part of the layer is only 158m (as indicated by the black arrows). FIG. 2 is a diagram illustrating a parameter v extracted by applying a conventional method for depth domain imaging of an anisotropic medium according to an embodiment of the present inventionP0The left ordinate D/(m) in fig. 2 is the subsurface depth; the abscissa Xline is a tie line number indicating a change in the lateral position; color code v on the rightP0/(m·s-1) Is the longitudinal wave vertical velocity; t is3x1、T1j22、T1f1、P2l and P1l is the name of geological horizon, which is the bottom of the xu family river group, the bottom of the two sections of the Jialing river group, the bottom boundary of the Feixian group, the bottom boundary of the Longtan group and the bottom boundary of the Liangshan group. It can be seen that v of the first layer is thinner at the formationP0The parameters were abnormal, with a minimum of only 1200m/s (as indicated by the black arrows), and vP0The transverse change of the anisotropic medium is obviously related to the layer thickness, so that the anisotropic medium parameter model obtained by inversion of the traditional method can be proved to be distorted, and the result of anisotropic medium depth domain imaging is inaccurate.
The embodiment of the invention provides an anisotropic medium depth domain imaging method, which can be used for improving the anisotropic medium depth domain imaging precision, and as shown in fig. 3, the method comprises the following steps:
step 301: explaining horizon data according to the isotropic depth domain of the work area, and calculating the thickness of the offset stratum of each horizon of the work area;
step 302: performing interpolation calculation on each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position, and determining a well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum;
step 303: according to the well seismic error curved surface after each layer of the work area is combined with the stratum thickness, modeling is carried out to obtain an anisotropic medium parameter model of the work area;
step 304: and carrying out anisotropic prestack depth migration processing on the anisotropic medium parameter model of the work area, and determining anisotropic medium depth domain imaging data of the work area.
In the embodiment of the invention, the horizon data are explained according to the isotropic depth domain of the work area, and the offset stratum thickness of each horizon of the work area is calculated; performing interpolation calculation on each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position, and determining a well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum; according to the well seismic error curved surface after each layer of the work area is combined with the stratum thickness, modeling is carried out to obtain an anisotropic medium parameter model of the work area; and performing anisotropic prestack depth migration processing on the anisotropic medium parameter model of the work area, and determining the anisotropic medium depth domain imaging data of the work area, so that the anisotropic medium parameter model can be accurately modeled in combination with the thickness of a migration stratum during the anisotropic medium parameter model modeling, and further the anisotropic medium depth domain imaging precision is improved. Compared with the prior art, the method not only considers the well seismic error at the well position, but also considers the influence of the thickness of the offset stratum on modeling, and avoids the condition that the obtained anisotropic medium parameter model is distorted to cause inaccurate imaging result of the anisotropic medium depth domain.
In specific implementation, the horizon data are explained according to the isotropic depth domain of the work area, and the offset stratum thickness of each horizon of the work area is calculated.
In an embodiment, the method for calculating the thickness of the offset stratum of each horizon of the work area according to the isotropic depth domain interpretation horizon data of the work area includes a plurality of methods, for example, the method may include: explaining horizon data according to an isotropic depth domain of the work area, and determining the offset stratum depth of each horizon of the work area; and calculating the thickness of the offset stratum of each layer of the work area according to the depth of the offset stratum of each layer of the work area.
In the above embodiments, there are various methods for calculating the thickness of the offset formation of each layer of the work area according to the depth of the offset formation of each layer of the work area, for example, the method may include: the thickness of the offset stratum of each layer of the work area can be calculated according to the offset stratum depth of each layer of the work area according to the following formula:
wherein the content of the first and second substances,
the thickness of the offset stratum of the ith layer of the work area;
the offset stratum depth of the ith horizon of the work area is, when i is equal to 1,
in the embodiment, the calculation of the thickness of the offset stratum of each layer of the work area is beneficial to the modeling of the anisotropic medium parameter model in the subsequent steps.
In specific implementation, after the offset stratum thickness of each layer of the work area is calculated, interpolation calculation is carried out on each layer of the work area according to the offset stratum thickness of each layer of the work area and well seismic errors at well positions, and a well seismic error curved surface formed by combining each layer of the work area with the stratum thickness is determined.
In one embodiment, the method for determining the well-seismic error surface after each layer of the work area is combined with the thickness of the formation may include: determining the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position; and performing interpolation calculation on the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area, and determining the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum.
In the embodiment, the information of the formation thickness can be introduced into the flow of anisotropic parameter inversion, so that the inversion result is more accurate. The inventor finds that the well-seismic error interpolation method in the traditional technology does not consider the influence of the change of the formation thickness on the extraction result of the Thomsen parameters of the anisotropic medium. The interpolation of the well seismic error is essentially to solve the trend of the change of anisotropic parameters, and v is obtained by inversion under the condition of severe change of the thickness of the stratumP0Distortion of the parameters and delta parameters, vP0The values of the parameters and the delta parameters are obviously related to the change of the thickness of the stratum. Therefore, the value of the Thomsen parameter of the anisotropic medium is introduced into the interpolation process in the embodiment of the invention.
In the embodiment, the well-seismic error curved surface after each layer of the work area is combined with the stratum thickness is determined, anisotropic medium parameter model modeling in the subsequent steps is facilitated, compared with the prior art, the influence of the stratum thickness on the subsequent anisotropic medium parameter model modeling value can be avoided, the distortion of the inversion result is avoided, and meanwhile, the value of the anisotropic parameter model obtained through inversion is more accurate and reliable, so that the stratum attitude of the offset section is more reliable.
In specific implementation, the method for imaging the depth domain of the anisotropic medium provided by the embodiment of the present invention may further include: and determining the well seismic error of each layer of the work area at each well position after combining the thickness of the stratum according to the well seismic error curved surface after combining each layer of the work area with the thickness of the stratum. In an embodiment, there are various methods for determining the ratio parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium for each layer of the work area, for example, the method may include: determining the proportion parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area according to the following formula:
wherein the content of the first and second substances,
a proportional parameter of a delta parameter of the ith horizon of the work area based on the Thomsen parameter of the anisotropic medium; mi(xw,yw) The well seismic error is the well seismic error at the w well position of the i position of the work area;
when the well seismic error at the w-th well position after the thickness of the stratum is combined at the i-1 th position of the work area, i is 1,
is the thickness of the offset formation at the ith horizon of the work area.
In an embodiment, the proportional parameter at w well position of delta parameter of Thomsen parameter based on anisotropic medium at i layer of work area can be calculated firstly
And the proportional parameter at the w well position of the delta parameter of the Thomsen parameter based on the anisotropic medium of the i layer position of the work area
Interpolating and extrapolating to the full work area to obtain the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of the ith layer of the work area
In the embodiment, the well seismic error curved surface of each layer of the work area combined with the thickness of the stratum is determined according to the following formula:
wherein the content of the first and second substances,
the surface is a well seismic error curved surface of the ith position of the work area after the thickness of the stratum is combined;
a proportional parameter of a delta parameter of the ith horizon of the work area based on the Thomsen parameter of the anisotropic medium;
the curved surface is the well seismic error curved surface after the thickness of the stratum is combined at the i-1 th position of the work area, when i is 1,
is the thickness of the offset formation at the ith horizon of the work area.
In the above embodiment, the calculation formula of the well-seismic error surface after each layer of the work area is combined with the thickness of the formation can be derived as follows:
(1) determining the following formula as an estimation formula of the delta parameter of the Thomsen parameter of the anisotropic medium:
wherein the content of the first and second substances,
the thickness of the offset stratum of the ith layer of the work area;
the actual stratum thickness of the ith horizon of the work area;
(2) then, acquiring logging hierarchical data of the work area; and determining the logging layering depth of each layer of the work area according to the logging layering data of the work area. And the actual formation thickness may be determined by the following equation:
wherein the content of the first and second substances,
logging layered depth of the ith layer of the work area;
the actual stratum thickness of the ith horizon of the work area;
(3) because the logging layering depth has the following relation with the thickness of the offset stratum:
wherein the content of the first and second substances,
the offset stratum depth of the ith horizon of the work area is, when i is equal to 1,
the surface is a well seismic error curved surface of the ith position of the work area after the thickness of the stratum is combined;
(4) substituting the formulas in the steps (2) and (3) into the formula in the step (1) to obtain:
from the above calculation results, it can be seen that the change of the delta parameter with position is only related to
It is related. Embodiments of the present invention may first determine the ratio of the delta parameters of the Thomsen parameters based on anisotropic media for each horizon of the work areaExample parameters
Then to the ratio parameter
Interpolation is carried out, and finally, a well seismic error curved surface (Mitie curved surface) of the ith position of the work area after the thickness of the stratum is combined is calculated by the following formula
Wherein, when i is 1,
obtained by the method based on the embodiment of the invention
Thickness of offset formation from ith horizon
Correlation, the inverted anisotropic medium parameter model is not affected by the formation thickness variation. When i > 1, in
Performing interpolation calculation on the object to obtain the final product
Both with the offset formation thickness of the corresponding horizon
Correlation, also with the Mitie information of the previous layer
And (4) correlating. Therefore, according to the method provided by the embodiment of the invention, the influence of the thickness change of the stratum can be effectively avoided from shallow to deep. Fig. 4 is a schematic flow chart of an example of a method for imaging an anisotropic medium depth domain in an embodiment of the present invention, which shows a flow of interpolation calculation, and the general idea is to obtain a mixture curved surface layer by layer from top to bottom, and input the mixture curved surface of the layer into the obtaining process of the next layer. Introducing this flow into the subsequent modeling flow may form an optimized anisotropic parametric modeling flow:
(1) firstly, obtaining logging hierarchical data and isotropic depth domain interpretation horizon data of a work area, and when a horizon serial number i is 1, using logging hierarchical information at a well position in the logging hierarchical data of the work area and a depth domain interpretation horizon in the isotropic depth domain interpretation horizon data to calculate a well seismic error at a w well position after the thickness of a combined stratum of a 1 st horizon of the work area
(2) Combining the 1 st horizon of the work area with the borehole seismic error at the w-th horizon
Offset formation thickness at w-well location from work zone 1-level
Dividing to obtain a proportional parameter of a delta parameter of the Thomsen parameter based on the anisotropic medium at the w well position of the 1 st layer of the work area
(3) Proportional parameter of delta parameter of Thomsen parameter based on anisotropic medium at w well position of 1 st layer of work area
Extrapolating the interpolation to the full work area to obtain the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of the 1 st layer of the work area
(4) Using the formula
Well seismic error curved surface after calculating the thickness of the stratum combined with the 1 st layer of the work area
(5) And processing the next layer of bits of the work area, namely i-i + 1. The well-seismic error curved surface after the combination of the ith layer of the work area and the stratum thickness is solved by using the well logging layered information at the well position in the well logging layered data of the work area and the depth domain interpretation layer in the isotropic depth domain interpretation layer data
(6) Calculating the well seismic error M at the w well position of the i position of the work areai(xw,yw) And the well seismic error at the w well position after the combination of the i-1 position of the work area and the thickness of the stratum
Offset formation thickness at well location from i-th horizon of work area
Substituting into formula
Obtaining a proportional parameter at the w well position of the delta parameter of the Thomsen parameter based on the anisotropic medium at the i layer position of the work area
(7) Proportional parameter at w well position of delta parameter of Thomsen parameter based on anisotropic medium of i layer position of work area
Interpolating and extrapolating to the full work area to obtain the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of the ith layer of the work area
(8) Using the formula
Well seismic error curved surface after calculating thickness of combination stratum of ith position of work area
(9) And (4) judging whether the ith layer of the work area is the last layer of the stratum of the work area, if so, ending the process, otherwise, skipping to the step (5) to continue the calculation.
In the embodiment, the well-seismic error curved surface of each layer of the work area can be accurately calculated by combining the thickness of the offset stratum, so that the imaging precision of the depth domain of the anisotropic medium obtained by subsequent modeling is improved. Compared with the prior art, the method not only considers the well seismic error at the well position, but also considers the influence of the thickness of the offset stratum on modeling, and avoids the distortion of the obtained anisotropic medium parameter model.
In specific implementation, after determining the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum, the modeling is performed according to the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum to obtain an anisotropic medium parameter model of the work area, which may include: acquiring the isotropic depth domain speed of a work area; extracting a stratigraphic dip angle and an azimuth angle from the isotropic depth domain interpretation horizon data of the work area to obtain dip angle and azimuth angle attribute data of the work area; according to the dip angle and azimuth angle attribute data of the work area, explaining layer data of each layer by combining a well seismic error curved surface, an isotropic depth domain speed and an isotropic depth domain after the formation thickness, and establishing a grid chromatography inversion matrix of the work area; carrying out grid chromatography inversion on the established grid chromatography inversion matrix to obtain the anisotropic medium of the work areaA parametric model; the anisotropic medium parameter model comprises the following parameters: v. ofP0Parameter, delta parameter and epsilon parameter, where vP0The parameter is the vertical velocity of longitudinal wave, and the delta parameter and the epsilon parameter are Thomsen parameters of the anisotropic medium.
In the embodiment, the position data is interpreted by each layer according to the dip angle and azimuth angle attribute data of the work area, the well seismic error curved surface after the thickness of the stratum is combined, the isotropic depth domain speed and the isotropic depth domain to obtain the anisotropic medium parameter model of the work area, the anisotropic medium parameter model can be accurately modeled, and the imaging precision of the anisotropic medium depth domain is further improved.
In an embodiment, the method for imaging an anisotropic medium depth domain provided in the embodiment of the present invention may further include: acquiring logging hierarchical data of a work area; v in the parametric model of anisotropic media as a function of work areaP0Calculating anisotropic medium depth domain interpretation horizon data of the work area according to the parameters and the delta parameters; carrying out error comparison on the anisotropic medium depth domain interpretation horizon data of the work area and the logging hierarchical data; and when the anisotropic medium depth domain interpreted horizon data has errors with the logging layering data, replacing the isotropic depth domain interpreted horizon data of the work area with the anisotropic medium depth domain interpreted horizon data generated by calculation, re-modeling to obtain an anisotropic medium parameter model of the work area, and repeating the steps until the calculated anisotropic medium depth domain interpreted horizon data has no errors with the logging layering data according to the anisotropic medium parameter model of the work area obtained by re-modeling.
In the embodiment, the anisotropic medium parameter model of the work area can be accurately obtained through error comparison of the anisotropic medium depth domain interpretation horizon data of the work area and the logging hierarchical data.
In specific implementation, the anisotropic pre-stack depth migration processing is performed on the anisotropic medium parameter model of the work area, and the determining of the anisotropic medium depth domain imaging data of the work area may include: anisotropic prestack depth migration processing is carried out on the anisotropic medium parameter model, the inclination angle and the azimuth angle attribute data of the work area, and a common imaging point gather is obtained; and determining the anisotropic medium depth domain imaging data of the work area according to the common imaging point gather.
In an embodiment, the method for imaging an anisotropic medium depth domain provided in the embodiment of the present invention may further include: when the in-phase axis in the common imaging point gather is not straight, picking up residual delay information from the common imaging point gather; performing chromatographic inversion on the residual delay information, the anisotropic medium depth domain interpretation horizon data and the anisotropic medium parameter model to obtain an epsilon parameter; and replacing the epsilon parameter in the anisotropic medium parameter model by the obtained epsilon parameter, replacing the anisotropic medium parameter model with the epsilon parameter, performing anisotropic prestack depth migration again to obtain a common imaging point gather, and repeatedly executing the steps until the in-phase axis in the common imaging point gather is straight.
In the embodiment, the common imaging point gather with the straight in-phase axis can be finally obtained by repeatedly executing the steps, so that the anisotropic medium depth domain imaging data of the work area can be determined, and the anisotropic medium depth domain imaging precision is improved.
Based on the above embodiments, the present invention provides the following embodiment to explain a detailed flow of the method for imaging the depth domain of the anisotropic medium, fig. 5 is a schematic flow chart of an example of the method for imaging the depth domain of the anisotropic medium in the embodiment of the present invention, as shown in fig. 5, generally, a first inversion v is adoptedP0And delta and independently inverting the epsilon, wherein the detailed steps are as follows:
(1) firstly, inputting isotropic depth domain interpretation horizon data (depth domain interpretation horizon for short), and extracting a formation dip angle and an azimuth angle from the isotropic depth domain interpretation horizon data to obtain dip angle and azimuth angle attribute data volumes (dip angle and azimuth angle volumes for short);
(2) inputting logging hierarchical data (which can be called a well Marker value), and solving a depth error (well point Mitie) between the logging hierarchical data and a depth domain interpretation horizon so as to determine a well seismic error at the well location;
(3) according to the well seismic errors at the well positions, performing layer-by-layer interpolation calculation on each layer position of the work area, and solving a well seismic error curved surface (such as a Mitie curved surface in the figure) of each layer position in the whole work area range;
(4) collecting a well seismic error surface, a depth domain interpretation horizon and pre-acquired isotropic depth domain speed data to a database, and combining the data with an inclination angle body and an azimuth angle body to establish a grid chromatography inversion matrix;
(5) under the preset Thomsen parameter relationship of the anisotropic medium, if epsilon is assumed to be 2 delta, obtaining an anisotropic medium parameter model through grid chromatography inversion according to the established grid chromatography inversion matrix, wherein the anisotropic medium parameter model comprises the following parameters: v. ofP0δ and ε, wherein vP0Is the vertical velocity of longitudinal wave, delta and epsilon are Thomsen parameters of anisotropic medium;
(6) according to vP0And delta, obtaining the calculated anisotropic medium depth domain interpreted horizon data, replacing the isotropic depth domain interpreted horizon data in the updating step (1) with the calculated anisotropic medium depth domain interpreted horizon data, and comparing the calculated anisotropic medium depth domain interpreted horizon data with the input well Marker value; and if the comparison result indicates that the calculated anisotropic medium depth domain interpreted horizon data has errors with the input well Marker value, the steps are carried out again to obtain the calculated anisotropic medium depth domain interpreted horizon data for iterative calculation until the comparison result indicates that the calculated anisotropic medium depth domain interpreted horizon data has no errors with the input well Marker value, namely the well information is completely matched with the seismic stratification. After determining that the anisotropic medium depth domain interpretation horizon data has no error with the input well Marker value, utilizing the dip angle, the azimuth angle body and the vP0And delta and epsilon are subjected to TTI anisotropic prestack depth migration to obtain a Common Imaging Point (CIP) gather. When the in-phase axis in the CIP gather is determined to be not straight, the operation in the step (7) is carried out; when the homophase axis in the CIP gather is determined to be straight, the operation in the step (8) is carried out;
(7) picking up residual delay along layers on the CIP gather, inputting residual delay information, a depth domain interpretation layer and an anisotropic medium parameter model into a chromatography inversion program, and determining epsilon; replacing epsilon in the anisotropic medium parameter model determined in the step by the determined epsilon, carrying out TTI anisotropic depth migration again by using the replaced anisotropic medium parameter model to obtain a Common Imaging Point (CIP) gather, and repeating the steps until a phase axis in the CIP gather is straight; then carrying out step (8);
(8) and outputting the anisotropic medium depth domain imaging profile according to the obtained Common Imaging Point (CIP) gather with the straight in-phase axis.
A specific example is given below to illustrate the specific application of the method proposed by the present invention.
FIG. 6 is a schematic diagram of an example of a method for imaging an anisotropic medium depth domain according to an embodiment of the present invention, in which v in a parametric model of an anisotropic medium obtained by inversion based on the isotropic velocity model of FIG. 1 using the method for imaging an anisotropic medium depth domain according to an embodiment of the present inventionP0And (4) parameters. The left ordinate D/(m) in FIG. 6 is the subsurface depth; the abscissa Xline is a tie line number indicating a change in the lateral position; color code v on the rightP0/(m·s-1) Is the longitudinal wave vertical velocity; t is3x1、T1j22、T1f1、P2l and P1l is the name of geological horizon, which is the bottom of the xu family river group, the bottom of the two sections of the Jialing river group, the bottom boundary of the Feixian group, the bottom boundary of the Longtan group and the bottom boundary of the Liangshan group. From fig. 6 it can be seen that the first layer velocity distortion in fig. 2 has been eliminated. In addition, the velocity in the layer is more uniform. At this point, the change in intrazone velocity is related only to the isotropic zone velocity trend and the interwell seismic error trend, and not to the thickness of the zone and overburden. Therefore, the model obtained by applying the method for imaging the depth domain of the anisotropic medium provided by the embodiment of the invention is more accurate and reliable. Fig. 7 and 8 are schematic flow charts of an example of a method for imaging an anisotropic medium depth domain in an embodiment of the present invention, and fig. 7 and 8 respectively use a conventional method for imaging an anisotropic medium depth domain and a method for imaging an anisotropic medium depth domain provided by an embodiment of the present invention to form a contrast of an offset profile. In FIGS. 7 and 8The ordinate D/(m) is the underground depth; the abscissa Xline is a tie line number indicating a change in the lateral position; t is3x1、T1j22、T1f1、P2l and P1l is the name of geological horizon, which is the bottom of the xu family river group, the bottom of the two sections of the Jialing river group, the bottom boundary of the Feixian group, the bottom boundary of the Longtan group and the bottom boundary of the Liangshan group. The arrows in fig. 7 indicate the locations where the distortion in fig. 1 occurs. The anisotropic medium parameter model obtained by the traditional anisotropic medium depth domain imaging method is easily influenced by the variation trend of the thickness of the stratum, and although the offset profile of the model is matched with a drilled well at a well position, the stratum attitude is not reasonable. As shown in fig. 8, by applying the method for anisotropic medium depth domain imaging provided by the embodiment of the present invention, an interpolation method based on a proportional parameter of a delta parameter of a Thomsen parameter of an anisotropic medium is used when an anisotropic parameter model is obtained, so that the influence of the layer thickness on the value of the anisotropic medium parameter model is effectively avoided, and the formation shape of the offset profile is more reliable.
The embodiment of the invention also provides a device for imaging the depth domain of the anisotropic medium, which is described in the following embodiment. Because the principle of the device for solving the problems is similar to the method for imaging the depth domain of the anisotropic medium, the implementation of the device can refer to the implementation of the method for imaging the depth domain of the anisotropic medium, and repeated details are not repeated.
The device for imaging the depth domain of the anisotropic medium provided by the embodiment of the invention can be used for improving the imaging precision of the depth domain of the anisotropic medium, and as shown in fig. 9, the device can comprise:
the stratum thickness obtaining module 01 is used for explaining horizon data according to the isotropic depth domain of the work area and calculating the thickness of the offset stratum of each horizon of the work area;
the well-seismic error curved surface determining module 02 is used for performing interpolation calculation on each layer of the work area according to the offset stratum thickness of each layer of the work area and the well-seismic error at the well position, and determining a well-seismic error curved surface after each layer of the work area is combined with the stratum thickness;
the anisotropic medium parameter model modeling module 03 is used for modeling according to the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum to obtain an anisotropic medium parameter model of the work area;
and the depth migration processing module 04 is configured to perform anisotropic prestack depth migration processing on the anisotropic medium parameter model of the work area, and determine anisotropic medium depth domain imaging data of the work area.
In one embodiment, the formation thickness acquisition module is specifically operable to: explaining horizon data according to an isotropic depth domain of the work area, and determining the offset stratum depth of each horizon of the work area; and calculating the thickness of the offset stratum of each layer of the work area according to the depth of the offset stratum of each layer of the work area.
In one embodiment, the offset stratigraphic thickness for each horizon of the work area may be calculated from the offset stratigraphic depth for each horizon of the work area as follows:
wherein the content of the first and second substances,
the thickness of the offset stratum of the ith layer of the work area;
the offset stratum depth of the ith horizon of the work area is, when i is equal to 0,
in one embodiment, the borehole error surface determination module is specifically operable to: determining the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position; and performing interpolation calculation on the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area, and determining the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum.
In one embodiment, the apparatus for anisotropic media depth domain imaging provided by the embodiment of the present invention may further include: a borehole seismic error determination module at each well location for: determining well seismic errors at each well position after each layer of the work area is combined with the thickness of the stratum according to the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum; the well-seismic error surface determination module can be specifically used for: determining the proportion parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area according to the following formula:
wherein the content of the first and second substances,
a proportional parameter of a delta parameter of the ith horizon of the work area based on the Thomsen parameter of the anisotropic medium; mi(xw,yw) The well seismic error is the well seismic error at the w well position of the i position of the work area;
the well seismic error at the w well position after the thickness of the stratum is combined at the i-1 position of the work area, when i is 1, M isi(xw,yw)＝0；
Is the thickness of the offset formation at the ith horizon of the work area.
In one embodiment, the borehole error surface determination module is specifically operable to: determining a well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum according to the following formula:
wherein the content of the first and second substances,
junction at i-th layer of work areaCombining the well seismic error curved surface after the formation thickness;
a proportional parameter of a delta parameter of the ith horizon of the work area based on the Thomsen parameter of the anisotropic medium;
the curved surface is the well seismic error curved surface after the thickness of the stratum is combined at the i-1 th position of the work area, when i is 1,
is the thickness of the offset formation at the ith horizon of the work area.
In one embodiment, the anisotropic media parametric model modeling module is specifically configured to: acquiring the isotropic depth domain speed of a work area; extracting a stratigraphic dip angle and an azimuth angle from the isotropic depth domain interpretation horizon data of the work area to obtain dip angle and azimuth angle attribute data of the work area; according to the dip angle and azimuth angle attribute data of the work area, explaining layer data of each layer by combining a well seismic error curved surface, an isotropic depth domain speed and an isotropic depth domain after the formation thickness, and establishing a grid chromatography inversion matrix of the work area; carrying out grid chromatographic inversion on the established grid chromatographic inversion matrix to obtain an anisotropic medium parameter model of the work area; the anisotropic medium parameter model comprises the following parameters: v. ofP0Parameter, delta parameter and epsilon parameter, where vP0The parameter is the vertical velocity of longitudinal wave, and the delta parameter and the epsilon parameter are Thomsen parameters of the anisotropic medium.
In one embodiment, the apparatus for anisotropic media depth domain imaging provided by the embodiment of the present invention may further include: an error comparison module to: acquiring logging hierarchical data of a work area; v in the parametric model of anisotropic media as a function of work areaP0Calculating anisotropic medium depth domain interpretation horizon data of the work area according to the parameters and the delta parameters; interpreting the depth domain of the anisotropic medium of the work area to layer numberComparing errors with the logging hierarchical data; and when the anisotropic medium depth domain interpreted horizon data has errors with the logging layering data, replacing the isotropic depth domain interpreted horizon data of the work area with the anisotropic medium depth domain interpreted horizon data generated by calculation, re-modeling to obtain an anisotropic medium parameter model of the work area, and repeating the steps until the calculated anisotropic medium depth domain interpreted horizon data has no errors with the logging layering data according to the anisotropic medium parameter model of the work area obtained by re-modeling.
In one embodiment, the depth offset processing module is specifically configured to: anisotropic prestack depth migration processing is carried out on the anisotropic medium parameter model, the inclination angle and the azimuth angle attribute data of the work area, and a common imaging point gather is obtained; and determining the anisotropic medium depth domain imaging data of the work area according to the common imaging point gather.
In one embodiment, the apparatus for anisotropic media depth domain imaging provided by the embodiment of the present invention may further include: an in-phase axis correction module to: when the in-phase axis in the common imaging point gather is not straight, picking up residual delay information from the common imaging point gather; performing chromatographic inversion on the residual delay information, the anisotropic medium depth domain interpretation horizon data and the anisotropic medium parameter model to obtain an epsilon parameter; and replacing the epsilon parameter in the anisotropic medium parameter model by the obtained epsilon parameter, replacing the anisotropic medium parameter model with the epsilon parameter, performing anisotropic prestack depth migration again to obtain a common imaging point gather, and repeatedly executing the steps until the in-phase axis in the common imaging point gather is straight.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the processor implements the method for imaging the depth domain of the anisotropic medium.
Embodiments of the present invention further provide a computer-readable storage medium, which stores a computer program for executing the above method for anisotropic media depth domain imaging.
In the embodiment of the invention, the horizon data are explained according to the isotropic depth domain of the work area, and the offset stratum thickness of each horizon of the work area is calculated; performing interpolation calculation on each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position, and determining a well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum; according to the well seismic error curved surface after each layer of the work area is combined with the stratum thickness, modeling is carried out to obtain an anisotropic medium parameter model of the work area; and performing anisotropic prestack depth migration processing on the anisotropic medium parameter model of the work area, and determining the anisotropic medium depth domain imaging data of the work area, so that the anisotropic medium parameter model can be accurately modeled in combination with the thickness of a migration stratum during the anisotropic medium parameter model modeling, and further the anisotropic medium depth domain imaging precision is improved. Compared with the prior art, the method not only considers the well seismic error at the well position, but also considers the influence of the thickness of the offset stratum on modeling, and avoids the condition that the obtained anisotropic medium parameter model is distorted to cause inaccurate imaging result of the anisotropic medium depth domain.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (17)
1. A method of anisotropic media depth domain imaging, comprising:
explaining horizon data according to the isotropic depth domain of the work area, and calculating the thickness of the offset stratum of each horizon of the work area;
performing interpolation calculation on each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position, and determining a well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum;
according to the well seismic error curved surface after each layer of the work area is combined with the stratum thickness, modeling is carried out to obtain an anisotropic medium parameter model of the work area;
and carrying out anisotropic prestack depth migration processing on the anisotropic medium parameter model of the work area, and determining anisotropic medium depth domain imaging data of the work area.
2. The method of claim 1, wherein computing the offset stratigraphic thickness for each horizon of the work area from the isotropic depth domain interpreted horizon data for the work area comprises:
explaining horizon data according to an isotropic depth domain of the work area, and determining the offset stratum depth of each horizon of the work area;
and calculating the thickness of the offset stratum of each layer of the work area according to the depth of the offset stratum of each layer of the work area.
3. The method of claim 2, wherein the offset formation thickness for each horizon of the work area is calculated from the offset formation depth for each horizon of the work area according to the following formula:
wherein the content of the first and second substances,
the thickness of the offset stratum of the ith layer of the work area;
the offset stratum depth of the ith horizon of the work area is, when i is equal to 1,
4. the method of claim 1, wherein interpolating each horizon of the work area based on the thickness of the offset formation and the borehole error at the borehole location for each horizon of the work area to determine a borehole error surface for each horizon of the work area in combination with the thickness of the formation comprises:
determining the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position;
and performing interpolation calculation on the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area, and determining the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum.
5. The method of claim 4, further comprising: determining the well seismic error at each well position after each layer of the work area is combined with the thickness of the stratum according to the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum;
determining the proportion parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area according to the following formula:
wherein the content of the first and second substances,
a proportional parameter of a delta parameter of the ith horizon of the work area based on the Thomsen parameter of the anisotropic medium; mi(xw,yw) The well seismic error is the well seismic error at the w well position of the i position of the work area;
when the well seismic error at the w-th well position after the thickness of the stratum is combined at the i-1 th position of the work area, i is 1,
is the thickness of the offset formation at the ith horizon of the work area.
6. The method of claim 4, wherein the borehole error surface for each horizon in the work area in combination with the formation thickness is determined as follows:
wherein the content of the first and second substances,
the surface is a well seismic error curved surface of the ith position of the work area after the thickness of the stratum is combined;
a proportional parameter of a delta parameter of the ith horizon of the work area based on the Thomsen parameter of the anisotropic medium;
the curved surface is the well seismic error curved surface after the thickness of the stratum is combined at the i-1 th position of the work area, when i is 1,
is the thickness of the offset formation at the ith horizon of the work area.
7. The method of claim 1, wherein the step of obtaining the anisotropic medium parameter model of the work area by modeling according to the well-seismic error surface after each layer of the work area is combined with the thickness of the stratum comprises:
acquiring the isotropic depth domain speed of a work area;
extracting a stratigraphic dip angle and an azimuth angle from the isotropic depth domain interpretation horizon data of the work area to obtain dip angle and azimuth angle attribute data of the work area;
according to the dip angle and azimuth angle attribute data of the work area, explaining layer data of each layer by combining a well seismic error curved surface, an isotropic depth domain speed and an isotropic depth domain after the formation thickness, and establishing a grid chromatography inversion matrix of the work area;
carrying out grid chromatographic inversion on the established grid chromatographic inversion matrix to obtain an anisotropic medium parameter model of the work area; the anisotropic medium parameter model comprises the following parameters: v. ofP0Parameter, delta parameter and epsilon parameter, where vP0The parameter is the vertical velocity of longitudinal wave, and the delta parameter and the epsilon parameter are Thomsen parameters of the anisotropic medium.
8. The method of claim 7, further comprising:
acquiring logging hierarchical data of a work area;
v in the parametric model of anisotropic media as a function of work areaP0Calculating anisotropic medium depth domain interpretation horizon data of the work area according to the parameters and the delta parameters;
carrying out error comparison on the anisotropic medium depth domain interpretation horizon data of the work area and the logging hierarchical data; and when the anisotropic medium depth domain interpreted horizon data has errors with the logging layering data, replacing the isotropic depth domain interpreted horizon data of the work area with the anisotropic medium depth domain interpreted horizon data generated by calculation, re-modeling to obtain an anisotropic medium parameter model of the work area, and repeating the steps until the calculated anisotropic medium depth domain interpreted horizon data has no errors with the logging layering data according to the anisotropic medium parameter model of the work area obtained by re-modeling.
9. The method of claim 8, wherein performing anisotropic prestack depth migration on the parametric model of the anisotropic medium for the work area to determine anisotropic medium depth domain imaging data for the work area comprises:
anisotropic prestack depth migration processing is carried out on the anisotropic medium parameter model, the inclination angle and the azimuth angle attribute data of the work area, and a common imaging point gather is obtained;
and determining the anisotropic medium depth domain imaging data of the work area according to the common imaging point gather.
10. The method of claim 9, further comprising:
when the in-phase axis in the common imaging point gather is not straight, picking up residual delay information from the common imaging point gather;
performing chromatographic inversion on the residual delay information, the anisotropic medium depth domain interpretation horizon data and the anisotropic medium parameter model to obtain an epsilon parameter;
and replacing the epsilon parameter in the anisotropic medium parameter model by the obtained epsilon parameter, replacing the anisotropic medium parameter model with the epsilon parameter, performing anisotropic prestack depth migration again to obtain a common imaging point gather, and repeatedly executing the steps until the in-phase axis in the common imaging point gather is straight.
11. An apparatus for depth domain imaging of anisotropic media, comprising:
the stratum thickness acquisition module is used for explaining horizon data according to the isotropic depth domain of the work area and calculating the thickness of the offset stratum of each horizon of the work area;
the well-seismic error curved surface determining module is used for carrying out interpolation calculation on each layer of the work area according to the offset stratum thickness of each layer of the work area and the well-seismic error at the well position, and determining a well-seismic error curved surface after each layer of the work area is combined with the stratum thickness;
the anisotropic medium parameter model modeling module is used for modeling according to the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum to obtain an anisotropic medium parameter model of the work area;
and the depth migration processing module is used for performing anisotropic prestack depth migration processing on the anisotropic medium parameter model of the work area and determining the anisotropic medium depth domain imaging data of the work area.
12. The apparatus of claim 11, wherein the formation thickness acquisition module is specifically configured to:
explaining horizon data according to an isotropic depth domain of the work area, and determining the offset stratum depth of each horizon of the work area;
and calculating the thickness of the offset stratum of each layer of the work area according to the depth of the offset stratum of each layer of the work area.
13. The apparatus of claim 11, wherein the borehole error surface determination module is specifically configured to: determining the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area according to the thickness of the offset stratum of each layer of the work area and the well seismic error at the well position;
and performing interpolation calculation on the proportional parameter of the delta parameter of the Thomsen parameter based on the anisotropic medium of each layer of the work area, and determining the well seismic error curved surface after each layer of the work area is combined with the thickness of the stratum.
14. The apparatus of claim 11, wherein the anisotropic media parametric model modeling module is specifically configured to:
acquiring the isotropic depth domain speed of a work area;
extracting a stratigraphic dip angle and an azimuth angle from the isotropic depth domain interpretation horizon data of the work area to obtain dip angle and azimuth angle attribute data of the work area;
according to the dip angle and azimuth angle attribute data of the work area, explaining layer data of each layer by combining a well seismic error curved surface, an isotropic depth domain speed and an isotropic depth domain after the formation thickness, and establishing a grid chromatography inversion matrix of the work area;
carrying out grid chromatographic inversion on the established grid chromatographic inversion matrix to obtain an anisotropic medium parameter model of the work area; the anisotropic medium parameter model comprises the following parameters: v. ofP0Parameter, delta parameter and epsilon parameter, where vP0The parameter being vertical velocity of longitudinal waveThe degree, delta parameter and epsilon parameter are the Thomsen parameters of the anisotropic medium.
15. The apparatus of claim 14, wherein the depth migration processing module is specifically configured to:
anisotropic prestack depth migration processing is carried out on the anisotropic medium parameter model, the inclination angle and the azimuth angle attribute data of the work area, and a common imaging point gather is obtained;
and determining the anisotropic medium depth domain imaging data of the work area according to the common imaging point gather.
16. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 10 when executing the computer program.
17. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 10.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011000068.1A CN112305589A (en) | 2020-09-22 | 2020-09-22 | Method and device for imaging depth domain of anisotropic medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011000068.1A CN112305589A (en) | 2020-09-22 | 2020-09-22 | Method and device for imaging depth domain of anisotropic medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112305589A true CN112305589A (en) | 2021-02-02 |
Family
ID=74488641
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011000068.1A Pending CN112305589A (en) | 2020-09-22 | 2020-09-22 | Method and device for imaging depth domain of anisotropic medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112305589A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6785612B1 (en) * | 2003-05-29 | 2004-08-31 | Pgs Americas, Inc. | Seismic velocity update for anisotropic depth migration |
| US20100135115A1 (en) * | 2008-12-03 | 2010-06-03 | Chevron U.S.A. Inc. | Multiple anisotropic parameter inversion for a tti earth model |
| US20100312534A1 (en) * | 2008-02-28 | 2010-12-09 | Shiyu Xu | Rock Physics Model For Simulating Seismic Response In Layered Fractured Rocks |
| CN103149588A (en) * | 2013-02-20 | 2013-06-12 | 中国石油天然气股份有限公司 | Method and system for calculating VTI (Velocity Time Integral) anisotropic parameter by utilizing well seismic calibration |
| CN104360385A (en) * | 2014-11-26 | 2015-02-18 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for comprehensively establishing initial depth interval velocity model by combining seismogeology understanding |
| CN109696704A (en) * | 2019-01-28 | 2019-04-30 | 中国海洋石油集团有限公司 | A kind of seismic aeolotropy δ modeling method based on p-wave impedance constraint |
-
2020
- 2020-09-22 CN CN202011000068.1A patent/CN112305589A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6785612B1 (en) * | 2003-05-29 | 2004-08-31 | Pgs Americas, Inc. | Seismic velocity update for anisotropic depth migration |
| US20100312534A1 (en) * | 2008-02-28 | 2010-12-09 | Shiyu Xu | Rock Physics Model For Simulating Seismic Response In Layered Fractured Rocks |
| US20100135115A1 (en) * | 2008-12-03 | 2010-06-03 | Chevron U.S.A. Inc. | Multiple anisotropic parameter inversion for a tti earth model |
| CN103149588A (en) * | 2013-02-20 | 2013-06-12 | 中国石油天然气股份有限公司 | Method and system for calculating VTI (Velocity Time Integral) anisotropic parameter by utilizing well seismic calibration |
| CN104360385A (en) * | 2014-11-26 | 2015-02-18 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for comprehensively establishing initial depth interval velocity model by combining seismogeology understanding |
| CN109696704A (en) * | 2019-01-28 | 2019-04-30 | 中国海洋石油集团有限公司 | A kind of seismic aeolotropy δ modeling method based on p-wave impedance constraint |
Non-Patent Citations (2)
| Title |
|---|
| 杨宗青;李宏伟;欧居刚;常梅斯;林杨;: "TTI介质各向异性参数优化提取方法", 石油地球物理勘探, vol. 55, no. 01, pages 111 - 114 * |
| 韩令贺;胡自多;冯会元;刘威;杨哲;王艳香;: "井震联合网格层析各向异性速度建模研究及应用", 岩性油气藏, no. 04, pages 91 - 97 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106405651B (en) | Full waveform inversion initial velocity model construction method based on logging matching | |
| CN106597533A (en) | Depth domain velocity modeling method for piedmont zone seismic data processing | |
| CN104459782B (en) | Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit | |
| CN107329171A (en) | Depth Domain reservoir seismic inversion method and device | |
| CN113552625B (en) | Multi-scale full waveform inversion method for conventional land-domain seismic data | |
| CN109581501B (en) | Method for modeling depth domain speed of desert area | |
| CN109696704B (en) | Seismic anisotropy delta modeling method based on longitudinal wave impedance constraint | |
| US8531913B2 (en) | Estimating subsurface elastic parameters | |
| CN105626056A (en) | Method and device for determining three-dimensional stratum pore pressure of target area | |
| CN104459798A (en) | Velocity modeling method based on RTM imaging | |
| EP3929630B1 (en) | Computer implemented method for correcting a reservoir model of a reservoir geological formation based on seismic images | |
| CN109725354B (en) | Anisotropic speed modeling method and system | |
| CN111260791B (en) | Method for updating geosteering model | |
| CN112305589A (en) | Method and device for imaging depth domain of anisotropic medium | |
| CN111208558A (en) | Method and device for establishing ultra-deep low-amplitude three-dimensional geological structure | |
| CN110646840A (en) | Angle gather extraction method and system | |
| CN114876454A (en) | Method and system for identifying complex lithologic water flooded layer of highly deviated well | |
| Walker et al. | Stochastic inversion for facies: A case study on the Schiehallion field | |
| US10338248B2 (en) | Method for processing seismic images | |
| CN109581521B (en) | Local chromatographic method and system for TTI anisotropy | |
| CN113933897B (en) | Two-dimensional prestack data-based method and system for predicting gaogite pinch-out line | |
| CN114076982B (en) | Karst paleo-landform restoration method and device based on waveform characteristic difference | |
| CN113267810B (en) | Seismic exploration full-depth velocity modeling method and device | |
| CN110967757A (en) | Real-time pre-stack seismic inversion method and system while drilling | |
| CN113960658B (en) | Logging constraint speed modeling method and device based on geological seismic model |
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 |