CN112649843B - Initial model construction method and system based on layer speed constraint - Google Patents
Initial model construction method and system based on layer speed constraint Download PDFInfo
- Publication number
- CN112649843B CN112649843B CN201910959464.8A CN201910959464A CN112649843B CN 112649843 B CN112649843 B CN 112649843B CN 201910959464 A CN201910959464 A CN 201910959464A CN 112649843 B CN112649843 B CN 112649843B
- Authority
- CN
- China
- Prior art keywords
- wave impedance
- composite
- density
- synthesized
- well
- 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
- 238000010276 construction Methods 0.000 title claims abstract description 16
- 239000002131 composite material Substances 0.000 claims abstract description 156
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 42
- 230000009466 transformation Effects 0.000 claims abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims description 20
- 230000002194 synthesizing effect Effects 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims 2
- 230000000452 restraining effect Effects 0.000 abstract description 6
- 238000005755 formation reaction Methods 0.000 description 30
- 238000010586 diagram Methods 0.000 description 4
- 238000013507 mapping Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000001010 compromised effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
An initial model construction method and system based on layer speed constraint are disclosed. The method comprises the following steps: acquiring a parawell seismic channel for each stratum of each well in a work area; calculating a formation transformation factor between a seismic record at a formation of the well and a parawell seismic trace; for each well in a work area, calculating a correlation coefficient between the seismic record and the seismic record of each well side seismic channel for each stratum of the well, and further calculating a weight corresponding to the correlation coefficient; calculating the synthesized longitudinal wave impedance, the synthesized transverse wave impedance and the synthesized density; obtaining normalized composite longitudinal wave impedance, normalized composite transverse wave impedance, normalized composite density and normalized layer speed; obtaining constraint synthesized longitudinal wave impedance, constraint synthesized transverse wave impedance and constraint synthesized density; and obtaining final synthesized longitudinal wave impedance, final synthesized transverse wave impedance and final synthesized density, and obtaining an initial model. The invention can obtain a more accurate initial model by restraining the layer speed, and improve the inversion accuracy.
Description
Technical Field
The invention relates to the field of geophysics in petrochemical industry, in particular to an initial model construction method and system based on layer speed constraint.
Background
The importance of the initial model is that since the seismic data is band limited and lacks low frequency components, it is generally impossible to recover information in the full frequency band, so in order to improve the accuracy of inversion, the low frequency components need to be compensated in, and the mode of obtaining the low frequency is the initial model, which is also called a low frequency model. The low frequency model establishes a large structure of the whole inversion result, and if the low frequency model is inaccurate, the accuracy of the whole inversion result is greatly compromised, so that the first step of obtaining the accurate inversion result is to obtain the accurate low frequency model.
At present, in the field of seismic inversion, an initial model is mainly created by interpolation or convolution, and the interpolation method mainly comprises the following steps: the deconvolution method mainly comprises a stratum transformation factor-based method, such as a deconvolution method, a kriging method, a neural network method and the like, and has a common point that extrapolation is performed by starting from well logging data of well points, and the difference is that no seismic data is added in an interpolation method, and the deconvolution method comprises the step of calculating the seismic data. Neither interpolation nor convolution methods have data constraints.
The velocity data can reflect general trends in subsurface formations, particularly layer velocities, but it is not directly applicable in inversion. The current initial modeling method is basically based on interpolation, and starts from a well, and has no constraint, but the purpose of an initial model is to acquire underground general structural information, interpolation is carried out under the structural framework, and interpolation is inaccurate without the framework. Therefore, it is necessary to develop an initial model construction method and system based on layer speed constraint.
The information disclosed in the background section of the invention is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention provides an initial model construction method and system based on layer speed constraint, which can obtain a more accurate initial model by carrying out constraint on layer speed and improve inversion accuracy.
According to one aspect of the invention, an initial model construction method based on layer speed constraint is provided. The method may include: performing horizon calibration on each stratum of each well in the work area and on the synthetic seismic records to obtain a parawell seismic channel; calculating a formation transformation factor between a seismic record at the formation of the well and a parawell seismic trace; for each well in a work area, calculating a correlation coefficient between the seismic record and the seismic record of each well side seismic channel for each stratum of the well, and further calculating a weight corresponding to the correlation coefficient; calculating the synthesized longitudinal wave impedance, the synthesized transverse wave impedance and the synthesized density; obtaining normalized composite longitudinal wave impedance, normalized composite transverse wave impedance, normalized composite density and normalized layer speed; obtaining constraint synthesized longitudinal wave impedance, constraint synthesized transverse wave impedance and constraint synthesized density; and obtaining final synthesized longitudinal wave impedance, final synthesized transverse wave impedance and final synthesized density, and further obtaining an initial model.
Preferably, the side-of-well seismic trace is calculated by equation (1):
s(t) il ＝w(t) il *r(t) il (1)
wherein s (t) il Is the side seismic trace of the ith well and the first stratum, w (t) il The seismic wavelet of the ith well and the first stratum, r (t) il Is the reflection coefficient sequence of the ith well and the first stratum.
Preferably, the formation transformation factor is calculated by equation (2):
wherein,formation transformation factor, seis (t), for the ith well, the first formation il Seismic recording for the ith well, the first formation, s (t) il Is the well side seismic channel of the ith well and the first stratum.
Preferably, the correlation coefficient is calculated by the formula (3):
wherein x is ijl Sample point value, y of jth seismic record of ith well and ith stratum ijl Sample point value and gamma of seismic record of jth side-of-well seismic channel of ith well and jth stratum il And n is the number of sampling points for the seismic record.
Preferably, the weight corresponding to the correlation coefficient is calculated by the formula (4):
wherein ε il Is the weight of the correlation coefficient, gamma il Is a correlation coefficient.
Preferably, the synthetic longitudinal wave impedance is calculated by equation (5):
calculating the synthetic shear wave impedance by equation (6):
calculating the composite density by formula (7):
wherein AI (t) * To synthesize longitudinal wave impedance, SI (t) * To synthesize transverse wave impedance, DEN (t) * Is the synthetic density.
Preferably, the normalized composite longitudinal wave impedance is calculated by equation (8):
calculating the normalized composite transverse wave impedance by equation (9):
calculating the normalized composite density by formula (10):
calculating the normalized layer velocity by formula (11):
wherein,for the kth normalized composite longitudinal wave impedance>For the kth normalized composite transverse wave impedance, a ∈>For the kth normalized composite density, +.>The kth normalized layer speed.
Preferably, the constrained composite longitudinal wave impedance is calculated by equation (12):
calculating the constrained composite shear wave impedance by equation (13):
calculating the constrained composite density by equation (14):
wherein,synthesizing longitudinal wave impedance for kth constraint, < +.>Synthesizing transverse wave impedance for the kth constraint, < +.>For the kth constraint composite density, weight is the adjustment factor.
Preferably, the final synthetic longitudinal wave impedance is calculated by equation (15):
calculating the final synthetic shear wave impedance by equation (16):
calculating the final composite density by equation (17):
wherein,for the kth final composite longitudinal wave impedance, < +.>For the kth final composite transverse wave impedance, < >>For the kth final composite density, AI (t)' is the composite longitudinal wave impedance calculation parameter,/-)>SI (t)' is the synthetic transverse wave impedance calculation parameter, < +.>DEN (t)' is a composite density calculation parameter, +.>
According to another aspect of the present invention, an initial model building system based on layer speed constraint is provided, which is characterized in that the system includes: a memory storing computer executable instructions; a processor executing computer executable instructions in the memory, the processor performing the steps of: performing horizon calibration on each stratum of each well in the work area and on the synthetic seismic records to obtain a parawell seismic channel; calculating a formation transformation factor between a seismic record at the formation of the well and a parawell seismic trace; for each well in a work area, calculating a correlation coefficient between the seismic record and the seismic record of each well side seismic channel for each stratum of the well, and further calculating a weight corresponding to the correlation coefficient; calculating the synthesized longitudinal wave impedance, the synthesized transverse wave impedance and the synthesized density; obtaining normalized composite longitudinal wave impedance, normalized composite transverse wave impedance, normalized composite density and normalized layer speed; obtaining constraint synthesized longitudinal wave impedance, constraint synthesized transverse wave impedance and constraint synthesized density; and obtaining final synthesized longitudinal wave impedance, final synthesized transverse wave impedance and final synthesized density, and further obtaining an initial model.
Preferably, the side-of-well seismic trace is calculated by equation (1):
s(t) il ＝w(t) il *r(t) il (1)
wherein s (t) il Is the side seismic trace of the ith well and the first stratum, w (t) il The seismic wavelet of the ith well and the first stratum, r (t) il Is the reflection coefficient sequence of the ith well and the first stratum.
Preferably, the formation transformation factor is calculated by equation (2):
wherein,formation transformation factor, seis (t), for the ith well, the first formation il Seismic recording for the ith well, the first formation, s (t) il Is the well side seismic channel of the ith well and the first stratum.
Preferably, the correlation coefficient is calculated by the formula (3):
wherein x is ijl Sample point value, y of jth seismic record of ith well and ith stratum ijl Sample point value and gamma of seismic record of jth side-of-well seismic channel of ith well and jth stratum il And n is the number of sampling points for the seismic record.
Preferably, the weight corresponding to the correlation coefficient is calculated by the formula (4):
wherein ε il Is the weight of the correlation coefficient, gamma il Is a correlation coefficient.
Preferably, the synthetic longitudinal wave impedance is calculated by equation (5):
calculating the synthetic shear wave impedance by equation (6):
calculating the composite density by formula (7):
wherein AI (t) * To synthesize longitudinal wave impedance, SI (t) * To synthesize transverse wave impedance, DEN (t) * Is the synthetic density.
Preferably, the normalized composite longitudinal wave impedance is calculated by equation (8):
calculating the normalized composite transverse wave impedance by equation (9):
calculating the normalized composite density by formula (10):
calculating the normalized layer velocity by formula (11):
wherein,for the kth normalized composite longitudinal wave impedance>For the kth normalized composite transverse wave impedance, a ∈>For the kth normalized composite density, +.>The kth normalized layer speed.
Preferably, the constrained composite longitudinal wave impedance is calculated by equation (12):
calculating the constrained composite shear wave impedance by equation (13):
calculating the constrained composite density by equation (14):
wherein,synthesizing longitudinal wave impedance for kth constraint, < +.>Synthesizing a tie for the kth constraintWave impedance, < >>For the kth constraint composite density, weight is the adjustment factor.
Preferably, the final synthetic longitudinal wave impedance is calculated by equation (15):
calculating the final synthetic shear wave impedance by equation (16):
calculating the final composite density by equation (17):
wherein,for the kth final composite longitudinal wave impedance, < +.>For the kth final composite transverse wave impedance, < >>For the kth final composite density, AI (t)' is the composite longitudinal wave impedance calculation parameter,/-)>SI (t)' is the synthetic transverse wave impedance calculation parameter, < +.>DEN (t)' is a composite density calculation parameter, +.>
The method and apparatus of the present invention have other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the present invention.
Drawings
The foregoing and other objects, features and advantages of the invention will be apparent from the following more particular descriptions of exemplary embodiments of the invention as illustrated in the accompanying drawings wherein like reference numbers generally represent like parts throughout the exemplary embodiments of the invention.
Fig. 1 shows a flow chart of the steps of an initial model building method based on layer speed constraints according to the present invention.
Fig. 2 shows a schematic diagram of layer speed data according to an embodiment of the invention.
FIG. 3 shows a schematic diagram of an initial model according to one embodiment of the invention.
Detailed Description
The invention will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present invention are illustrated in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
Fig. 1 shows a flow chart of the steps of an initial model building method based on layer speed constraints according to the present invention.
In this embodiment, the initial model construction method based on the layer speed constraint according to the present invention may include: step 101: performing horizon calibration on each stratum of each well in the work area and on the synthetic seismic records to obtain a parawell seismic channel; step 102: calculating a formation transformation factor between a seismic record at a formation of the well and a parawell seismic trace; step 103: for each well in a work area, calculating a correlation coefficient between the seismic record and the seismic record of each well side seismic channel for each stratum of the well, and further calculating a weight corresponding to the correlation coefficient; step 104: calculating the synthesized longitudinal wave impedance, the synthesized transverse wave impedance and the synthesized density; step 105: obtaining normalized composite longitudinal wave impedance, normalized composite transverse wave impedance, normalized composite density and normalized layer speed; step 106: obtaining constraint synthesized longitudinal wave impedance, constraint synthesized transverse wave impedance and constraint synthesized density; step 107: and obtaining final synthesized longitudinal wave impedance, final synthesized transverse wave impedance and final synthesized density, and further obtaining an initial model.
In one example, a parawell seismic trace is calculated by equation (1):
s(t) il ＝w(t) il *r(t) il (1)
wherein s (t) il Is the side seismic trace of the ith well and the first stratum, w (t) il The seismic wavelet of the ith well and the first stratum, r (t) il Is the reflection coefficient sequence of the ith well and the first stratum.
In one example, the formation transformation factor is calculated by equation (2):
wherein,formation transformation factor, seis (t), for the ith well, the first formation il Seismic recording for the ith well, the first formation, s (t) il Is the well side seismic channel of the ith well and the first stratum.
In one example, the correlation coefficient is calculated by equation (3):
wherein x is ijl Sample point value, y of jth seismic record of ith well and ith stratum ijl Sample point value and gamma of seismic record of jth side-of-well seismic channel of ith well and jth stratum il And n is the number of sampling points for the seismic record.
In one example, the weight corresponding to the correlation coefficient is calculated by equation (4):
wherein ε il Is the weight of the correlation coefficient, gamma il Is a correlation coefficient.
In one example, the composite longitudinal wave impedance is calculated by equation (5):
the synthetic transverse wave impedance is calculated by equation (6):
the composite density is calculated by equation (7):
wherein AI (t) * To synthesize longitudinal wave impedance, SI (t) * To synthesize transverse wave impedance, DEN (t) * Is the synthetic density.
In one example, the normalized composite longitudinal wave impedance is calculated by equation (8):
calculating normalized composite transverse wave impedance by equation (9):
the normalized composite density is calculated by equation (10):
calculating a normalized layer speed by formula (11):
wherein,for the kth normalized composite longitudinal wave impedance>For the kth normalized composite transverse wave impedance, a ∈>For the kth normalized composite density, +.>The kth normalized layer speed.
In one example, the constrained composite longitudinal wave impedance is calculated by equation (12):
the constrained composite transverse wave impedance is calculated by equation (13):
the constrained composite density is calculated by equation (14):
wherein,synthesizing longitudinal wave impedance for kth constraint, < +.>Synthesizing transverse wave impedance for the kth constraint, < +.>For the kth constraint composite density, weight is the adjustment factor.
In one example, the final composite longitudinal wave impedance is calculated by equation (15):
the final composite transverse wave impedance is calculated by equation (16):
the final composite density is calculated by equation (17):
wherein,for the kth final composite longitudinal wave impedance, < +.>For the kth final composite transverse wave impedance, < >>For the kth final composite density, AI (t)' is the composite longitudinal wave impedance calculation parameter,/-)>SI (t)' is the synthetic transverse wave impedance calculation parameter, < +.>DEN (t)' is a composite density calculation parameter, +.>
Specifically, the initial model construction method based on the layer speed constraint according to the present invention may include:
for each stratum of each well in a work area, performing horizon calibration for synthetic seismic records, and obtaining a side-of-well seismic trace through a formula (1), wherein the reflection coefficient r (t) of one interface is obtained il The wave impedance of the upper layer and the lower layer is obtained, and the expression is as follows:
wherein r (t) il Is the reflection coefficient ρ il 、ρ i(l+1) For the density of the upper and lower layers, v il 、v i(l+1) For the velocities of the upper and lower layers, data is obtained from the log.
And (5) obtaining the relation between the side-well seismic channel and the longitudinal wave impedance, the transverse wave impedance and the density of the well logging along the layer.
The seismic data itself is a data carrier containing a spatially varying relationship, so that there is a correlation between the parawell seismic trace and the log longitudinal wave impedance, transverse wave impedance and density data, which relationship can be expressed by a convolution formula:
AI(t) il ＝aiw(t) il *s(t) il (19)
SI(t) il ＝siw(t) il *s(t) il (20)
DEN(t) il ＝denw(t) il *s(t) il (21)
wherein at the ith well at the first layer, AI (t) il For logging longitudinal wave impedance, SI (t) il For logging shear wave impedance, DEN (t) il Aiw (t) is the logging density il For longitudinal wave matching factor, siw (t) il As transversal wave matching factor, new (t) il Is a density matching factor. Deconvolution of formulas (19) - (21) can be obtained:
aiw(t) il ＝AI(t) il *s(t) il -1 (22)
siw(t) il ＝SI(t) il *s(t) il -1 (23)
denw(t) il ＝DEN(t) il *s(t) il -1 (24)
the relation between the seismic record and the well side seismic channel is obtained along the layer, and the relation can be obtained through a convolution model:
transforming equation (25) may yield equation (2), by which a formation transformation factor between the seismic record at the formation of the well and the parawell seismic trace is calculated.
For each well in the work area, calculating a correlation coefficient between the seismic record and the seismic record of each well side seismic channel according to a formula (3) for each stratum of the well, and further calculating a weight corresponding to the correlation coefficient according to a formula (4).
Because the spatial variation relationship of the seismic data is fixed in the same work area, the seismic trace and the side seismic trace have the same spatial variation relationship, and can be expressed by the formula:
wherein AI (t) * 、SI(t) * 、DEN(t) * The composite compressional impedance, composite shear impedance, and composite density of seis (t) at the desired seismic record.
Substituting the formula (22) and the formula (2) into the formula (26) to obtain the synthesized longitudinal wave impedance as the formula (5), and similarly obtaining the formula (6) and the formula (7), wherein the synthesized longitudinal wave impedance, the synthesized transverse wave impedance and the synthesized density can be calculated through the formulas (5) - (7) respectively.
Since the layer velocity is the velocity of seismic wave propagation in the layered formation, it directly reflects the lithology of the formation and can be used to divide the formation, it itself contains low frequency information and can be used as constraint data to participate in calculations. Normalizing the composite longitudinal wave impedance, the composite transverse wave impedance, the composite density and the layer speed to be between 0 and 1, and calculating normalized composite longitudinal wave impedance, normalized composite transverse wave impedance, normalized composite density and normalized layer speed through formulas (8) - (11).
Setting an adjusting factor weight to be a pure decimal between 0 and 1, setting different values to adjust specific gravity according to actual needs during calculation, restraining normalized synthesized longitudinal wave impedance, normalized synthesized transverse wave impedance and normalized synthesized density according to the adjusting factor and normalized layer speed, and calculating restrained synthesized longitudinal wave impedance, restrained synthesized transverse wave impedance and restrained synthesized density through formulas (12) - (14).
And (3) mapping the constrained composite longitudinal wave impedance, the constrained composite transverse wave impedance and the constrained composite density to value ranges respectively, and calculating the final composite longitudinal wave impedance, the final composite transverse wave impedance and the final composite density through formulas (15) - (17).
The method can obtain a more accurate initial model by restraining the layer speed, and improves the inversion accuracy.
Application example
In order to facilitate understanding of the solution and the effects of the embodiments of the present invention, a specific application example is given below. It will be understood by those of ordinary skill in the art that the examples are for ease of understanding only and that any particular details thereof are not intended to limit the present invention in any way.
The initial model construction method based on the layer speed constraint according to the present invention may include:
for each stratum of each well in a work area, performing horizon calibration for synthetic seismic records, and obtaining a side-of-well seismic trace through a formula (1), wherein the reflection coefficient r (t) of one interface is obtained il The wave impedance of the upper layer and the lower layer is obtained, and the expression is as follows:
wherein r (t) il Is the reflection coefficient ρ il 、ρ i(l+1) For the density of the upper and lower layers, v il 、v i(l+1) For the velocities of the upper and lower layers, data is obtained from the log.
And (4) solving the relation between the side-well seismic channel and the longitudinal wave impedance, the transverse wave impedance and the density of the well logging along the layer, wherein the longitudinal wave matching factor, the transverse wave matching factor and the density matching factor can be obtained by adopting the formulas (22) - (24) respectively.
And (3) obtaining the relation between the seismic records and the well side seismic traces along the layers, and calculating stratum transformation factors between the seismic records at the stratum of the well and the well side seismic traces through a formula (2).
For each well in the work area, calculating a correlation coefficient between the seismic record and the seismic record of each well side seismic channel according to a formula (3) for each stratum of the well, and further calculating a weight corresponding to the correlation coefficient according to a formula (4). And then calculating the synthesized longitudinal wave impedance, the synthesized transverse wave impedance and the synthesized density through formulas (5) - (7), respectively.
Fig. 2 shows a schematic diagram of layer speed data according to an embodiment of the invention.
Normalizing the composite longitudinal wave impedance, the composite transverse wave impedance, the composite density and the layer velocity data as shown in fig. 2 to between 0 and 1, and calculating normalized composite longitudinal wave impedance, normalized composite transverse wave impedance, normalized composite density and normalized layer velocity by formulas (8) - (11).
Setting the adjustment factor weight to be a pure decimal between 0 and 1, restraining the normalized composite longitudinal wave impedance, the normalized composite transverse wave impedance and the normalized composite density according to the adjustment factor and the normalized layer speed, and calculating the restrained composite longitudinal wave impedance, the restrained composite transverse wave impedance and the restrained composite density through formulas (12) - (14).
FIG. 3 shows a schematic diagram of an initial model according to one embodiment of the invention.
And mapping the constrained composite longitudinal wave impedance, the constrained composite transverse wave impedance and the constrained composite density to value ranges respectively, and calculating the final composite longitudinal wave impedance, the final composite transverse wave impedance and the final composite density through formulas (15) - (17), so as to obtain an initial model, as shown in fig. 3.
The method uses the layer speed as constraint data to create an initial model, so that the trend of the underground structure can be added into the initial model, and the method can add the trend into three elastic parameters to perform prestack inversion.
In summary, the invention can obtain a more accurate initial model by restraining the layer speed, and improve the inversion accuracy.
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention has been given for the purpose of illustrating the benefits of embodiments of the invention only and is not intended to limit embodiments of the invention to any examples given.
According to an embodiment of the present invention, there is provided an initial model building system based on layer speed constraint, characterized in that the system includes: a memory storing computer executable instructions; a processor executing computer executable instructions in the memory, the processor performing the steps of: performing horizon calibration on each stratum of each well in the work area and on the synthetic seismic records to obtain a parawell seismic channel; calculating a formation transformation factor between a seismic record at a formation of the well and a parawell seismic trace; for each well in a work area, calculating a correlation coefficient between the seismic record and the seismic record of each well side seismic channel for each stratum of the well, and further calculating a weight corresponding to the correlation coefficient; calculating the synthesized longitudinal wave impedance, the synthesized transverse wave impedance and the synthesized density; obtaining normalized composite longitudinal wave impedance, normalized composite transverse wave impedance, normalized composite density and normalized layer speed; obtaining constraint synthesized longitudinal wave impedance, constraint synthesized transverse wave impedance and constraint synthesized density; and obtaining final synthesized longitudinal wave impedance, final synthesized transverse wave impedance and final synthesized density, and further obtaining an initial model.
In one example, a parawell seismic trace is calculated by equation (1):
s(t) il ＝w(t) il *r(t) il (1)
wherein s (t) il Is the side seismic trace of the ith well and the first stratum, w (t) il The seismic wavelet of the ith well and the first stratum, r (t) il Is the reflection coefficient sequence of the ith well and the first stratum.
In one example, the formation transformation factor is calculated by equation (2):
wherein,formation transformation factor, seis (t), for the ith well, the first formation il Seismic recording for the ith well, the first formation, s (t) il Is the well side seismic channel of the ith well and the first stratum.
In one example, the correlation coefficient is calculated by equation (3):
wherein x is ijl Sample point value, y of jth seismic record of ith well and ith stratum ijl Sample point value and gamma of seismic record of jth side-of-well seismic channel of ith well and jth stratum il And n is the number of sampling points for the seismic record.
In one example, the weight corresponding to the correlation coefficient is calculated by equation (4):
wherein ε il Is the weight of the correlation coefficient, gamma il Is a correlation coefficient.
In one example, the composite longitudinal wave impedance is calculated by equation (5):
the synthetic transverse wave impedance is calculated by equation (6):
the composite density is calculated by equation (7):
wherein AI (t) * To synthesize longitudinal wave impedance, SI (t) * To synthesize transverse wave impedance, DEN (t) * Is the synthetic density.
In one example, the normalized composite longitudinal wave impedance is calculated by equation (8):
calculating normalized composite transverse wave impedance by equation (9):
the normalized composite density is calculated by equation (10):
calculating a normalized layer speed by formula (11):
wherein,for the kth normalized composite longitudinal wave impedance>For the kth normalized composite transverse wave impedance, a ∈>For the kth normalized composite density, +.>The kth normalized layer speed.
In one example, the constrained composite longitudinal wave impedance is calculated by equation (12):
the constrained composite transverse wave impedance is calculated by equation (13):
the constrained composite density is calculated by equation (14):
wherein,synthesizing longitudinal wave impedance for kth constraint, < +.>Synthesizing transverse wave impedance for the kth constraint, < +.>For the kth constraint composite density, weight is the adjustment factor.
In one example, the final composite longitudinal wave impedance is calculated by equation (15):
the final composite transverse wave impedance is calculated by equation (16):
the final composite density is calculated by equation (17):
wherein,for the kth final composite longitudinal wave impedance, < +.>For the kth final composite transverse wave impedance, < >>For the kth final composite density, AI (t)' is the composite longitudinal wave impedance calculation parameter,/-)>SI (t)' is the synthetic transverse wave impedance calculation parameter, < +.>DEN (t)' is a composite density calculation parameter, +.>
The system can obtain a more accurate initial model by restraining the layer speed, and improves the inversion accuracy.
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention has been given for the purpose of illustrating the benefits of embodiments of the invention only and is not intended to limit embodiments of the invention to any examples given.
The foregoing description of embodiments of the invention has been presented for purposes of illustration and description, and is not intended to be exhaustive or limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the various embodiments described.
Claims (8)
1. An initial model construction method based on layer speed constraint is characterized by comprising the following steps:
performing horizon calibration on each stratum of each well in the work area and on the synthetic seismic records to obtain a parawell seismic channel;
calculating a formation transformation factor between a seismic record at the formation of the well and a parawell seismic trace;
for each well in a work area, calculating a correlation coefficient between the seismic record and the seismic record of each well side seismic channel for each stratum of the well, and further calculating a weight corresponding to the correlation coefficient;
calculating the synthesized longitudinal wave impedance, the synthesized transverse wave impedance and the synthesized density;
obtaining normalized composite longitudinal wave impedance, normalized composite transverse wave impedance, normalized composite density and normalized layer speed;
obtaining constraint synthesized longitudinal wave impedance, constraint synthesized transverse wave impedance and constraint synthesized density;
obtaining final synthesized longitudinal wave impedance, final synthesized transverse wave impedance and final synthesized density, and further obtaining an initial model;
wherein the constrained composite longitudinal wave impedance is calculated by equation (12):
calculating the constrained composite shear wave impedance by equation (13):
calculating the constrained composite density by equation (14):
wherein,synthesizing longitudinal wave impedance for kth constraint, < +.>Synthesizing transverse wave impedance for the kth constraint, < +.>For the kth constraint synthesis density, weight is the regulator, ++>For the kth normalized composite longitudinal wave impedance>For the kth normalized composite transverse wave impedance, a ∈>For the kth normalized composite density, +.>Normalized layer speed for kth;
wherein the final synthetic longitudinal wave impedance is calculated by equation (15):
calculating the final synthetic shear wave impedance by equation (16):
calculating the final composite density by equation (17):
wherein,for the kth final composite longitudinal wave impedance, < +.>For the kth final composite transverse wave impedance, < >>For the kth final composite density, AI (t)' is the composite longitudinal wave impedance calculation parameter,/-)>SI (t)' is the synthetic transverse wave impedance calculation parameter, < +.>DEN (t)' is a composite density calculation parameter, +.>
2. The initial model construction method based on layer velocity constraints of claim 1, wherein the parawell seismic trace is calculated by equation (1):
s(t) il ＝w(t) il *r(t) il (1)
wherein s (t) il Is the side seismic trace of the ith well and the first stratum, w (t) il The seismic wavelet of the ith well and the first stratum, r (t) il Is the reflection coefficient sequence of the ith well and the first stratum.
3. The initial model construction method based on layer velocity constraints according to claim 1, wherein the formation transformation factor is calculated by formula (2):
wherein,formation transformation factor, seis (t), for the ith well, the first formation il Seismic recording for the ith well, the first formation, s (t) il Is the well side seismic channel of the ith well and the first stratum.
4. The initial model construction method based on layer speed constraint according to claim 3, wherein the correlation coefficient is calculated by formula (3):
wherein x is ijl Sample point value, y of jth seismic record of ith well and ith stratum ijl Sample point value and gamma of seismic record of jth side-of-well seismic channel of ith well and jth stratum il And n is the number of sampling points of the seismic record.
5. The initial model construction method based on layer speed constraint according to claim 4, wherein the weight corresponding to the correlation coefficient is calculated by formula (4):
wherein ε il Is the weight of the correlation coefficient, gamma il Is a correlation coefficient.
6. The initial model construction method based on layer velocity constraints according to claim 5, wherein the synthetic longitudinal wave impedance is calculated by formula (5):
calculating the synthetic shear wave impedance by equation (6):
calculating the composite density by formula (7):
wherein AI (t) * To synthesize longitudinal wave impedance, SI (t) * To synthesize transverse wave impedance, DEN (t) * To synthesize density, at the ith well at the first layer, AI (t) il For logging longitudinal wave impedance, SI (t) il For logging shear wave impedance, DEN (t) il Is the log density.
7. The initial model building method based on layer velocity constraints according to claim 1, wherein the normalized composite longitudinal wave impedance is calculated by formula (8):
calculating the normalized composite transverse wave impedance by equation (9):
calculating the normalized composite density by formula (10):
calculating the normalized layer velocity by formula (11):
8. an initial model building system based on layer speed constraints, the system comprising:
a memory storing computer executable instructions;
a processor executing computer executable instructions in the memory, the processor performing the steps of:
performing horizon calibration on each stratum of each well in the work area and on the synthetic seismic records to obtain a parawell seismic channel;
calculating a formation transformation factor between a seismic record at the formation of the well and a parawell seismic trace;
for each well in a work area, calculating a correlation coefficient between the seismic record and the seismic record of each well side seismic channel for each stratum of the well, and further calculating a weight corresponding to the correlation coefficient;
calculating the synthesized longitudinal wave impedance, the synthesized transverse wave impedance and the synthesized density;
obtaining normalized composite longitudinal wave impedance, normalized composite transverse wave impedance, normalized composite density and normalized layer speed;
obtaining constraint synthesized longitudinal wave impedance, constraint synthesized transverse wave impedance and constraint synthesized density;
obtaining final synthesized longitudinal wave impedance, final synthesized transverse wave impedance and final synthesized density, and further obtaining an initial model;
wherein the constrained composite longitudinal wave impedance is calculated by equation (12):
calculating the constrained composite shear wave impedance by equation (13):
calculating the constrained composite density by equation (14):
wherein,synthesizing longitudinal wave impedance for kth constraint, < +.>The transverse wave impedance is synthesized for the kth constraint,for the kth constraint synthesis density, weight is the regulator, ++>For the kth normalized composite longitudinal wave impedance>For the kth normalized composite transverse wave impedance, a ∈>For the kth normalized composite density,normalized layer speed for kth;
wherein the final synthetic longitudinal wave impedance is calculated by equation (15):
calculating the final synthetic shear wave impedance by equation (16):
calculating the final composite density by equation (17):
wherein,for the kth final composite longitudinal wave impedance, < +.>For the kth final composite transverse wave impedance, < >>For the kth final composite density, AI (t)' is the composite longitudinal wave impedance calculation parameter,/-)>SI (t)' is the synthetic transverse wave impedance calculation parameter, < +.>DEN (t)' is a composite density calculation parameter, +.>
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910959464.8A CN112649843B (en) | 2019-10-10 | 2019-10-10 | Initial model construction method and system based on layer speed constraint |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910959464.8A CN112649843B (en) | 2019-10-10 | 2019-10-10 | Initial model construction method and system based on layer speed constraint |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112649843A CN112649843A (en) | 2021-04-13 |
| CN112649843B true CN112649843B (en) | 2024-01-23 |
Family
ID=75343449
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910959464.8A Active CN112649843B (en) | 2019-10-10 | 2019-10-10 | Initial model construction method and system based on layer speed constraint |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112649843B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101231346A (en) * | 2008-02-26 | 2008-07-30 | 中国矿业大学(北京) | Method for estimating coal, rock mass physical mechanics parameter through seismic wave velocity |
| CN103487835A (en) * | 2012-06-12 | 2014-01-01 | 中国石油化工股份有限公司 | Multi-resolution wave impedance inversion method based on model constraints |
| WO2016135505A1 (en) * | 2015-02-26 | 2016-09-01 | Ireson Richard | A method of producing and utilising high resolution impedance logs derived from vsp data for use in assessing an oilfield subterranean formation |
| CN107272056A (en) * | 2016-04-06 | 2017-10-20 | 中国石油化工股份有限公司 | A kind of method that initial model is set up based on Duo Jing stratum transformation factor |
| CN107561580A (en) * | 2016-07-01 | 2018-01-09 | 中国石油化工股份有限公司 | The method for establishing initial geological model based on the more wells of multilayer |
| CN107561581A (en) * | 2016-07-01 | 2018-01-09 | 中国石油化工股份有限公司 | The method that more well models are established based on correlation coefficient process |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2895091B1 (en) * | 2005-12-21 | 2008-02-22 | Inst Francais Du Petrole | METHOD FOR UPDATING A GEOLOGICAL MODEL WITH SEISMIC DATA |
| US11487036B2 (en) * | 2017-01-12 | 2022-11-01 | Cgg Services Sas | Reflection full waveform inversion methods with density and velocity models updated separately |
| US10739480B2 (en) * | 2017-03-24 | 2020-08-11 | Exxonmobil Upstream Research Company | Full wavefield inversion with reflected seismic data starting from a poor velocity model |
-
2019
- 2019-10-10 CN CN201910959464.8A patent/CN112649843B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101231346A (en) * | 2008-02-26 | 2008-07-30 | 中国矿业大学(北京) | Method for estimating coal, rock mass physical mechanics parameter through seismic wave velocity |
| CN103487835A (en) * | 2012-06-12 | 2014-01-01 | 中国石油化工股份有限公司 | Multi-resolution wave impedance inversion method based on model constraints |
| WO2016135505A1 (en) * | 2015-02-26 | 2016-09-01 | Ireson Richard | A method of producing and utilising high resolution impedance logs derived from vsp data for use in assessing an oilfield subterranean formation |
| CN107272056A (en) * | 2016-04-06 | 2017-10-20 | 中国石油化工股份有限公司 | A kind of method that initial model is set up based on Duo Jing stratum transformation factor |
| CN107561580A (en) * | 2016-07-01 | 2018-01-09 | 中国石油化工股份有限公司 | The method for establishing initial geological model based on the more wells of multilayer |
| CN107561581A (en) * | 2016-07-01 | 2018-01-09 | 中国石油化工股份有限公司 | The method that more well models are established based on correlation coefficient process |
Non-Patent Citations (2)
| Title |
|---|
| "Waveform-Preserving Processing Flow of Multichannel Seismic Reflection Data for Adjoint-State Full-Waveform Inversion of Ocean Thermohaline Structure";Daniel Dagnino;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;第56卷(第3期);全文 * |
| "基于地震谱反演的地层速度估算方法及应用";刘洁;《石油地球物理勘探》;第51卷(第5期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112649843A (en) | 2021-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102129084B (en) | Method and device for acquiring seismic thin reservoir speed through well control | |
| US7840394B2 (en) | Method for generating a 3D earth model | |
| EP3028071B1 (en) | Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media | |
| US11333783B2 (en) | Integrated method for estimation of seismic wavelets and synthesis of seismic records in depth domain | |
| US7525874B2 (en) | Method for converting seismic data from the time domain to the depth domain | |
| US20140350861A1 (en) | Multi-parameter inversion through offset dependent elastic fwi | |
| US20020013661A1 (en) | Method of estimating elastic and compositional parameters from seismic and echo-acoustic data | |
| US20230296799A1 (en) | Work flow based acoustic processing system and method | |
| US11493658B2 (en) | Computer-implemented method and system employing nonlinear direct prestack seismic inversion for poisson impedance | |
| US20090120634A1 (en) | Method for Tomographic Inversion by Matrix Transformation | |
| UA76697C2 (en) | Method for predicting characteristics of rock or fluid by seismic exploration data | |
| Landa et al. | Interpretation of velocity estimates from coherency inversion | |
| EP2073041A1 (en) | Method to estimate a seismic ray parameter for a seismogram | |
| CN113031068A (en) | Reflection coefficient accurate base tracking prestack seismic inversion method | |
| US7382683B1 (en) | Computing an absorption parameter for a mode-converted seismic wave | |
| CN112649843B (en) | Initial model construction method and system based on layer speed constraint | |
| EP3203269B1 (en) | Multi-stack (broadband) seismic wavelet estimation method | |
| CN107561581B (en) | The method for establishing more well models based on correlation coefficient process | |
| US6778907B1 (en) | Method for estimation of propagation paths for seismic signals | |
| CN109143335B (en) | Method, system, medium and apparatus for producing synthetic seismic record | |
| CN112147700A (en) | Low-frequency model construction method and system for speed abnormal area | |
| Gao et al. | Multicomponent seismic data registration by nonlinear optimization | |
| CN107561580B (en) | The method for establishing initial geological model based on the more wells of multilayer | |
| WO2014065785A1 (en) | Data double-searching apparatus, methods, and systems | |
| CN112147681A (en) | Pre-stack inversion method and system based on gamma _ Zoeppritz equation |
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 |