CN106556870B - A kind of well logging modeling method and system based on holding edge - Google Patents
A kind of well logging modeling method and system based on holding edge Download PDFInfo
- Publication number
- CN106556870B CN106556870B CN201510617079.7A CN201510617079A CN106556870B CN 106556870 B CN106556870 B CN 106556870B CN 201510617079 A CN201510617079 A CN 201510617079A CN 106556870 B CN106556870 B CN 106556870B
- Authority
- CN
- China
- Prior art keywords
- log
- edge
- thickness
- thin layer
- sliding window
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Include based on the well logging modeling method for keeping edge and system, this method the invention discloses a kind of:Log is inputted, the thickness of well logging most thin layer is obtained according to existing result of log interpretation;According to the length for selecting sliding window by the thickness of the most thin layer acquired in result of log interpretation, edge holding processing is carried out to log;According to the thickness by the most thin layer acquired in result of log interpretation, to keeping that treated through edge, log carries out 1D modelings.This method keeps algorithm to pre-process log using edge, not only smooth log but also the boundary for remaining layer position, and the thickness results for having well log interpretation most thin layer are combined, select best sliding window so that electrobed and modeling are more reasonable, accurate.
Description
Technical field
It is specifically a kind of based on the well logging modeling side for keeping edge the invention belongs to geophysical prospecting for oil technical field
Method and system.
Background technology
Conventional well logging modeling method is according to well logging measured curve be layered after directly sampling or standard are smooth
And modeling.Since there is log the feature of " high frequency " shake, directly sampling or smooth post-sampling will produce certain
Error, meanwhile, when underground medium thin interbed is relatively developed, standard will smoothly scabble the amplitude of two side datas, and be unfavorable for thin layer
It is fine division and modeling.
For the feature of log " high frequency " shake and the development of underground thin interbed, it is desirable to provide a kind of raising well logging point
The method of layer and modeling accuracy.
Invention content
The first technical problem to be solved by the present invention is to need to provide a kind of well logging modeling method based on holding edge,
This method can improve the reliability and accuracy of electrobed and modeling.In addition, the present invention also provides one kind based on holding
The well logging modeling at edge.
In order to solve the above-mentioned technical problem, embodiments herein provide firstly it is a kind of based on keep edge well logging build
Mould method, including:Log is inputted, the thickness of well logging most thin layer is obtained according to existing result of log interpretation;According to by logging well
The thickness of most thin layer acquired in explanation results selectes the length of sliding window, and edge holding processing is carried out to log;Root
According to the thickness by the most thin layer acquired in result of log interpretation, to keeping that treated through edge, log carries out 1D modelings.
Preferably, in carrying out the step of edge keeps processing to log, point is calculated for each of log,
The mean value of each sliding window comprising the calculating point and corresponding variance are calculated, the flat of the sliding window of variance yields minimum is taken
Value of the mean value as the calculating point.
Preferably, according to the thickness by the most thin layer acquired in result of log interpretation, to being kept through edge, treated
Log carried out in the step of 1D modelings, using the thickness of the most thin layer acquired in result of log interpretation as the sampling interval, to warp
Edge keeps treated log to carry out longitudinal dispersion sampling, and carries out uniform assignment in the horizontal, and then establishes stratiform
Model.
Preferably, the mean value D of sliding window m is calculated using following expressionmAnd variances sigmam,
Wherein, DmFor sample value D in sliding window miAverage value, n be sample value number.
Preferably, the log includes sound wave, density, neutron, resistivity and velocity logging curve.
Another aspect according to the ... of the embodiment of the present invention additionally provides a kind of based on the well logging modeling for keeping edge, packet
It includes:Thickness acquisition module inputs log, and the thickness of well logging most thin layer is obtained according to existing result of log interpretation;Edge
Processing module is kept, according to the length for selecting sliding window by the thickness of the most thin layer acquired in result of log interpretation, to surveying
Well curve carries out edge holding processing;Modeling module, according to the thickness by the most thin layer acquired in result of log interpretation, to warp
Edge keeps treated log to carry out 1D modelings.
Preferably, the edge keeps processing module further, and point is calculated for each of log, and it includes to be somebody's turn to do to calculate
The mean value of each sliding window of point and corresponding variance are calculated, taking the average value of the sliding window of variance yields minimum to be used as should
Calculate the value of point.
Preferably, the modeling module is further, is between sampling with the thickness of the most thin layer acquired in result of log interpretation
Every to keeping treated log to carry out longitudinal dispersion sampling through edge, and carrying out uniform assignment, Jin Erjian in the horizontal
Vertical stratified model.
Preferably, the edge holding processing module calculates the mean value D of sliding window m using following expressionmThe side and
Poor σm,
Wherein, DmFor sample value D in sliding window miAverage value, n be sample value number.
Compared with prior art, one or more of said program embodiment can have the following advantages that or beneficial to effect
Fruit.
Relative to existing technology, the method for the embodiment of the present invention keeps algorithm to locate log in advance using edge
Reason, not only smooth log remains the boundary of layer position again, and combines the thickness results for having well log interpretation most thin layer, selected
Best sliding window so that electrobed and modeling are more rationally, accurately.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
It obtains it is clear that being understood by implementing technical scheme of the present invention.The purpose of the present invention and other advantages can by
Specifically noted structure and/or flow are realized and are obtained in specification, claims and attached drawing.
Description of the drawings
Attached drawing is used for providing to the technical solution of the application or further understanding for the prior art, and constitution instruction
A part.Wherein, technical side of the attached drawing of the embodiment of the present application together with embodiments herein for explaining the application is expressed
Case, but do not constitute the limitation to technical scheme.
Fig. 1 is the flow diagram based on the well logging modeling method for keeping edge of the embodiment of the present invention.
Fig. 2 is certain work area practical logging curve of one example of the present invention and the curve pretreatment knot of different smooth manners
Fruit.
Fig. 3 is the 1D well logging modeling results of one example of the present invention.
Fig. 4 is the structural schematic diagram based on the well logging modeling for keeping edge of the embodiment of the present invention
Specific implementation mode
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, how to be applied to the present invention whereby
Technological means solves technical problem, and the realization process for reaching relevant art effect can fully understand and implement.This Shen
Each feature that please be in embodiment and embodiment, can be combined with each other under the premise of not colliding, be formed by technical solution
Within protection scope of the present invention.
In addition, the flow of attached drawing can be in the computer system of such as a group of computer-executable instructions the step of illustrating
Middle execution.Also, although logical order is shown in flow charts, and it in some cases, can be with different from herein
Sequence executes shown or described step.
(embodiment one)
Conventional electrobed and modeling method be directly sampled based on log or standard it is smooth after be layered.For
Improve the precision of electrobed and modeling, the embodiment of the present invention is a kind of based on the well logging modeling method for keeping edge, the party
Method can improve the reliability and accuracy of electrobed and modeling.This method carries out holding edge to log first
Reason, then according to treated, log carries out 1D modelings, to realize that the high-precision based on log models.
Fig. 1 is the flow diagram based on the well logging modeling method for keeping edge of the embodiment of the present application.Below with reference to Fig. 1
Carry out each step of process in detail.
In step s 110, log is inputted, the thickness of well logging most thin layer is obtained according to existing result of log interpretation, with
Selection for sliding window when keeping edge treated.
Specifically, log in the present embodiment is this field Logging Curves, Ke Yishi:Sound wave, density, in
Son, resistivity or velocity logging curve etc..And it can be reservoir thickness explanation results to have result of log interpretation, according to reservoir
Thickness interpretation result can directly acquire the thickness of well logging most thin layer.
In the step s 120, according to the length for selecting sliding window by the thickness of the most thin layer acquired in result of log interpretation
Degree carries out edge holding processing to log.
The basic principle kept the following detailed description of edge.
Refer to a log (this method is suitable for all logs) for 1D logs, sets length of window
And point is each calculated for the 1D logs, calculate the mean value and corresponding under the sliding window for including the calculating point
Variance takes value of the average value of the sliding window of variance yields minimum as the calculating point.Since the value of the calculating point includes only side
The data of edge side, therefore both remained the presence at edge and be not extended, substituted with hourly value realizes log again
Smoothly.
Steps are as follows for the specific implementation that edge is kept:
(1) by taking central point i as an example, choosing multiwindow, (such as 4 points, using 4 points as window, slided, and all includes
Including heart point i)
window-1:(Di-3,Di-2,Di-1,Di), m=1, to calculate three sample values that calculation window includes on the left of point i;
window-2:(Di-2,Di-1,Di,Di+1), m=2, to calculate two sample values that calculation window includes on the left of point i
A sample value for including with right side calculation window;
window-3:(Di-1,Di,Di+1,Di+2), m=3, to calculate two sample values that calculation window includes on the right side of point i
A sample value for including with left side calculation window;
window-4:(Di,Di+1,Di+2,Di+3), m=4, to calculate three sample values that calculation window includes on the right side of point i;
Wherein each window must include to calculate point Di, m expression window serial numbers;Di+1、Di+2、Di+3、Di-1、Di-2、Di-3
Indicate i points or so corresponding neighborhood calculation point.
(2) mean value and variance of each window are calculated
Wherein, DmFor the middle sample value D of sliding window miAverage value, n be sample value number.The expression formula is namely
The average value D of log data is calculated separately on each window (Window-m, m=1,2,3,4)mAnd variances sigmam;
(3) use the average value of the window with minimum variance in multiwindow as the value of the calculating point, i.e.,:
Di={ Dm, work as σmFor minimum duration } (2)
In step s 130, according to the thickness by the most thin layer acquired in result of log interpretation, to being handled through edge holding
Log afterwards carries out 1D modelings.
Specifically, using the thickness of the most thin layer acquired in result of log interpretation as the sampling interval, to being handled through edge holding
Log afterwards carries out longitudinal dispersion sampling, and carries out uniform assignment in the horizontal, and then establishes stratified model.
Relative to existing technology, the method for the embodiment of the present invention keeps algorithm to locate log in advance using edge
Reason, not only smooth log remains the boundary of layer position again, and combines the thickness results for having well log interpretation most thin layer, selected
Best sliding window so that electrobed and modeling are more rationally, accurately.
Example:
Certain work area practical logging curve is modeled below with the method for the present invention to illustrate to calculate based on holding edge
The validity of method.
Input log (lighter curve in Fig. 2) is by taking speed as an example first, and collects result of log interpretation, includes mainly
The thickness of well log interpretation most thin layer.The work area target zone thin interbed is relatively developed, and most thickness of thin layer is 6m.Then 6m conducts are selected
Sliding window is carried out the pretreatment kept based on arithmetic mean, harmonic average, edge to rate curve, obtains different pretreatments
Rate curve after mode.
Lighter curve is primitive curve in Fig. 2, and black curve is pretreated curve, and specifically, black is bent in left figure
Line is arithmetic mean sharpening result, and black curve is harmonic average sharpening result in middle figure, and black curve is to keep side in right figure
Edge sharpening result.It can be seen that the smooth manner relative to arithmetic mean and harmonic average, keeps edge algorithms both effectively to keep
The edge of layer position, while the data of both sides of edges have been carried out smoothly, therefore, effectively increase the essence of electrobed and modeling
Degree.
Next according to well logging holding edge filter as a result, choosing 6m as the sampling interval, to log progress 1D
Modeling, from figure 3, it can be seen that in interval of interest 3500m or so, it is more clear, accurate to portraying for thin sand-mud interbed.Therefore originally
The 1D well logging modeling results that inventive method obtains have higher precision.
(embodiment two)
Fig. 4 is the structural schematic diagram based on the well logging modeling for keeping edge of the embodiment of the present invention, below with reference to Fig. 4
To illustrate each composed structure of this system.
In the present embodiment, well logging modeling includes following components:Thickness acquisition module 41, edge keep processing mould
Block 42 and modeling module 43.Wherein, thickness acquisition module 41 keeps processing module 42 to connect with edge, and edge keeps processing module
42 connect with modeling module 43.
Referring to FIG. 4, when being modeled to log, at the thickness acquisition module 41 of the present embodiment, edge holding
Reason module 42 and modeling module 43 execute the operation of the step 110, step 120 and step 130 of embodiment one respectively, herein no longer
Explanation is developed in details.
Those skilled in the art should be understood that each of the device that the above-mentioned embodiment of the present invention is provided and/or system
Each step in component part and method, they can be concentrated on a single computing device, or be distributed in multiple calculating
On the network that device is formed.Optionally, they can be realized with the program code that computing device can perform.It is thus possible to
It is stored in storage device and is performed by computing device, or they are fabricated to each integrated circuit modules,
Or by them multiple modules or step be fabricated to single integrated circuit module and realize.In this way, the present invention is not limited to
Any specific hardware and software combines.
Although disclosed herein embodiment it is as above, the content only technical solution of the present invention for ease of understanding
And the embodiment used, it is not limited to the present invention.Technical staff in any fields of the present invention is not departing from this
Under the premise of the disclosed spirit and scope of invention, any modification and variation can be carried out in the form and details of implementation,
But the scope of patent protection of the present invention, still should be subject to the scope of the claims as defined in the appended claims.
Claims (5)
1. it is a kind of based on the well logging modeling method for keeping edge, including:
Log is inputted, the thickness of well logging most thin layer is obtained according to existing result of log interpretation;
According to the length for selecting sliding window by the thickness of the most thin layer acquired in result of log interpretation, side is carried out to log
Edge holding is handled;
According to the thickness by the most thin layer acquired in result of log interpretation, to keeping that treated through edge, log carries out 1D
Modeling;
Wherein, in carrying out the step of edge keeps processing to log, point is calculated for each of log, calculates packet
The mean value and corresponding variance of each sliding window containing the calculating point take the average value of the sliding window of variance yields minimum to make
For the value of the calculating point;
According to the thickness by the most thin layer acquired in result of log interpretation, to keeping treated log to carry out through edge
In the step of 1D is modeled, using the thickness of the most thin layer acquired in result of log interpretation as the sampling interval, to being handled through edge holding
Log afterwards carries out longitudinal dispersion sampling, and carries out uniform assignment in the horizontal, and then establishes stratified model.
2. according to the method described in claim 1, it is characterized in that, calculating the mean value D of sliding window m using following expressionm
And variances sigmam,
Wherein, DmFor sample value D in sliding window miAverage value, n be sample value number.
3. according to the method described in claim 1, it is characterized in that,
The log includes sound wave, density, neutron, resistivity and velocity logging curve.
4. it is a kind of based on the well logging modeling for keeping edge, including:
Thickness acquisition module inputs log, and the thickness of well logging most thin layer is obtained according to existing result of log interpretation;
Edge keeps processing module, according to the length for selecting sliding window by the thickness of the most thin layer acquired in result of log interpretation
Degree carries out edge holding processing to log;
Modeling module, according to the thickness by the most thin layer acquired in result of log interpretation, to keeping treated to survey through edge
Well curve carries out 1D modelings;
Wherein, the edge keeps processing module further, and point is calculated for each of log, and it includes the calculating point to calculate
Each sliding window mean value and corresponding variance, take the average value of the sliding window of variance yields minimum as the calculating point
Value;
The modeling module is further, using the thickness of the most thin layer acquired in result of log interpretation as the sampling interval, to through edge
It keeps treated log to carry out longitudinal dispersion sampling, and carries out uniform assignment in the horizontal, and then establish stratified model.
5. system according to claim 4, which is characterized in that the edge keep processing module using following expression come
Calculate the mean value D of sliding window mmAnd variances sigmam,
Wherein, DmFor sample value D in sliding window miAverage value, n be sample value number.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510617079.7A CN106556870B (en) | 2015-09-24 | 2015-09-24 | A kind of well logging modeling method and system based on holding edge |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510617079.7A CN106556870B (en) | 2015-09-24 | 2015-09-24 | A kind of well logging modeling method and system based on holding edge |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106556870A CN106556870A (en) | 2017-04-05 |
CN106556870B true CN106556870B (en) | 2018-08-31 |
Family
ID=58415552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510617079.7A Active CN106556870B (en) | 2015-09-24 | 2015-09-24 | A kind of well logging modeling method and system based on holding edge |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106556870B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108825223B (en) * | 2018-06-22 | 2021-10-22 | 西南石油大学 | Method for extracting bedding characteristics of shale stratum |
CN109085646A (en) * | 2018-10-18 | 2018-12-25 | 中国海洋石油集团有限公司 | A kind of stacked sand body recognition methods of the delta facies based on EPS phase body attribute |
-
2015
- 2015-09-24 CN CN201510617079.7A patent/CN106556870B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN106556870A (en) | 2017-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104297785B (en) | Petrofacies constraint reservoir physical parameter inversion method and device | |
CN105425289B (en) | The method and apparatus for determining low frequency wave impedance | |
CN112784980B (en) | Intelligent logging horizon dividing method | |
CN108897042A (en) | Content of organic matter earthquake prediction method and device | |
CN107894618B (en) | A kind of full waveform inversion gradient preprocess method based on model smoothing algorithm | |
CN107783185A (en) | A kind of processing method and processing device of tomographic statics | |
CN113326852A (en) | Model training method, device, equipment, storage medium and program product | |
CN106556870B (en) | A kind of well logging modeling method and system based on holding edge | |
CN110297273A (en) | The prestack inversion method and device of reservoir | |
CN110879412A (en) | Underground transverse wave velocity inversion method, device, computing equipment and storage medium | |
CN106291682A (en) | A kind of poststack sound impedance inversion method based on base method for tracing | |
CN105986819B (en) | The method and apparatus with integrated interpretation are automatically processed for well-log information | |
CN105353409B (en) | A kind of method and system for full waveform inversion focus to be inhibited to encode cross-talk noise | |
CN109655890A (en) | A kind of shallow mid-deep strata joint chromatography inversion speed modeling method of Depth Domain and system | |
CN107870361A (en) | A kind of earthquake diving Wave chromatography imaging method, device and terminal device | |
CN112946746B (en) | Method and device for improving AVO inversion accuracy of thin coal seam | |
CN109490954B (en) | Wave field forward modeling method and device | |
CN106950601A (en) | Static correcting method and device | |
CN109633781B (en) | Geological property acquisition method and device, electronic equipment and storage medium | |
CN112444850A (en) | Seismic data velocity modeling method, storage medium, and computing device | |
CN113866827B (en) | Interpretation velocity modeling seismic imaging method, system, medium and equipment | |
CN115659773A (en) | Full waveform inversion acceleration method based on depth network and related device | |
CN112016956A (en) | BP neural network-based ore grade estimation method and device | |
CN105137480A (en) | Method and apparatus for determining pre-stack time migration velocity of VTI medium | |
CN114444555A (en) | U-net neural network-based seismic horizon interpretation method and device |
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 |