CN106094051B - A kind of granitic formation Division and contrast method - Google Patents
A kind of granitic formation Division and contrast method Download PDFInfo
- Publication number
- CN106094051B CN106094051B CN201610403178.XA CN201610403178A CN106094051B CN 106094051 B CN106094051 B CN 106094051B CN 201610403178 A CN201610403178 A CN 201610403178A CN 106094051 B CN106094051 B CN 106094051B
- Authority
- CN
- China
- Prior art keywords
- well
- curve
- mrow
- point
- spectrum
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a kind of granitic formation Division and contrast method, it comprises the following steps：1) borehole log data prepares；2) log maximum entropy spectrum analysis；3) data difference curve is obtained；4) it is spectrum curve by data difference Curve transform；5) processing is filtered to the spectrum curve of reconstruct；6) abnormity point is determined；7) abnormity point is numbered；8) classification and correlation is carried out：Each well key similar differences point on key section is sequentially connected, that is, completes granitic formation division and contrast.The present invention is to carry out com-parison and analysis to the abnormity point in each well, compared with traditional similar stratum control methods, for stratum small to this lithology difference of granite, vertical variations are small, cyclicity does not protrude, there is stronger operability, more conventional Strata Comparison method has more practicality.The present invention is as a result of maximum entropy spectrum analysis method, therefore Strata Comparison quantitative extent is strengthened, and operability is stronger, and the efficiency of Strata Comparison greatly reinforces.
Description
Technical field
The present invention relates to a kind of granitic formation Division and contrast method, belong to oil field prospecting field.
Background technology
Classification and correlation is oil field prospecting and a link essential in development, and follow-up geology is built
Mould, straticulate exploration division and well net deployment etc. are all carried out on the basis of accurate Stratigraphic Division And Correlation.It is currently based on
It is a variety of that conventional crumb rock and carbonate strata have explored such as seismic stratigraphy, sequence stratigraphy, chronostratigraphy etc.
Stratigraphic Division And Correlation method, and form corresponding theoretical system and technological means.But granitic formation is divided but to rarely have and related to
And it is primarily present two aspect reasons, be on the one hand existing theoretical system by being then based on time domain, it is emphasised that waiting
Shi Linian, it is closely related with cycle of sedimentation concept, and granitic formation is most of more single without cycle of sedimentation concept, lithology,
Stratum inside story Seismic reflection character is generally unclear, the vertical electrical property feature of individual well is mostly more uniform, therefore is not suitable for magma and invades
Enter and subsequent construction, slacking formed reservoir；On the other hand it is that the discovery of granites types oil field is less, therefore grinds accordingly
Study carefully that technological means is relatively deficient, Stratigraphic Division And Correlation thinking during tradition etc. is not suitable for such reservoir formation Division and contrast work
Make.
The content of the invention
In view of the above-mentioned problems, it is an object of the invention to provide a kind of granitic formation Division and contrast based on abnormal point analysis
Method.
To achieve the above object, the present invention takes following technical scheme：A kind of granitic formation Division and contrast method, it is wrapped
Include following steps：
1) borehole log data prepares：By the negative of the start-stop position of the log of each well in oil field to be measured and stitching position
Exceptional value is changed into 0；
2) log maximum entropy spectrum analysis：At borehole log data by each well of maximum entropy spectrum analysis function pair
Reason, obtain the well logging spectrum curve of each well；
3) data difference curve is obtained：The maximum entropy of each depth point of the well logging spectrum curve of each well is calculated, and
According to the real data of log corresponding to each depth point, the data difference curve of each well is obtained；
4) it is spectrum curve by data difference Curve transform：The data difference curve negotiating iterative Simpson of each well is integrated
Formula carries out mathematic integral conversion, reconstructs the spectrum curve of each well；
5) processing is filtered to the spectrum curve of reconstruct：According to every mouthful of well well logging and sidewall coring lithology in oil field to be measured
Statistical result obtains the upper limit, intermediate value and the lower limit of weathered stratum thickness, and is closed according to the positive correlation of EWL and formation thickness
System obtains the EWL of the upper limit of weathered stratum thickness, intermediate value and lower limit, and the scope according to EWL chooses filtering parameter,
So as to be filtered processing to the spectrum curve of reconstruct, three filter curves of the upper limit of each well, intermediate value and lower limit are obtained；
6) abnormity point is determined：According to three filter curves of each well obtained after filtering process, choose in three curves
In the presence of large change amplitude point as discrepancy；
7) abnormity point is numbered：The abnormity point of each well to being determined in step 6) carries out vertical numbering from top to bottom；
8) classification and correlation is carried out：Each well key similar differences point on key section is sequentially connected, that is, completes flower hilllock
Rock classification and correlation.
In the step 2), function and maximum entropy SX(ω) and correlation function rx(m) relation between is made by imperial order by wiener-pungent
Reason determination, i.e.,
Wherein, rx(m) it is to sample obtained discrete-time series, j is imaginary unit, and ω is power spectrum signal density, and m is
Discrete variable.
In the step 3), the calculation formula of data difference curve is：Data difference=real data-maximum entropy.
In the step 5), the calculation formula of EWL is：
W=k × h
Wherein, W is EWL, and k is coefficient correlation, and h is formation thickness, k=0.8.
In the step 5), retain the EWL in spectrum curve, invalid low ripple, the middle high trapezoidal shape of ripple are filtered
Go, so that curve becomes more smooth；Wherein 0.5WEWL~1.5WEWLFor EWL section, 0.25WEWL~
0.5WEWLAnd 1.5WEWL~1.75WEWLFor invalid low ripple and middle high ripple section.
In the step 6), the depth segment of abnormity point development, upper limit filtering are determined with the upper limit filter curve of each well
The depth segment in the range of 5 meters as abnormity point development up and down of depth point corresponding at the peak swing of curve to the left, then with
The amplitude maximum position of medium filtering curve corresponding to the depth segment of abnormity point development determines abnormity point.
For the present invention due to taking above technical scheme, it has advantages below：1st, the present invention is to the abnormity point in each well
Carry out com-parison and analysis, compared with traditional similar stratum control methods, it is small to this lithology difference of granite, vertical variations are small,
For the stratum that cyclicity does not protrude, there is stronger operability, more conventional Strata Comparison method has more practicality.2nd, originally
Invention is as a result of maximum entropy spectrum analysis method, therefore Strata Comparison quantitative extent is strengthened, and operability is stronger, Strata Comparison
Efficiency greatly reinforce.It is 3rd, of the invention due to using spectrum curve of the way of contrast by data difference Curve transform for spatial domain,
Therefore the characteristics of being suitable for this special reservoir of granite.
Brief description of the drawings
Fig. 1 is the schematic diagram on Division and contrast stratum of the present invention.
Embodiment
The present invention is described in detail with reference to the accompanying drawings and examples.
A kind of granitic formation Division and contrast method provided by the invention, it comprises the following steps：
1) borehole log data prepares
The curve that usual well logging personnel provide had done corresponding pretreatment, but due to processing method and the difference of custom
Fixed constant value -999 or -999.25 can be assigned in the start-stop position of curve and stitching position.And this negative anomaly can be serious
Influence the subsequent processing steps of the present invention, it is therefore desirable to first by the negative of the start-stop position of the log of each well and stitching position
Exceptional value is changed into 0；
2) log maximum entropy spectrum analysis
Handled by maximum entropy spectrum analysis function pair borehole log data, obtain the well logging spectrum curve of each well.Its
In, function and maximum entropy SX(ω) and correlation function rx(m) relation between determines by wiener-khintchine's theorem, i.e.,
Wherein, rx(m) it is to sample obtained discrete-time series, j is imaginary unit, and ω is power spectrum signal density, and m is
Discrete variable；
In the present embodiment, pretreated log data in step 1) is input to business software with las forms
In cyclolog softwares, previously given interval of interest top bottom depth and window width, maximum entropy spectrum analysis function pair is then called
Log data is handled, you can obtains the well logging spectrum curve of each well.
3) data difference curve is obtained
The maximum entropy of each depth point of the well logging spectrum curve of each well is calculated, and according to corresponding to each depth point
The real data of log, obtain the data difference curve of each well；Wherein, data difference=real data-maximum entropy.
4) it is spectrum curve by data difference Curve transform
The data difference curve negotiating iterative Simpson integral formula of each well in step 3) is subjected to mathematic integral conversion,
Reconstruct the spectrum curve of new each well.
Wherein, the process can borrow the Inpefa modules of business software cyclolog softwares to realize, comprise the following steps that：
Right mouse button is clicked on data sub-window, selects Processing → Spectral Attributes → INPEFA；Predicting
The interval of interest (default value is whole mouth well) of error filtering window selection analysis, and provide the size of prediction window (default value is
10 meters), you can complete conversion of the data difference curve to spectrum curve.
5) processing is filtered to the spectrum curve of reconstruct
According to every mouthful of well well logging in oil field to be measured and sidewall coring lithology statistical result obtain weathered stratum thickness the upper limit,
Intermediate value and lower limit, and according to the positive correlation of EWL and formation thickness obtain the upper limit of weathered stratum thickness, intermediate value and
The EWL of lower limit, i.e.,：WEWL=kCoefficient correlation×hFormation thickness, wherein, kCoefficient correlation=0.8；Selected according to the scope of EWL
Filtering parameter is taken, so as to be filtered processing to the spectrum curve of reconstruct, the EWL in reservation spectrum curve will be invalid
Low ripple, the trapezoidal shape of middle high ripple filter off, so as to obtain the three of the upper limit of each well, intermediate value and lower limit smooth filter curves.
Generally by (0.5WEWL~1.5WEWL) section is as pervasive EWL section, (0.25WEWL~0.5WEWL) and
(1.5WEWL~1.75WEWL) section is as low ripple and middle high ripple section.
In the present embodiment, the frequency spectrum using the ladder passband filter module in business software cyclolog softwares to reconstruct
Curve is filtered processing, and weathered stratum thickness is provided according to every mouthful of P oil fields well well logging and sidewall coring lithology statistical result
The upper limit is 160m, intermediate value 100m, lower limit 70m, shows how to set filter by taking weathered stratum thickness intermediate value 100m as an example below
Parameter in ripple module：W is asked for successively according to the formula in step 5)EWL=80, low ripple section is (20,40), EWL
Section is (40,120), and middle high ripple section is (120,140), low-cut wavelength in such filtration module
Amplitude 0 and the parameter items of Amplitude 100 choose 20m, 40m, i.e. Amplitude respectively<20 curve values whole mistake
Filter, 20<Amplitude<Curve values between 40 proportionally partly retain, and 40<Amplitude<120 whole reservations,
High-cut wavelength Ampl itude 0 and the parameter items of Amplitude 100 choose 120,140, i.e., 120 respectively<
Amplitude<Curve values between 140 proportionally partly retain, Amplitude>140 whole filterings.By that analogy, may be used
Three representative curve T_DT_ maximums, T_DT_ intermediate values and T_DT_ minimum values are obtained, reflects stratum longitudinal difference change.
6) abnormity point is determined
According to three filter curves of each well obtained after step 5) filtering process, exist in three curves of selection larger
The point of amplitude of variation comprises the following steps that as discrepancy：
1. the depth segment of abnormity point development is determined with " T_DT_ maximums " curve.At the peak swing of curve to the left pair
The depth point answered, in the range of upper and lower 5 meters, the depth segment as abnormity point development.Such as P1 wells, from shallow to depth, first is to the left
Depth point corresponding to the peak swing of side is 1290 meters, then the depth segment that 1285~1295 meters of depth segments are developed as abnormity point；
2. abnormity point is determined with " T_DT_ intermediate values " curve.According to the fine-characterization of curve, amplitude maximum position is determined, i.e.,
For abnormity point.Abnormity point in above-mentioned 1285~1295 meters of depth segments is 1294 meters.
7) abnormity point is numbered
The abnormity point of each well to determining in step 6) carries out vertical numbering from top to bottom, be designated as respectively T1, T2,
T3……。
8) classification and correlation is carried out
Each well key similar differences point on key section is sequentially connected, that is, completes granitic formation division and contrast.
In the present embodiment, by P1T1、P2T1It is connected, as first bed boundary；By P1T2、P2T2It is connected, as second
Bed boundary.So far, P oil fields are divided into more set stratum.
The various embodiments described above are merely to illustrate the present invention, wherein the structure of each part, connected mode etc. are all can be
Change, every equivalents carried out on the basis of technical solution of the present invention and improvement, it should not exclude the present invention's
Outside protection domain.
Claims (5)
1. a kind of granitic formation Division and contrast method, it comprises the following steps：
1) borehole log data prepares：By the start-stop position of log of each well in oil field to be measured and the negative anomaly of stitching position
Value is changed into 0；
2) log maximum entropy spectrum analysis：Handled by the borehole log data of each well of maximum entropy spectrum analysis function pair,
Obtain the well logging spectrum curve of each well；
3) data difference curve is obtained：Calculate each well well logging spectrum curve each depth point maximum entropy, and according to
The real data of log corresponding to each depth point, obtain the data difference curve of each well；
4) it is spectrum curve by data difference Curve transform：By the data difference curve negotiating iterative Simpson integral formula of each well
Mathematic integral conversion is carried out, reconstructs the spectrum curve of each well；
5) processing is filtered to the spectrum curve of reconstruct：Counted according to every mouthful of well well logging and sidewall coring lithology in oil field to be measured
As a result the upper limit, intermediate value and the lower limit of weathered stratum thickness are obtained, and is obtained according to the positive correlation of EWL and formation thickness
The EWL of the upper limit of weathered stratum thickness, intermediate value and lower limit is taken, the scope according to EWL chooses filtering parameter, so as to
Processing is filtered to the spectrum curve of reconstruct, obtains three filter curves of the upper limit of each well, intermediate value and lower limit；
6) abnormity point is determined：According to three filter curves of each well obtained after filtering process, choose in three filter curves
In the presence of large change amplitude point as abnormity point；
7) abnormity point is numbered：The abnormity point of each well to being determined in step 6) carries out vertical numbering from top to bottom；
8) classification and correlation is carried out：Each well key similar differences point on key section is sequentially connected, that is, with completing granite
Layer division and contrast.
A kind of 2. granitic formation Division and contrast method as claimed in claim 1, it is characterised in that：In the step 2),
Function and maximum entropy SXThe discrete-time series r that (ω) and sampling obtainx(m) relation between determines by wiener-khintchine's theorem, i.e.,
<mrow>
<msub>
<mi>S</mi>
<mi>X</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mi>&infin;</mi>
</msubsup>
<msub>
<mi>r</mi>
<mi>x</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>&omega;</mi>
<mi>m</mi>
</mrow>
</msup>
</mrow>
Wherein, j is imaginary unit, and ω is power spectrum signal density, and m is discrete variable.
A kind of 3. granitic formation Division and contrast method as claimed in claim 1, it is characterised in that：In the step 3),
The calculation formula of data difference curve is：Data difference=real data-maximum entropy.
A kind of 4. granitic formation Division and contrast method as claimed in claim 1, it is characterised in that：In the step 5),
The calculation formula of EWL is：
W=k × h
Wherein, W is EWL, and k is coefficient correlation, and h is formation thickness, k=0.8.
A kind of 5. granitic formation Division and contrast method as claimed in claim 1, it is characterised in that：In the step 6),
The depth segment of abnormity point development is determined with the upper limit filter curve of each well, at the peak swing of upper limit filter curve to the left pair
The depth segment in the range of 5 meters as abnormity point development up and down for the depth point answered, then with corresponding to the depth segment of abnormity point development
The amplitude maximum position of medium filtering curve determine abnormity point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610403178.XA CN106094051B (en) | 2016-06-08 | 2016-06-08 | A kind of granitic formation Division and contrast method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610403178.XA CN106094051B (en) | 2016-06-08 | 2016-06-08 | A kind of granitic formation Division and contrast method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106094051A CN106094051A (en) | 2016-11-09 |
| CN106094051B true CN106094051B (en) | 2017-12-15 |
Family
ID=57229078
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610403178.XA Active CN106094051B (en) | 2016-06-08 | 2016-06-08 | A kind of granitic formation Division and contrast method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106094051B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112302623B (en) * | 2020-10-31 | 2023-07-28 | 中国海洋石油集团有限公司 | Granite down-the-hole weathering crust structure division method based on logging data |
| CN114215513B (en) | 2022-02-21 | 2022-05-10 | 中海油研究总院有限责任公司 | Quantitative discrimination method, device, medium and equipment for buried hill reservoir mode |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1997039367A1 (en) * | 1996-04-12 | 1997-10-23 | Amoco Corporation | Method and apparatus for seismic signal processing and exploration |
| US6092026A (en) * | 1998-01-22 | 2000-07-18 | Bp Amoco Corporation | Seismic signal processing and exploration |
| CN103529474B (en) * | 2012-07-02 | 2016-05-04 | 中国石油化工股份有限公司 | Adopt lithology segmentation to realize the method for petrofacies fine description |
| CN104360386B (en) * | 2014-12-04 | 2016-09-21 | 中国海洋石油总公司 | A kind of panel method for granitic formation Division and contrast |
-
2016
- 2016-06-08 CN CN201610403178.XA patent/CN106094051B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106094051A (en) | 2016-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113759424B (en) | Karst reservoir filling analysis method and system based on spectral decomposition and machine learning | |
| CN109709607B (en) | Method and device for predicting thickness of thin sandstone reservoir | |
| CN102129084B (en) | Method and device for acquiring seismic thin reservoir speed through well control | |
| CN106842307B (en) | Reservoir stratum fine prediction method based on waveform classification and retrieval under forward constraint | |
| CN105863628B (en) | A kind of phase of development subtle hydrocarbon reservoir method | |
| CN109425900A (en) | A kind of Seismic Reservoir Prediction method | |
| CN111983683B (en) | Prediction method and system for lake-facies limestone reservoir under low-well condition | |
| CN105301647B (en) | The method for distinguishing grey matter mud stone and sandstone | |
| CN107462924A (en) | A kind of absolute wave impedance inversion method independent of well-log information | |
| CN103439740A (en) | Method and device for predicting relative impedance based on dipole seismic wavelet multiple integral | |
| CN117150929A (en) | Deep oil gas accurate navigation well earthquake combined clastic rock lithology model updating system | |
| CN106094051B (en) | A kind of granitic formation Division and contrast method | |
| CN106405644A (en) | Crack determination method and apparatus | |
| CN105093288B (en) | A kind of diffracted wave separation method based on kinematics wave field attributes | |
| CN114482995B (en) | Fine determination method for clay content of fine sediment | |
| CN113419274B (en) | Three-dimensional seismic slice attribute body extraction method based on high-precision sequence grid model | |
| CN111060967A (en) | Well control velocity field modeling method | |
| CN117251669A (en) | A method to analyze paleoclimate evolution trends based on lithological cycles | |
| CN112835098B (en) | Method and device for predicting energy storage coefficient of weathered-crust karst reservoir | |
| CN114428288B (en) | A quantitative prediction method for thin interbedded sandy beach reservoirs | |
| CN115267937A (en) | Braided river delta sand body configuration prediction method and system under offshore thin well pattern condition | |
| CN116104470A (en) | Geosteering formation identification and prediction method, system, equipment and storage medium | |
| CN110389381B (en) | Sand reservoir prediction method and device based on seismic attributes | |
| CN105891887A (en) | Velocity longitudinal and transverse high density analysis method based on stacked data | |
| CN115343781B (en) | A complex structural velocity modeling method based on structural constraints using well-seismic combined methods |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC Research Institute Patentee before: China National Offshore Oil Corporation |
|
| CP01 | Change in the name or title of a patent holder |