CN107561581A - The method that more well models are established based on correlation coefficient process - Google Patents
The method that more well models are established based on correlation coefficient process Download PDFInfo
- Publication number
- CN107561581A CN107561581A CN201610512842.4A CN201610512842A CN107561581A CN 107561581 A CN107561581 A CN 107561581A CN 201610512842 A CN201610512842 A CN 201610512842A CN 107561581 A CN107561581 A CN 107561581A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- well
- mouthful
- msubsup
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
A kind of method that more well models are established based on correlation coefficient process, is comprised the following steps:For every mouthful of well in work area, synthetic seismogram is made, carries out horizon calibration, obtains seismic trace near well;For every mouthful of well in work area, the stratum transformation factor between the earthquake record of the Jing Chu and seismic trace near well is asked for;For every mouthful of well in work area, the coefficient correlation between the earthquake record of the earthquake record of required opening position and opening position of logging well in work area is calculated;For every mouthful of well in work area, the weight of the coefficient correlation of the well is calculated;According to the weight of every mouthful of well in work area, stratum transformation factor, well logging p-wave impedance, well logging S-wave impedance, well logging density, ask for synthesizing p-wave impedance, synthesis S-wave impedance and integral density.This method can establish the initial geological model that both accurate and cans meet actual work area needs in the case of the complex formation of work area.
Description
Technical field
It is more particularly to a kind of to be based on phase relation the present invention relates to the electric digital data processing method of field of geophysical exploration
The method that number method establishes more well models.
Background technology
Geological data is the data medium for including spatial variations relation, and this relation can be that amplitude or other earthquakes are special
Property information, it simultaneously also reflect subsurface reservoir physical property spatial variations.Therefore, on the premise of known geological data, with ground
Shake information is main body, using well information as condition, can obtain the Spatial Variation of geological data, and then establishes one with storage
The initial geological model of sheaf space variation characteristic.
At present, the mode of more well informations is utilized in the method that initial address model is established based on Duo Jing stratum transformation factor
It is to use inverse distance-weighting algorithm, more well informations is introduced by the algorithm, the weight between each well is considered and is allocated.But
There is the limitation of practicality in the algorithm, in actual work area when work area situation complexity, between closest well and required road
Stratigraphic structure it is entirely different, result of calculation mistake can be caused.Therefore the algorithm can take in the case where work area construction is single
Preferable effect is obtained, in the case of complicated work area and is not applied to.
The content of the invention
It is an object of the invention to provide a kind of method that more well models are established based on correlation coefficient process, it enters to existing algorithm
Row improves, and can be applied to a variety of formation conditions.
The present invention uses solution below:
A kind of method that more well models are established based on correlation coefficient process, is comprised the following steps:
Step 1:For every mouthful of well in work area, synthetic seismogram is made, carries out horizon calibration, obtains earthquake by well
Road;
Step 2:For every mouthful of well in work area, the stratum between the earthquake record of the Jing Chu and seismic trace near well is asked for
Transformation factor;
Step 3:For every mouthful of well in work area, the earthquake record of required opening position and well logging opening position in work area are calculated
Coefficient correlation between earthquake record;
Step 4:For every mouthful of well in work area, the weight of the coefficient correlation of the well is calculated;
Step 5:According to the weight of every mouthful of well in work area, stratum transformation factor and well logging p-wave impedance, well logging
S-wave impedance, well logging density, ask for synthesizing p-wave impedance, synthesis S-wave impedance and integral density.
Preferably, synthetic seismogram is made by below equation (1):
s(t)i=w (t)i*r(t)i (1)
Wherein, s (t)iRepresent i-th mouthful of Jing Chu synthetic seismogram, r (t)iRepresent i-th mouthful of Jing Chu reflectance factor sequence
Row, w (t)iRepresent i-th mouthful of Jing Chu seismic wavelet.
Preferably, reflectance factor is calculated by below equation (2):
Wherein, ρ1i、ρ2iThe density that upper and lower two layers of expression interface, v1i、v2iThe speed that upper and lower two layers of expression interface.
Preferably, stratum transformation factor is calculated by below equation (10):
Wherein,Represent i-th mouthful of Jing Chu stratum transformation factor, seis (t)iRepresent i-th mouthful of Jing Chu earthquake record.
Preferably, the coefficient correlation is calculated by below equation (11):
Wherein, x is the earthquake record sample value of required opening position, and y is the earthquake record sample value of well logging opening position, and i is
Well number counting mark, j be earthquake record sampling number counting identify, n be earthquake record total sampling number, xijFor for
J-th of earthquake record sample value of opening position required by i-th mouthful of well, yijRemember for j-th of earthquake for i-th mouthful of borehole logging tool opening position
Record sample value, γiPhase relation between the earthquake record of the well logging opening position of earthquake record and i-th mouthful of well for required opening position
Number.
Preferably, the weight is calculated by below equation (12) using normalization algorithm:
Wherein, εiFor the weight of the coefficient correlation of i-th mouthful of well.
Preferably, synthesis p-wave impedance AI (t) is calculated by below equation (16) to (18)*, synthesis S-wave impedance SI (t)*
With integral density DEN (t)*:
Wherein, AI (t)i、SI(t)i、DEN(t)iI-th mouthful of Jing Chu well logging p-wave impedance, well logging shear wave resistance is represented respectively
Anti-, well logging density.
Compared with prior art, the beneficial effects of the present invention are improve existing building based on inverse distance-weighting
Mould method, required correlation between earthquake record and each borehole-side seismic data is calculated using correlation coefficient process, and with this
To distribute weight, meet the initial of actual work area needs so as to establish both accurate and cans in the case of the complex formation of work area
Geological model.
Brief description of the drawings
Disclosure exemplary embodiment is described in more detail in conjunction with the accompanying drawings, the disclosure it is above-mentioned and other
Purpose, feature and advantage will be apparent.
Fig. 1 shows the flow chart of the method that more well models are established based on correlation coefficient process according to exemplary embodiment.
Embodiment
Preferred embodiment of the present disclosure is more fully described below with reference to accompanying drawings.Although the disclosure is shown in accompanying drawing
Preferred embodiment, however, it is to be appreciated that may be realized in various forms the disclosure without should be limited by embodiments set forth here
System.On the contrary, these embodiments are provided so that the disclosure is more thorough and complete, and can be complete by the scope of the present disclosure
Ground is communicated to those skilled in the art.
Fig. 1 shows the flow chart of the method that more well models are established based on correlation coefficient process according to exemplary embodiment, its
Comprise the following steps:
Step 1:For every mouthful of well in work area, synthetic seismogram is made, carries out horizon calibration, obtains seismic trace near well
The synthetic seismogram s (t) of i-th mouthful of Jing Chu in work areaiAs shown in below equation (1):
s(t)i=w (t)i*r(t)i (1)
Wherein, r (t)iFor i-th mouthful of Jing Chu reflection coefficient sequence, w (t)iFor i-th mouthful of Jing Chu seismic wavelet.
The reflectance factor r (t) at one interfaceiIt can be calculated by the bilevel wave impedance in interface, its expression formula is such as
Shown in below equation (2):
Wherein, ρ1i、ρ2iThe density that upper and lower two layers of expression interface, v1i、v2iThe speed that upper and lower two layers of expression interface, speed
Data can obtain from well-log information.
After generating synthetic seismogram, synthetic seismogram and seismic trace near well are contrasted, to carry out layer position mark
It is fixed.Seismic trace near well is consistent with synthetic seismogram, is also denoted as s (t)i。
By horizon calibration, demarcation and mapping of the well logging layer position to seismic horizon can be achieved so that well logging layer position and earthquake
The matching of layer position.Horizon calibration is carried out by synthetic seismogram and belongs to prior art, will not be repeated here.
Step 2:For every mouthful of well in work area, the stratum between the earthquake record of the Jing Chu and seismic trace near well is asked for
Transformation factor
Relation between earthquake record and seismic trace near well can be expressed by convolution model, as shown in below equation (9):
Wherein, seis (t)iWithI-th mouthful of Jing Chu earthquake record and stratum transformation factor is represented respectively.Earthquake record
seis(t)iIt is by measuring obtained original earthquake data.
Formula (9) gives the relation between earthquake record and seismic trace near well, and earthquake can be obtained by the formula (9)
Stratum transformation factor between record and seismic trace near well, as shown in below equation (10):
Step 3:For every mouthful of well in work area, the earthquake record of required opening position and well logging opening position in work area are calculated
Coefficient correlation between earthquake record
For every mouthful of well in work area, the earthquake record of required opening position and survey in work area are calculated according to below equation (11)
Coefficient correlation between the earthquake record of well opening position:
Wherein, x is the earthquake record sample value of required opening position, and y is the earthquake record sample value of well logging opening position, and i is
Well number counting mark, j be earthquake record sampling number counting identify, n be earthquake record total sampling number, xijFor for
J-th of earthquake record sample value of opening position required by i-th mouthful of well, yijRemember for j-th of earthquake for i-th mouthful of borehole logging tool opening position
Record sample value, γiPhase relation between the earthquake record of the well logging opening position of earthquake record and i-th mouthful of well for required opening position
Number.
Required position can be determined according to the area and other geological conditions in work area by mesh generation.Every earthquake note
Record includes the sample value of all earthquake record sampled points, has both included the earthquake record sample value of required opening position, also comprising well logging
The earthquake record sample value of opening position.
Step 4:For every mouthful of well in work area, the weight of the coefficient correlation of this mouthful of well is calculated
For every mouthful of well, the weight of the coefficient correlation of this mouthful of well is calculated using normalization algorithm, such as below equation (12) institute
Show:
Wherein, εiFor the weight of the coefficient correlation of i-th mouthful of well.
Step 5:According to the weight of every mouthful of well in work area, stratum transformation factor and, well logging p-wave impedance, well logging shear wave
Impedance, well logging density, ask for synthesizing p-wave impedance, synthesis S-wave impedance and integral density
Geological data is the data medium for including spatial variations relation, thus seismic trace near well and well logging p-wave impedance, horizontal stroke
There is association between wave impedance and density data, this association can be expressed by convolution formula, such as below equation (3) extremely
(5) shown in:
AI(t)i=aiw (t)i*s(t)i; (3)
SI(t)i=siw (t)i*s(t)i; (4)
DEN(t)i=denw (t)i*s(t)i; (5)
Wherein, AI (t)i、SI(t)i、DEN(t)i、aiw(t)i、siw(t)i、denw(t)iThe survey of i-th mouthful of well is represented respectively
Well p-wave impedance, well logging S-wave impedance, well logging density, compressional wave matching attribute, shear wave matching attribute, the density matching factor.
Deconvolution conversion is carried out to formula (3) to (5) formula, compressional wave matching attribute aiw (t) can be obtainedi, shear wave matching because
Sub- siw (t)i, density matching factor denw (t)iExpression formula, as shown in below equation (6) to (8):
aiw(t)i=AI (t)i*s(t)i -1 (6)
siw(t)i=SI (t)i*s(t)i -1 (7)
denw(t)i=DEN (t)i*s(t)i -1 (8)
Further, since the spatial variations relation of geological data is certain in same work area, therefore by seismic channel and well
Seismic channel has an identical spatial variations relation, the synthesis p-wave impedance AI (t) at the required earthquake record taken*, synthesis it is horizontal
Wave impedance SI (t)*With integral density DEN (t)*It can be represented respectively by formula (13) to (15):
AI(t)*=aiw (t) * seis (t); (13)
SI(t)*=siw (t) * seis (t); (14)
DEN(t)*=denw (t) * seis (t); (15)
Based on aforementioned formula (3)-(8), formula (10) and formula (13)-(15), can ask for synthesizing p-wave impedance AI
(t)*, synthesis S-wave impedance SI (t)*With integral density DEN (t)*, as shown in below equation (16) to (18):
By formula (16) as can be seen that the synthesis p-wave impedance at the required seismic channel taken be well logging p-wave impedance with
The convolution of stratum transformation factor, and stratum transformation factor then carries out deconvolution acquisition by earthquake record and seismic trace near well.It is similar
Ground, synthesis S-wave impedance are the convolutions of well logging S-wave impedance and stratum transformation factor, and integral density is that well logging density becomes with stratum
Change the convolution of the factor.
Embodiment 1
In embodiment, by synthesize p-wave impedance ask for process exemplified by explanation according to exemplary embodiment based on correlation
The method that Y-factor method Y establishes more well models, synthesis S-wave impedance are similar with the process of asking for of integral density.
First, synthetic seismogram is made, horizon calibration is carried out to log data and geological data, calibrated data exist
Matching is reached on time and depth;
Then, for every mouthful of well in work area, earthquake record seis (t) is asked foriWith seismic trace near well s (t)iBetween ground
Layer transformation factor
Then, for every mouthful of well in work area, the earthquake record of required opening position and well logging opening position in work area are calculated
Coefficient correlation between earthquake record, using normalization algorithm, the coefficient correlation weight of every mouthful of well of distribution;
Finally, according to ask for weight, stratum transformation factor, well logging p-wave impedance, so that it may obtain earthquake record at this
Synthesize p-wave impedance model.
It can ask for synthesizing S-wave impedance and integral density accordingly using same flow.
By asking for stratum transformation factorAnd weight ε is distributed, due toItself it has been by earthquake record and well
What other earthquake record was asked for, obtaining stratum transformation factorAfterwards by itself and well logging p-wave impedance, well logging S-wave impedance, well logging
Density is calculated, and can obtain corresponding each model, i.e., by formula (16), (17), (18), can use simplified side
Formula obtains each model.
Above-mentioned technical proposal is a kind of embodiment of the present invention, for those skilled in the art, in this hair
On the basis of bright principle disclosed, it is easy to make various types of improvement or deformation, it is above-mentioned specific to be not limited solely to the present invention
The description of embodiment, therefore description above is simply preferable, and not restrictive meaning.
Claims (7)
1. a kind of method that more well models are established based on correlation coefficient process, is comprised the following steps:
Step 1:For every mouthful of well in work area, synthetic seismogram is made, carries out horizon calibration, obtains seismic trace near well;
Step 2:For every mouthful of well in work area, the stratum conversion between the earthquake record of the Jing Chu and seismic trace near well is asked for
The factor;
Step 3:For every mouthful of well in work area, the earthquake of the earthquake record of required opening position and well logging opening position in work area is calculated
Coefficient correlation between record;
Step 4:For every mouthful of well in work area, the weight of the coefficient correlation of the well is calculated;
Step 5:According to the weight of every mouthful of well in work area, stratum transformation factor and well logging p-wave impedance, well logging shear wave
Impedance, well logging density, ask for synthesizing p-wave impedance, synthesis S-wave impedance and integral density.
2. the method according to claim 1 that more well models are established based on correlation coefficient process, wherein passing through below equation (1)
Make synthetic seismogram:
s(t)i=w (t)i*r(t)i (1)
Wherein, s (t)iRepresent i-th mouthful of Jing Chu synthetic seismogram, r (t)iRepresent i-th mouthful of Jing Chu reflection coefficient sequence, w
(t)iRepresent i-th mouthful of Jing Chu seismic wavelet.
3. the method according to claim 2 that more well models are established based on correlation coefficient process, wherein passing through below equation (2)
Calculate reflectance factor:
<mrow>
<mi>r</mi>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&rho;</mi>
<mrow>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&rho;</mi>
<mrow>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, ρ1i、ρ2iThe density that upper and lower two layers of expression interface, v1i、v2iThe speed that upper and lower two layers of expression interface.
4. the method according to claim 2 that more well models are established based on correlation coefficient process, wherein passing through below equation
(10) stratum transformation factor is calculated:
<mrow>
<mo>&part;</mo>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>s</mi>
<mi>e</mi>
<mi>i</mi>
<mi>s</mi>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>*</mo>
<mi>s</mi>
<msup>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,Represent i-th mouthful of Jing Chu stratum transformation factor, seis (t)iRepresent i-th mouthful of Jing Chu earthquake record.
5. the method according to claim 4 that more well models are established based on correlation coefficient process, wherein passing through below equation
(11) coefficient correlation is calculated:
<mrow>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>n</mi>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msqrt>
<mrow>
<mi>n</mi>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msubsup>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>*</mo>
<msqrt>
<mrow>
<mi>n</mi>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, x is the earthquake record sample value of required opening position, and y is the earthquake record sample value of well logging opening position, and i is well number
Counting mark, j be earthquake record sampling number counting identify, n be earthquake record total sampling number, xijFor for i-th
J-th of earthquake record sample value of opening position, y required by mouth wellijFor j-th of earthquake record for i-th mouthful of borehole logging tool opening position
Sample value, γiPhase relation between the earthquake record of the well logging opening position of earthquake record and i-th mouthful of well for required opening position
Number.
6. the method according to claim 5 that more well models are established based on correlation coefficient process, wherein using normalization algorithm
The weight is calculated by below equation (12):
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<mrow>
<msub>
<mi>&Sigma;&gamma;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, εiFor the weight of the coefficient correlation of i-th mouthful of well.
7. the method according to claim 6 that more well models are established based on correlation coefficient process, wherein passing through below equation
(16) synthesis p-wave impedance AI (t) is calculated to (18)*, synthesis S-wave impedance SI (t)*With integral density DEN (t)*:
<mrow>
<mi>A</mi>
<mi>I</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
</msup>
<mo>=</mo>
<msub>
<mi>&Sigma;&epsiv;</mi>
<mi>i</mi>
</msub>
<mo>*</mo>
<mi>A</mi>
<mi>I</mi>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>*</mo>
<mo>&part;</mo>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>S</mi>
<mi>I</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
</msup>
<mo>=</mo>
<msub>
<mi>&Sigma;&epsiv;</mi>
<mi>i</mi>
</msub>
<mo>*</mo>
<mi>S</mi>
<mi>I</mi>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>*</mo>
<mo>&part;</mo>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>D</mi>
<mi>E</mi>
<mi>N</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
</msup>
<mo>=</mo>
<msub>
<mi>&Sigma;&epsiv;</mi>
<mi>i</mi>
</msub>
<mo>*</mo>
<mi>D</mi>
<mi>E</mi>
<mi>N</mi>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>*</mo>
<mo>&part;</mo>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, AI (t)i、SI(t)i、DEN(t)iI-th mouthful of Jing Chu well logging p-wave impedance, well logging S-wave impedance, survey is represented respectively
Well density.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610512842.4A CN107561581B (en) | 2016-07-01 | 2016-07-01 | The method for establishing more well models based on correlation coefficient process |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610512842.4A CN107561581B (en) | 2016-07-01 | 2016-07-01 | The method for establishing more well models based on correlation coefficient process |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107561581A true CN107561581A (en) | 2018-01-09 |
CN107561581B CN107561581B (en) | 2019-03-12 |
Family
ID=60969944
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610512842.4A Active CN107561581B (en) | 2016-07-01 | 2016-07-01 | The method for establishing more well models based on correlation coefficient process |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107561581B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112162317A (en) * | 2020-09-28 | 2021-01-01 | 北京中恒利华石油技术研究所 | Method for predicting thin reservoir based on seismic waveform transverse difference |
CN112649843A (en) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | Initial model construction method and system based on layer velocity constraint |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101126815A (en) * | 2006-08-17 | 2008-02-20 | 中国石油天然气股份有限公司 | Method for oil gas detection using lithologic seismic factor and lithologic resistance |
CN103399346A (en) * | 2013-08-15 | 2013-11-20 | 电子科技大学 | Well-to-seismic integration initial wave impedance modeling method |
CN105223617A (en) * | 2015-09-21 | 2016-01-06 | 中国石油化工股份有限公司江汉油田分公司物探研究院 | Seismic properties is utilized to control modeling method and the system of impedance initial value |
-
2016
- 2016-07-01 CN CN201610512842.4A patent/CN107561581B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101126815A (en) * | 2006-08-17 | 2008-02-20 | 中国石油天然气股份有限公司 | Method for oil gas detection using lithologic seismic factor and lithologic resistance |
CN103399346A (en) * | 2013-08-15 | 2013-11-20 | 电子科技大学 | Well-to-seismic integration initial wave impedance modeling method |
CN105223617A (en) * | 2015-09-21 | 2016-01-06 | 中国石油化工股份有限公司江汉油田分公司物探研究院 | Seismic properties is utilized to control modeling method and the system of impedance initial value |
Non-Patent Citations (8)
Title |
---|
ALBERT TARANTOLA: "A strategy for nonlinear elastic inversion of seismic reflection data", 《GEOPHYSICS》 * |
EHSAN JAMALI HONDORI ET AL.: "A new initial model for full waveform inversion without cycle skipping", 《SEG NEW ORLEANS ANNUAL MEETING》 * |
LONG TENG ET AL.: "Probabilistic petrophysical properties estimation integrating viscoelastic AVA inversion with statistical rock physics", 《SEG LAS VEGAS 2012 ANNUAL MEETING》 * |
周单: "地震和相约束的建模方法", 《中国地球物理2009》 * |
夏洪瑞等: "波阻抗反演中的一种建模方法", 《石油物探》 * |
宋海渤: "井—震联合描述储层性质空间变化的方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
张玉芬: "求取波阻抗初值的一种新方法—井约束波阻抗扫描法", 《物探与化探》 * |
董泽斌等: "《地震资料数据处理》", 31 December 1992, 北京:石油工业出版社 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112649843A (en) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | Initial model construction method and system based on layer velocity constraint |
CN112649843B (en) * | 2019-10-10 | 2024-01-23 | 中国石油化工股份有限公司 | Initial model construction method and system based on layer speed constraint |
CN112162317A (en) * | 2020-09-28 | 2021-01-01 | 北京中恒利华石油技术研究所 | Method for predicting thin reservoir based on seismic waveform transverse difference |
Also Published As
Publication number | Publication date |
---|---|
CN107561581B (en) | 2019-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108957549B (en) | Braided river sediment heterogeneous compact sandstone gas reservoir geological modeling method | |
CN113759424B (en) | Karst reservoir filling analysis method and system based on spectral decomposition and machine learning | |
US7493241B2 (en) | 3D velocity modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved for curved-ray prestack time migration and for such migration followed by prestack depth migration | |
US7706981B2 (en) | Method for generating reservoir models utilizing synthetic stratigraphic columns | |
CA2920008C (en) | Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media | |
CN102109616B (en) | Wave impedance inversion method for deposition cycle restriction | |
CN108802812A (en) | A kind of formation lithology inversion method of well shake fusion | |
CN102937720B (en) | Well control improves the method for seismic data resolution | |
CN107329171A (en) | Depth Domain reservoir seismic inversion method and device | |
CN104820243B (en) | Three-dimensional post-stack seismic data amplitude compensation method and device | |
Pendrel | Seismic inversion–still the best tool for reservoir characterization | |
CN106199694A (en) | Synthetic record method based on deep varitron ripple | |
Sheng et al. | Predicting gas content in high-maturity marine shales using artificial intelligence based seismic multiple-attributes analysis: A case study from the lower Silurian Longmaxi Formation, Sichuan Basin, China | |
CN107561581B (en) | The method for establishing more well models based on correlation coefficient process | |
CN107272056B (en) | A method of initial model is established based on the stratum Duo Jing transformation factor | |
CN109283577A (en) | A kind of seismic layer labeling method | |
CN109143335B (en) | Method, system, medium and apparatus for producing synthetic seismic record | |
Chopra et al. | Seismic characterization of a Triassic-Jurassic deep geothermal sandstone reservoir, onshore Denmark, using unsupervised machine learning techniques | |
CN107561580B (en) | The method for establishing initial geological model based on the more wells of multilayer | |
CN106226814B (en) | Utilize converted shear wave seismic data inversion reservoir S-wave impedance and the method for density | |
CN110632660B (en) | Thin sand body characterization method and device based on seismic data body | |
Sang et al. | Prestack simultaneous inversion of P-wave impedance and gas saturation using multi-task residual networks | |
Li et al. | Well Log Prediction Using Deep Sequence Learning | |
CN109444958A (en) | The method of compact reservoir porosity quantitative forecast | |
Nair et al. | Seismic inversion and its applications in reservoir characterization |
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 |