CN110390702A - The change face element plane interpolation method of analysis is moved towards based on geology - Google Patents
The change face element plane interpolation method of analysis is moved towards based on geology Download PDFInfo
- Publication number
- CN110390702A CN110390702A CN201910618381.2A CN201910618381A CN110390702A CN 110390702 A CN110390702 A CN 110390702A CN 201910618381 A CN201910618381 A CN 201910618381A CN 110390702 A CN110390702 A CN 110390702A
- Authority
- CN
- China
- Prior art keywords
- geology
- analysis
- face element
- directional derivative
- interpolation
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
Abstract
The present invention provides a kind of change face element plane interpolation method that analysis is moved towards based on geology, comprising: step 1, input plane scatterplot data；Step 2, the planar point set data inputted for step 1 require defined analysis grid according to analysis precision, are determined as geology trend by calculating different direction minimum directional derivative in analysis grid；Step 3, it moves towards to design oval interpolation face element for long axis direction with geology obtained in step 2；Step 4, it is analysis face element with interpolation ellipse defined in step 3, carries out the planar point set data interpolating of geology direction constraint；Step 5, using interpolated data obtained in step 4, sketching isopleth map and carries out geological verification.This makes prediction result be more in line with geologic rule based on the change face element plane interpolation method that geology moves towards analysis, improves geological research personnel and further analyzes the accuracy for solving geological problem using prediction result.
Description
Technical field
The present invention relates to geology and geostatistical techniques field, especially relate to a kind of move towards analysis based on geology
Become face element plane interpolation method.
Background technique
In geostatistics, predict to use when any point geological property using known geological property limited in plane
To various plane interpolation methods, existing interpolation method, which usually passes through, defines interpolation radius R, and utilizing by the center of circle, R of interpolation point is half
Known point in the round interpolation face element of diameter is constraint condition, calculates the geological property value of interpolation point.This interpolation method lacks
Be trapped in seldom consider geological property it is generally existing walk tropism, i.e., have very difference greatly in the speed of different directions geological property value variation
It is different；Meanwhile any quantitative method no for the determination of interpolation radius R cannot abundant body often with certain randomness
The geology correlation between scatterplot data is known.
Thus we have invented a kind of new change face element plane interpolation method for moving towards analysis based on geology, solve above
Technical problem.
Summary of the invention
The object of the present invention is to provide a kind of can be improved to carry out geological property planar interpolation using known scatterplot data
Accuracy and the rational change face element plane interpolation method that analysis is moved towards based on geology.
The purpose of the present invention can be achieved by the following technical measures: the change face element planar interpolation of analysis is moved towards based on geology
Method, it includes: step 1, input plane scatterplot data that the change face element plane interpolation method of analysis should be moved towards based on geology；Step 2,
For the planar point set data that step 1 inputs, defined analysis grid is required according to analysis precision, passes through calculating in analysis grid
Different direction minimum directional derivative is determined as geology trend；Step 3, with the trend of geology obtained in step 2 for long axis side
To the oval interpolation face element of design；Step 4, it is analysis face element with interpolation ellipse defined in step 3, carries out geology direction about
The planar point set data interpolating of beam；Step 5, using interpolated data obtained in step 4, sketching isopleth map and carries out geology and tests
Card.
The purpose of the present invention can be also achieved by the following technical measures:
In step 1, according to tri- array structure input plane geological property scatterplot data of X, Y, A.
In step 2, according to required precision defined analysis grid, 0-180 ° of side of certain step size computation is pressed in analysis grid
The directional derivative of position input scatterplot data, records minimum directional derivative and maximum directional derivative size, orientation.
In step 2, the determination of geology trend is completed in accordance with the following methods:
If a=f (x, y) is defined in certain field of ground point P (x, y), direct north azimuth is 0 °；Injection line l,
It is with direct north angleThen P point is defined as along the directional derivative of ray l
For directional derivative expression geological property a from P point along the change rate in the direction ray l, minimum directional derivative corresponds to azimuthRepresent geology trend.
In step 3, orientation is corresponded to according to the minimum directional derivative being calculated in step 2 and determines transverse direction；
According to long axis direction directional derivative and short-axis direction directional derivative size relation, the ratio between transverse and short axle are determined；According to edge
The nearest extreme point of second derivative determines the size of long axis on the direction ray l.
In step 4, using oval interpolation face element as analysis face element, the method for defining ellipsometric analysis face element:
Oval direction definition is that minimum directional derivative corresponds to azimuth
Oval eccentricity is acquired using long axis direction directional derivative corresponding with short-axis direction:
If long axis and short-axis direction directional derivative are respectively a, b, enableThen oval eccentricity e=c/a；
For the plane geological property scatterplot data inputted in step 1, certain step size computation 0- is pressed in analysis grid
The directional derivative of 180 ° of orientation input scatterplot data, acquires minimum directional derivative；The grid number constituted for minimum directional derivative
According to moving towards to carry out secondary derivation along geology, first extreme point P ˊ (x ˊ, y ˊ) of second derivative in this direction is determined, point of origin P (x,
Y) transverse length a is defined as with the distance between extreme point.
In step 4, for the arbitrary point in analysis grid, the direction of ellipsometric analysis face element used in interpolation calculation,
Eccentricity and long axis size be it is variable, i.e., utmostly guarantee participate in interpolation calculation known scatterplot between have preferably
Qualitative correlation.
In steps of 5, for the arbitrary point on analysis grid defined in step 2, with ellipse obtained in step 3 point
Face element is analysed as constraint, planar interpolation is carried out, is sketched isopleth using interpolation result, and carry out geological verification.
The change face element plane interpolation method that analysis is moved towards based on geology in the present invention is existed by analyzing known scatterplot data
Variation Features in plane define the variable oval interpolation analysis face element of a size, direction and compression, are constrained pair with this
The known scatterplot of interpolation point geological property constraint, has fully considered the geology phase between the directionality and scatterplot data of geologic data
Guan Xing makes prediction result be more in line with geologic rule, improves geological research personnel using prediction result and further analyzes solution
The accuracy of geological problem.
This is based on the direction that the change face element plane interpolation method that geology moves towards analysis includes by Calculation Plane scatterplot data
Derivative analyzes the part trend of geologic objective；And move towards to determine the direction of oval interpolation face element according to geology, and according to geology
The proportionate relationship of trend and the directional derivative on perpendicular direction calculates the eccentricity of oval interpolation face element；It is ellipse according to interpolation
First extreme point position of circle long axis direction second derivative determines long axis size；It is carried out in oval interpolation face element based on this flat
The change face element planar interpolation of face geology scatterplot data generates grid data.This is inserted based on the change face element plane that geology moves towards analysis
Value method moves towards size, eccentricity and the long axis direction of determining oval interpolation face element by automatically analyzing geological object plane, makes
The grid data that interpolation obtains is more in line with geologic rule, and so that geology isopleth is more in line with geology reality at figure and answer with research
With the habit of personnel, " buphthalmos " phenomenon that general isoplethes drawing software is easy to appear at figure is effectively prevented, machine is overcome
The defect of lacking direction property of drawing isoline improves supporting function of the graph in geological research.
Detailed description of the invention
Fig. 1 is the process of a specific embodiment of the change face element plane interpolation method of the invention that analysis is moved towards based on geology
Figure；
Fig. 2 is the conventional isopleth map determining radius circular interpolation face element Kriging regression method side and sketching；
Fig. 3 is the isopleth map that geology provided by the invention moves towards that constraint change face element interpolation method is sketched.
Specific embodiment
To enable above and other objects, features and advantages of the invention to be clearer and more comprehensible, preferably implementation is cited below particularly out
Example, and cooperate shown in attached drawing, it is described in detail below.
Step 1, two-dimensional surface scatterplot data are inputted；
Step 2, defined analysis grid is required according to analysis precision, by calculating different direction minimum side in analysis grid
Determine that interpolation point geology is moved towards to derivative；According to tri- array structure input plane geological property scatterplot data of X, Y, A, wanted according to precision
Defined analysis grid is sought, by the directional derivative of 0-180 ° of orientation input scatterplot data of certain step size computation, note in analysis grid
Record minimum directional derivative orientation.The determination of geology trend is completed in accordance with the following methods:
If a=f (x, y) is defined in certain field of ground point P (x, y), direct north azimuth is 0 °；Injection line l,
It is with direct north angleThen P point is defined as along the directional derivative of ray l
For directional derivative expression geological property a from P point along the change rate in the direction ray l, minimum directional derivative corresponds to azimuthRepresent geology trend.
Step 3, it moves towards to design oval interpolation face element for long axis direction with geology obtained in step 2；
Move towards to design oval interpolation face element for long axis direction with geology obtained in step 2, according to long axis direction with it is short
Directional derivative size in axis direction determines long axis and short axle ratio, is determined according to the first extreme point of long axis direction second derivative
Long axis size；
Step 4, it is analysis face element with the interpolation ellipse determined in step 3, carries out the planar point set number of geology direction constraint
According to interpolation.
The method for defining ellipsometric analysis face element are as follows:
Oval direction definition is that minimum directional derivative corresponds to azimuth
Oval eccentricity is acquired using long axis direction directional derivative corresponding with short-axis direction:
If long axis and short-axis direction directional derivative are respectively a, b, enableThen oval eccentricity e=c/a.
For the plane geological property scatterplot data inputted in step 1, certain step size computation 0- is pressed in analysis grid
The directional derivative of 180 ° of orientation input scatterplot data, acquires minimum directional derivative；The grid number constituted for minimum directional derivative
According to moving towards to carry out secondary derivation along geology, first extreme point P ˊ (x ˊ, y ˊ) of second derivative in this direction is determined, point of origin P (x,
Y) transverse length a is defined as with the distance between extreme point.
For analyzing the arbitrary point in grid, direction, eccentricity and the length of ellipsometric analysis face element used in interpolation calculation
Axis size be it is variable, i.e., utmostly guarantee participate in interpolation calculation known scatterplot between have preferable geology correlation.
Step 5, using interpolated data obtained in step 4, sketching isopleth map and carries out geological verification.
For defined in step 2 analysis grid on arbitrary point, using ellipse conformal analysis face element obtained in step 3 as
Constraint carries out planar interpolation, is sketched isopleth using interpolation result, and carry out geological verification.
As shown in FIG. 1, FIG. 1 is the flow charts of the change face element plane interpolation method of the invention that analysis is moved towards based on geology.
In step 101, input plane scatterplot data.Process enters step 102.
In step 102, for the planar point set data that step 101 inputs, defined analysis grid is required according to analysis precision,
It is determined as geology trend by calculating different direction minimum directional derivative in analysis grid.Process enters step 103.
In step 103, move towards to design oval interpolation face element for long axis direction with geology obtained in step 102, according to
Directional derivative size on long axis direction and short-axis direction, determines long axis and short axle ratio, according to the of long axis direction second derivative
One extreme point determines long axis size.Process enters step 104.
In step 104, it is analysis face element with interpolation ellipse defined in step 103, carries out the plane of geology direction constraint
Scatterplot data interpolating.Process enters step 105.
In step 105, using interpolated data obtained in step 104, the isopleth map for sketching as shown in Figure 3 with such as Fig. 2 institute
The isopleth map comparison sketched by Kriging regression method shown, and carry out geological verification.
Claims (8)
1. moving towards the change face element plane interpolation method of analysis based on geology, which is characterized in that the change of analysis should be moved towards based on geology
Face element plane interpolation method includes:
Step 1, input plane scatterplot data；
Step 2, the planar point set data inputted for step 1 require defined analysis grid according to analysis precision, in analysis grid
It is interior to be determined as geology trend by calculating different direction minimum directional derivative；
Step 3, it moves towards to design oval interpolation face element for long axis direction with geology obtained in step 2；
It step 4, is analysis face element with interpolation ellipse defined in step 3, the planar point set data for carrying out the constraint of geology direction are inserted
Value；
Step 5, using interpolated data obtained in step 4, sketching isopleth map and carries out geological verification.
2. the change face element plane interpolation method according to claim 1 for moving towards analysis based on geology, which is characterized in that in step
In rapid 1, according to tri- array structure input plane geological property scatterplot data of X, Y, A.
3. the change face element plane interpolation method according to claim 1 for moving towards analysis based on geology, which is characterized in that in step
In rapid 2, according to required precision defined analysis grid, scatterplot number is inputted by 0-180 ° of orientation of certain step size computation in analysis grid
According to directional derivative, record minimum directional derivative and maximum directional derivative size, orientation.
4. the change face element plane interpolation method according to claim 3 for moving towards analysis based on geology, which is characterized in that in step
In rapid 2, the determination of geology trend is completed in accordance with the following methods:
If a=f (x, y) is defined in certain field of ground point P (x, y), direct north azimuth is 0 °；Injection line l, and just
Northern angular separation isThen P point is defined as along the directional derivative of ray l
For directional derivative expression geological property a from P point along the change rate in the direction ray l, minimum directional derivative corresponds to azimuthIt represents
Geology trend.
5. the change face element plane interpolation method according to claim 4 for moving towards analysis based on geology, which is characterized in that in step
In rapid 3, orientation is corresponded to according to the minimum directional derivative being calculated in step 2 and determines transverse direction；According to long axis direction
Directional derivative and short-axis direction directional derivative size relation, determine the ratio between transverse and short axle；According to two along the direction ray l
The nearest extreme point of subderivative determines the size of long axis.
6. the change face element plane interpolation method according to claim 1 for moving towards analysis based on geology, which is characterized in that in step
In rapid 4, using oval interpolation face element as analysis face element, the method for defining ellipsometric analysis face element:
Oval direction definition is that minimum directional derivative corresponds to azimuth
Oval eccentricity is acquired using long axis direction directional derivative corresponding with short-axis direction:
If long axis and short-axis direction directional derivative are respectively a, b, enableThen oval eccentricity e=c/a；
For the plane geological property scatterplot data inputted in step 1,0-180 ° of certain step size computation is pressed in analysis grid
Orientation inputs the directional derivative of scatterplot data, acquires minimum directional derivative；The grid data edge constituted for minimum directional derivative
Geology is moved towards to carry out secondary derivation, determines first extreme point P ˊ (x ˊ, y ˊ) of second derivative in this direction, point of origin P (x, y) with
The distance between extreme point is defined as transverse length a.
7. the change face element plane interpolation method according to claim 6 for moving towards analysis based on geology, which is characterized in that in step
In rapid 4, for analyzing the arbitrary point in grid, direction, eccentricity and the long axis of ellipsometric analysis face element used in interpolation calculation
Size be it is variable, i.e., utmostly guarantee participate in interpolation calculation known scatterplot between have preferable geology correlation.
8. the change face element plane interpolation method according to claim 1 for moving towards analysis based on geology, which is characterized in that in step
In rapid 5, for the arbitrary point on analysis grid defined in step 2, using ellipse conformal analysis face element obtained in step 3 as about
Beam carries out planar interpolation, is sketched isopleth using interpolation result, and carry out geological verification.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910618381.2A CN110390702B (en) | 2019-07-09 | 2019-07-09 | Variable surface element plane interpolation method based on geological trend analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910618381.2A CN110390702B (en) | 2019-07-09 | 2019-07-09 | Variable surface element plane interpolation method based on geological trend analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110390702A true CN110390702A (en) | 2019-10-29 |
| CN110390702B CN110390702B (en) | 2022-11-08 |
Family
ID=68286423
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910618381.2A Active CN110390702B (en) | 2019-07-09 | 2019-07-09 | Variable surface element plane interpolation method based on geological trend analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110390702B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111242857A (en) * | 2020-01-06 | 2020-06-05 | 中国石油化工股份有限公司 | Contour line generation optimization method with geological direction characteristics |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030132934A1 (en) * | 2001-12-12 | 2003-07-17 | Technoguide As | Three dimensional geological model construction |
| CN103353923A (en) * | 2013-06-26 | 2013-10-16 | 中山大学 | Self-adaption spatial interpolation method and system based on spatial feature analysis |
| CN107211592B (en) * | 2011-02-18 | 2014-11-05 | 核工业北京地质研究院 | A kind of geological information interpolation computing method based on Spatial variability |
| CN104678434A (en) * | 2013-12-02 | 2015-06-03 | 中国石油化工股份有限公司 | Method for predicting storage layer crack development parameters |
-
2019
- 2019-07-09 CN CN201910618381.2A patent/CN110390702B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030132934A1 (en) * | 2001-12-12 | 2003-07-17 | Technoguide As | Three dimensional geological model construction |
| CN107211592B (en) * | 2011-02-18 | 2014-11-05 | 核工业北京地质研究院 | A kind of geological information interpolation computing method based on Spatial variability |
| CN103353923A (en) * | 2013-06-26 | 2013-10-16 | 中山大学 | Self-adaption spatial interpolation method and system based on spatial feature analysis |
| CN104678434A (en) * | 2013-12-02 | 2015-06-03 | 中国石油化工股份有限公司 | Method for predicting storage layer crack development parameters |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111242857A (en) * | 2020-01-06 | 2020-06-05 | 中国石油化工股份有限公司 | Contour line generation optimization method with geological direction characteristics |
| CN111242857B (en) * | 2020-01-06 | 2023-03-24 | 中国石油化工股份有限公司 | Contour line generation optimization method with geological direction characteristics |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110390702B (en) | 2022-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Roy et al. | Generative topographic mapping for seismic facies estimation of a carbonate wash, Veracruz Basin, southern Mexico | |
| CN110031896B (en) | Seismic random inversion method and device based on multi-point geostatistics prior information | |
| US8009170B2 (en) | Systems and methods for computing a variogram model | |
| CN105425292A (en) | Oil and gas prediction method and oil and gas prediction device | |
| CN104297785A (en) | Lithofacies constrained reservoir physical property parameter inversion method and device | |
| WO2012135043A1 (en) | Detecting structural and stratigraphic information from seismic data | |
| US10884150B2 (en) | Systems and methods for using probabilities of lithologies in an inversion | |
| Deng et al. | Fast Bayesian inversion method for the generalized petrophysical and compositional interpretation of multiple well logs with uncertainty quantification | |
| CN103869362A (en) | Method and equipment for obtaining body curvature | |
| CN102866426A (en) | Method for analyzing oil gas information of rock body by applying amplitude versus offset (AVO) large-angle trace gathers | |
| CN107656312B (en) | River channel sand body prediction method and device based on sub-azimuth superposition | |
| Lin et al. | Noise filtering, trajectory compression and trajectory segmentation on GPS data | |
| CN110390702A (en) | The change face element plane interpolation method of analysis is moved towards based on geology | |
| CN104122584A (en) | Method and device for determining directionality according to seismic data | |
| CN105093300B (en) | A kind of boundary recognition of geological body method and device | |
| CN105277974A (en) | Stratum data interpolation method | |
| WO2021127382A1 (en) | Full waveform inversion in the midpoint-offset domain | |
| CN102830430B (en) | A kind of horizon velocity modeling method | |
| CN112684498A (en) | Reservoir fracture prediction method and system based on wide-azimuth seismic data | |
| US20090043549A1 (en) | Methods, apparatus, and products for seismic ray tracing | |
| US10082588B2 (en) | Adaptive structure-oriented operator | |
| CN109581496B (en) | Bidirectional construction dip angle constraint method and device in tomography inversion | |
| Gervais et al. | Identifying influence areas with connectivity analysis–application to the local perturbation of heterogeneity distribution for history matching | |
| CN111381279A (en) | Reservoir porosity quantitative prediction method and device | |
| CN105095634A (en) | Migration velocity model building method based on geological mass |
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 |