CN109633781A - Geological property acquisition methods and device, electronic equipment, storage medium - Google Patents
Geological property acquisition methods and device, electronic equipment, storage medium Download PDFInfo
- Publication number
- CN109633781A CN109633781A CN201910004452.XA CN201910004452A CN109633781A CN 109633781 A CN109633781 A CN 109633781A CN 201910004452 A CN201910004452 A CN 201910004452A CN 109633781 A CN109633781 A CN 109633781A
- Authority
- CN
- China
- Prior art keywords
- sample
- geological
- point
- sample point
- distribution
- 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
- 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)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
This application discloses a kind of geological property acquisition methods and devices, the described method includes: according to well log interpretation curve, the attribute information of each sample point in corresponding geological sample space of logging well is obtained, the attribute information of the sample point includes the coordinate value and geological property value where sample point in geology sample space；It is handled by the attribute information to each sample point, obtains the property distribution rules of the geological sample space towards sample point；According to the property distribution rules of the geological sample space towards sample point, the geological property value of each unknown point in the geological sample space is calculated；It is distributed the geological property value of the sample point and the unknown point as the geological property in the geological sample space.According to method provided herein, the geological property regularity of distribution of geological sample can be accurately obtained, provides important evidence for the geological prospecting to be carried out.
Description
Technical field
This application involves computer application technology, in particular to a kind of geological property acquisition methods and device, electronics
Equipment, computer readable storage medium.
Background technique
In geological prospecting, it usually needs sampled to the geologic province explored into well logging, obtain geological sample
Log data obtains the geological property distribution of geological sample, thus by ground then by being analyzed log data and being handled
The geological property distribution of matter sample is generalized in the exploration and research of entire geologic province.
In existing realization, the general well logging for carrying out geologic province by logging equipment is sampled, and obtains the survey of geological sample
Then well data carry out well log interpretation to geological sample according to log data, and generate well log interpretation curve according to explanation results,
To which the exploration for carrying out corresponding geologic province by well log interpretation curve is studied.But well log interpretation curve can only reflect sampling
The corresponding geological property of point can not be according to well log interpretation song for not carrying out the position of log data acquisition in geological sample
Line knows corresponding geological property, can not only reflect the geological property distribution of geological sample comprehensively by well log interpretation curve,
To influence the levels of precision of geological prospecting.
Therefore, the geological property distribution for how accurately acquiring geological sample, is currently still to have technical problem to be solved.
It should be noted that information is only used for reinforcing the reason to the background of the application disclosed in above-mentioned background technology part
Solution, therefore may include the information not constituted to the prior art known to persons of ordinary skill in the art.
Summary of the invention
Based on above-mentioned technical problem, this application provides a kind of geological property acquisition methods and device, electronic equipment, calculating
Machine readable storage medium storing program for executing.
Wherein, technical solution used by the application are as follows:
A kind of geological property acquisition methods, comprising: according to well log interpretation curve, obtain corresponding geological sample space of logging well
In each sample point attribute information, the attribute information of the sample point includes the coordinate value in geology sample space where sample point
With geological property value；It is handled by the attribute information to each sample point, obtains the geological sample space towards sample point
Property distribution rules；According to the property distribution rules of the geological sample space towards sample point, the geological sample is calculated
The geological property value of each unknown point in space；Using the geological property value of the sample point and the unknown point as the geology sample
The geological property in this space is distributed.
A kind of geological property acquisition device, comprising: sample point attribute information obtains module, for according to well log interpretation song
Line obtains the attribute information of each sample point in corresponding geological sample space of logging well, and the attribute information of the sample point includes sample
Coordinate value and geological property value where this point in geology sample space；Sample point property distribution rules obtain module, for leading to
It crosses and the attribute information of each sample point is handled, obtain the property distribution rules of the geological sample space towards sample point；
Unknown point geological property obtains module, for calculating according to the property distribution rules of the geological sample space towards sample point
The geological property value of each unknown point in the geological sample space；Geological property distributed acquisition module is used for the sample point
It is distributed with the geological property value of the unknown point as the geological property in the geological sample space.
A kind of electronic equipment, including processor and memory, the memory is for storing holding for the processor
Row instruction；Wherein, the processor is configured to execute above-mentioned geological property acquisition methods via the executable instruction is executed.
A kind of computer readable storage medium, is stored thereon with computer program, and the computer program is held by processor
Above-mentioned geological property acquisition methods are realized when row.
In the above-mentioned technical solutions, according to well log interpretation curve, each sample in corresponding geological sample space of logging well is obtained
The attribute information of point, and by the processing to each sample point attribute information, obtain attribute of the geological sample space towards sample point
The regularity of distribution, then the property distribution rules according to geology sample space towards each sample point, calculate each in geological sample space
The geological property value of unknown point, using the geological property value of sample point and unknown point as the geological property in the geological sample space
Distribution.
Technical solution provided herein has fully considered property distribution rules of the geological sample space towards sample point,
And construct and meet the geological property acquisition model containing a large amount of space sample points, so as to carry out ground according to constructed model
The geological property value of unknown point calculates in matter sample space, passes through the ground of sample point and unknown point in geological sample space as a result,
Matter attribute come accurately reflect geological sample geological property distribution, for the geological prospecting to be carried out provide important exploration according to
According to.
It should be understood that above general description and following detailed description be only it is exemplary and explanatory, not
The application can be limited.
Detailed description of the invention
The drawings herein are incorporated into the specification and forms part of this specification, and shows the implementation for meeting the application
Example, and in specification it is used to explain the principle of the application together.
Fig. 1 is a kind of flow chart of geological property acquisition methods shown according to an exemplary embodiment；
Fig. 2 is according to the flow chart that step 120 is described shown in Fig. 1 corresponding embodiment；
Fig. 3 is according to the flow chart that step 130 is described shown in Fig. 1 corresponding embodiment；
Fig. 4 is a kind of flow chart of the geological property acquisition methods shown according to another exemplary embodiment；
Fig. 5 is a kind of block diagram of geological property acquisition device shown according to an exemplary embodiment；
Fig. 6 is the hardware block diagram of a kind of electronic equipment shown according to an exemplary embodiment.
Through the above attached drawings, it has been shown that the specific embodiment of the application will be hereinafter described in more detail, these attached drawings
It is not intended to limit the range of the application design in any manner with verbal description, but is by referring to specific embodiments
Those skilled in the art illustrate the concept of the application.
Specific embodiment
Here will the description is performed on the exemplary embodiment in detail, the example is illustrated in the accompanying drawings.Following description is related to
When attached drawing, unless otherwise indicated, the same numbers in different drawings indicate the same or similar elements.Following exemplary embodiment
Described in embodiment do not represent all embodiments consistent with the application.On the contrary, they be only with it is such as appended
The example of the consistent device and method of some aspects be described in detail in claims, the application.
Fig. 1 is a kind of flow chart of geological property acquisition methods according to shown by an exemplary embodiment.Such as Fig. 1 institute
Show, this method at least may comprise steps of:
Step 110, according to well log interpretation curve, the attribute letter of each sample point in corresponding geological sample space of logging well is obtained
Breath.
Firstly the need of explanation, well log interpretation curve is according to the log data of obtained geological sample, to geology sample
This progress well log interpretation obtains, and the log data of geological sample is then that the well logging for carrying out geologic province by logging equipment is adopted
Sample is obtained.
It can be regarded as by the well logging sampling that logging equipment carries out geologic province: special logging equipment used, along well bore
The change in physical of stratigraphic section and pit shaft is measured, to obtain the corresponding resistivity of each measurement point, natural potential, sound wave speed
The physical parameters such as degree.Well site where well logging is known as geological sample, each measurement point is known as sample point, each measurement point is surveyed
Amount gained physical parameter is known as the log data of geological sample.
Illustratively, depth match, curve smoothing processing, environmental correction, number are carried out in the log data to geological sample
It, can be by carrying out porosity calculation, water saturation and hydrocarbon saturation to log data after the pretreatments such as value standardization
The modes such as calculating, formation absolute permeability calculating, carry out qualitative interpretation to log data, to obtain rock corresponding to each sample point
The geological properties value such as property, oiliness, permeability.According to qualitative interpretation geological property value obtained, and according to each sample point
The position (such as depth-logger) of place geology sample space is sequentially generated well log interpretation achievement curve, by generating well logging solution
Release the variation of geological property in curve reflection geological sample space.
The attribute information of each sample point includes the coordinate where each sample point in geology sample space in geological sample space
Value and geological property value, the specific location of geology sample space where the coordinate value of sample point is used to indicate the sample point.Example
Property, the coordinate value of sample point is made of depth-logger and azimuth, deviation angle, the corresponding depth-logger of each sample point, side
Parallactic angle, deviation angle are recorded in log data, and after carrying out qualitative interpretation to log data, by well log interpretation curve, institute is anti-
The well logging trajectory coordinates reflected, which directly convert, to be obtained.Each sample point in well log interpretation curve record geological sample space as a result,
Attribute information.
It should be noted that for the logging operation of carried out geological sample, to the resulting log data of well logging into
After row qualitative interpretation, well logging result can be generated according to explanation results and explains curve, therefore, the described well logging solution of the present embodiment
Release curve be it is known, geological property acquisition methods described in the present embodiment are realized based on known well log interpretation curve
's.
Step 120, it is handled, is obtained in geological sample space towards sample point by the attribute information to each sample point
Property distribution rules.
Wherein, to the processing of each sample point attribute information, including the processing to each sample point coordinate value and to each sample point
The processing of geological property value；Correspondingly, the property distribution rules in geological sample space towards sample point include each sample point institute
The address properties regularity of distribution that the coordinate attributes regularity of distribution of reflection and each sample point are reflected.
As shown in Fig. 2, in one exemplary embodiment, by the processing to each sampling point coordinate value, to obtain geological sample
The process of the coordinate attributes regularity of distribution in space towards sample point at least may comprise steps of:
Step 121, it according to the coordinate value of each sample point, calculates separately in geological sample space between any two sample point
Distance.
As previously mentioned, can get the coordinate value and geology category of each sample point in geological sample space according to well log interpretation curve
Property value, if setting the coordinate value of sample point as (Xi,Yi,Zi), wherein " Xi" indicate the corresponding depth-logger of sample point, " Yi" indicate sample
The corresponding azimuth of this point, " Zi" indicate the corresponding deviation angle of sample point, and the geological property value of sample point is set as " Pi", and
The quantity of the included sample point in geological sample space is " N " (N is natural number), then obtains the sample point set in geological sample space
It closes:
{(X1,Y1,Z1,P1), (X2,Y2,Z2,P2) ... ..., (XN,YN,ZN,PN)}
Illustratively, after the sample point set for obtaining geological sample space, also to the exceptional sample in sample point set
Point is cleaned, and the exceptional sample point that sample point overall trend is not met in sample point set is excluded, to be promoted over the ground
The accuracy that matter sample space is obtained towards the property distribution rules of sample point.The exceptional sample point cleaning carried out is optional
The methods of exceptional sample point screening is carried out with 3 δ principles, cluster or with the average geology attribute value of population sample point, this place is not
It is limited.
It should be noted that the position where being intended to indicate that sample point due to the coordinate value of sample point in geology sample space
It sets, the sample point overall trend in sample point set is embodied by the geological property value of sample point, therefore, above-mentioned to sample point set
The execution that middle exceptional sample point is cleaned depends on the geological property value of each sample point.
After carrying out the cleaning of exceptional sample point, according to the coordinate value of retained each sample point, calculate in geological sample space
The distance between any two sample point, to obtain the distance between sample point two-by-two.For example, in each sample point coordinate value known
Under conditions of, the distance between any two sample point is calculated according to the range formula of point-to-point transmission in space.
Illustratively, it can also be obtained according to the distance between any two sample point in gained geological sample space is calculated
Obtain the distance sample space of sample point
Step 122, it by the distribution statistics to distance between any two sample point, obtains between any two sample point
The distribution characteristics of distance.
Wherein, distribution statistics are carried out to distance between any two sample point in geology sample space, including respectively to upper
State the process that each distance is grouped and counts in the distance sample space of sample point.
Maximum distance and minimum range can be obtained respectively from the distance sample space of above-mentioned sample point, and calculate it is maximum away from
From the difference between minimum range, to calculate the grouping that each distance is grouped in sample space of adjusting the distance according to gained difference
Step-length.
Illustratively, if setting in the distance sample space of sample point maximum distance as " hmax", minimum range is " hmin", and
If being grouped step-length is " Δ h ", then according to formula Δ h=(hmax-hmin)/m, which is calculated, obtains grouping step-length.Wherein, set by " m " expression
The fixed number of packet that distance between any two points is grouped, and m > 0.
According to grouping step-length obtained, the distance between any two sample point in geology sample space is divided
Group will meet the same distance apart from section and be divided into same group.
As a result, by counting the quantity of distance between divided any two sample point in each grouping respectively, and count
The ratio of the sum of distance shared by the quantity of distance in each grouping obtains the spacing of any two sample point in geological sample space
From the distribution character met, to be fitted to obtain according to the distribution character met and be with distance between any two sample point
The probability statistical distribution function " g (h) " of sample space, and ∑ g (hi)=1.
Probability statistical distribution function " g (h) " obtained, the spacing of any two sample point as in geological sample space
From distribution characteristics, which reflects the coordinate attributes regularity of distribution of the geological sample space towards sample point.
Illustratively, the quantity of divided distance and each grouping of statistics are divided apart from institute in counting each grouping
After accounting for the ratio of distance sum, if obtaining the distance between any two sample point in the geological sample space meets normal state
Distribution then is fitted to obtain Density Function of Normal Distribution to be above-mentioned probability statistical distribution function " g (h) " based on least square method.
It should be noted that based on the distribution that distance is met between any two sample point in different geology sample spaces
Characteristic may also be different, and the concrete functional form for being fitted to obtain probability statistical distribution function " g (h) " may also be different.Citing
For, fitting gained probability statistical distribution function " g (h) " can also be cumulative distribution density function, t- distribution density function, Thailand
Forms, this place such as series expansion density function is strangled to be not limited.
In a further exemplary embodiment, it by the processing to sample point geological property value each in geology sample space, obtains
The process for obtaining the geological property regularity of distribution towards sample point in geological sample space, sits each sample point with above embodiment
The treatment process of scale value is similar, but not fully identical.
Illustratively, in the sample point set obtained after carrying out the cleaning of exceptional sample point, each sample point is obtained respectively
Geological property value, and to the geological property value of each sample point carry out distribution statistics, to obtain in geological sample space towards sample
The geological property regularity of distribution of this point.It includes the geology to sample point respectively that the geological property value of each sample point, which carries out distribution statistics,
The process that attribute value is grouped and counts.
By obtaining maximum geology attribute value and most respectively from the geological property value of sample point each in geological sample space
Small geological property value, and the difference between maximum geology attribute value and minimum geology attribute value is calculated, in terms of according to gained difference
Calculate the grouping step-length being grouped to the geological property value of sample point.By according to grouping step-length obtained, to geological sample
The geological property value of sample point is grouped in space, and the geological property value for meeting same geological property section is divided into together
One group.
As a result, by counting the quantity of divided geological property value in each grouping respectively, and with counting in each grouping
The ratio of geological property value sum shared by the quantity of matter attribute value obtains the geological property value institute of sample point in geological sample space
The distribution character of satisfaction, to be fitted to obtain the geological property value with sample point according to the distribution character that geology attribute value is met
For the probability statistical distribution function " f (p) " of sample space, and ∑ f (pi)=1.
Probability statistical distribution function " f (p) " obtained is each sample point geological property value in geological sample space
Distribution characteristics, the distribution characteristics reflect the coordinate attributes regularity of distribution of the geological sample space towards sample point.
It should be noted that the distribution that the geological property value based on each sample point in different geology sample spaces is met is special
Property may also be different, be fitted to obtain probability statistical distribution function " f (p) " concrete functional form may also difference, for example, intending
Close gained probability statistical distribution function " f (p) " can be Density Function of Normal Distribution, cumulative distribution density function, t- distribution it is close
Function, any one in taylor series expansion density function are spent, this place is not limited.
Step 130, it according to the property distribution rules of sample point each in geology sample space, calculates each in geological sample space
The geological property value of unknown point.
Wherein, each unknown point in geological sample space refers to, the difference of each sample point is different from geological sample space
Position.
Illustratively, according to the logging operation carried out, the spatial model in automation building geological sample space, by right
Gained spatial model carries out gridding processing, using the central point of each grid as the spatial point in geological sample space, and will
Spatial point in addition to sample point is known as unknown point.Geological sample space can be regarded as the sky as corresponding to each sample point as a result,
Between area of space corresponding to region and each unknown point constitute.
Due in the spatial model in constructed geological sample space the relative position of sample point and unknown point be it is fixed, can
To be reference according to the coordinate value of each sample point, the coordinate value of each unknown point is obtained.
As shown in figure 3, in one exemplary embodiment, the realization of above-mentioned steps 130 may comprise steps of:
Step 131, to each unknown point for being different from sample point, if the random screening dry sample from geological sample space respectively
This point solves sample space as the distance of unknown point.
Firstly the need of explanation, method described by the present embodiment is respectively for each unknown in geological sample space
What point proposed, in other words, to each unknown point in geology sample space, the content that can be described by executing the present embodiment,
The geological property value of each unknown point is obtained respectively.
Wherein, several sample points of random screening refer to from geological sample space, are distributed in not from geological sample space
Know and randomly select several sample points in the sample point of a surrounding, and makes the coordinate value of selected sample point with respect to unknown point coordinate
The sum of weight of value is equal to 1 or levels off to 1.Sample is solved using selected several sample points as the distance of unknown point as a result,
This space, also, according to the acquisition methods of sample point described above, can to each unknown point in geology sample space
The distance obtained corresponding thereto solves sample space.
It should be noted that distance solve sample space in, the coordinate value and geological property value of each sample point be by
Obtained by well log interpretation curve according to content described by step 110, in other words, distance solves each in sample space
The attribute value of sample point is distributed in the sample point set in geological sample space described by step 121.
Step 132, unknown point and the distance between arbitrary sample point in solution sample space are calculated.
As previously described, distance corresponding to the coordinate value of each unknown point and each unknown point in geological sample space
The coordinate value for solving arbitrary sample point in sample space is given value, therefore, can be asked according to the distance of point-to-point transmission in space
It solves formula and calculates the distance between arbitrary sample point in unknown point and respective distances solution sample space.
Illustratively, gained unknown point and the distance between arbitrary sample point expression in solution sample space will be calculated
For " h ' ", carried out unknown point and apart from the acquisition for solving the distribution characteristics of distance between arbitrary sample point in sample space below
Process is then to seek uniting using unknown point and apart from distance is solved in sample space between arbitrary sample point as the probability of sample space
Count the process of distribution function " g (h ') ".
Step 133, according to acquisition process corresponding to the distribution characteristics of distance between any two sample point, unknown point is obtained
The distribution characteristics of distance between arbitrary sample point in sample space is solved with distance.
Wherein, acquisition process corresponding to the distribution characteristics of distance is between any two sample point in geological sample space
Refer to, the described process that distribution statistics are carried out to the distance between any two sample point in geology sample space of step 122.
In other words, current embodiment require that by appointing to unknown point in geology sample space in solution sample space
Distance carries out distribution statistics between meaning sample point, to obtain the spacing of unknown point with the arbitrary sample point in solution sample space
From distribution characteristics.It below will be according to content described by step 122, to unknown point in geology sample space and apart from solution sample
The distribution characteristics acquisition process of distance is described in detail between arbitrary sample point in space:
In one exemplary embodiment, to unknown point and apart from solving in sample space, distance is carried out between arbitrary sample point
Distribution statistics still include that distance is grouped and counts between arbitrary sample point to unknown point and apart from solving in sample space
Process.
Wherein, it is walked to unknown point and apart from the grouping that distance is grouped between arbitrary sample point in sample space is solved
Grouping step-length " Δ h " used in long still optional step 122, and according to grouping step-length " Δ h " to unknown point and apart from solution
The distance between arbitrary sample point is grouped in sample space, will meet the same distance apart from section be divided into it is same
The quantity of group, gained grouping is consistent with number of packet " m " set in step 122.
By counting the quantity of divided distance in each grouping respectively, and count the number of divided distance in each grouping
The ratio for measuring shared distance sum, obtains unknown point and distance is met between the arbitrary sample point in solution sample space
Distribution character, to be fitted to obtain with unknown point and the arbitrary sample in solution sample space according to the distribution character met
Distance is the probability statistical distribution function " g (h ') " of sample space between point.Probability statistical distribution function " g (h ') " is geology
Unknown point and apart from the distribution characteristics for solving distance between arbitrary sample point in sample space in sample space.
In a further exemplary embodiment, since above-mentioned probability statistical distribution function " g (h ') " is according to probability statistics point
The acquisition process of cloth function " g (h) " is obtained, and it is empty that the difference between the two function is only embodied in the corresponding distance sample of fitting
Between difference, as a result, in geological sample space is calculated unknown point with apart from solve sample space in arbitrary sample point it
Between distance " h ' " after, will directly calculate resulting distance sample space and bring into probability statistical distribution function " g (h) " and count
It calculates, can accordingly obtain probability statistical distribution function " g (h ') ".
Step 134, by unknown point and apart from the distribution characteristics for solving distance between each sample point in sample space, as institute
The weight that unknown point relative distance solves each sample point geological property value in sample space is stated, by calculating each sample point geology category
Property value weight and, obtain the geological property value of unknown point.
Illustratively, if setting the coordinate value of any one unknown point in geological sample space as (X ', Y ', Z '), the unknown point
The calculation formula of corresponding geological property value is as follows:
Wherein, " λi" indicate that unknown point relative distance solves the weight of each sample point geological property value in sample space,
Substantially unknown point and apart from the distribution characteristics for solving distance between each sample point in sample space, that is, have λi=g (h '), " N ' "
Indicate that distance solves the quantity of sample point in sample space, " Pi" then indicate that distance solves the geology of each sample point in sample space
Attribute value.
Therefore, the coordinate category in present exemplary embodiment institute providing method, according to geology sample space towards sample point
Property the regularity of distribution, and combine geology sample space in each sample point geological property value, with calculate separately to obtain geological sample sky
Between in each unknown point geological property value.
Step 140, it is distributed the geological property value of sample point and unknown point as the geological property in geological sample space.
Wherein, the sky as corresponding to geological sample space area of space as corresponding to each sample point and each unknown point
Between region constitute, therefore, can be empty as geological sample using the geological property value of sample point each in geological sample space and unknown point
Between geological property distribution, carry out comprehensively and accurate embody so as to the geological property distribution to geology sample space.
In one exemplary embodiment, by can by the geological property value of sample point each in geological sample space and unknown point
It is shown in the spatial model in constructed geological sample space depending on changing ground, that is, may make that user is clear according to visual spatial model
Obtain to Chu the geological property distribution in geological sample space.
To sum up, in method provided herein, the above-mentioned calculating process to unknown point geological property value can be considered as institute
Building meets the geological property containing a large amount of space sample points and obtains model, is based on the computation model, can be according to known survey
Well explains that the geological property value of several unknown points is calculated in curve, to pass through sample point and unknown point in geological sample space
Geological property come accurately reflect geological sample geological property distribution, provide important survey for the geological prospecting to be carried out
Visit foundation.
Further, since the geological property value of above-mentioned unknown point is according to the coordinate category in geology sample space towards sample point
Property the regularity of distribution be calculated, enable unknown point obtained by the application geological property value and unknown point in geological sample sky
Between in true geological property distribution be adapted, the accuracy of geological property value corresponding to collateral security unknown point.
In a further exemplary embodiment, as shown in figure 4, after step 120, above-mentioned geological property acquisition methods also wrap
Include following steps:
Step 210, the geological property value according to each unknown point in geological sample space, statistics obtain geological sample space
The geological property regularity of distribution towards unknown point.
Illustratively, after the geological property value for obtaining each unknown point in geological sample space based on the above embodiment, according to
The coordinate value of each unknown point obtains the unknown point set in geological sample space:
{(X1′,Y1′,Z1′,P1'), (X2′,Y2′,Z2′,P2') ... ..., (XM′YM′,ZM′,PM′)}
Wherein, " M " indicates the quantity of unknown point in geological sample space, can be with sample point in geological sample space
Quantity is identical, can also be different, this place is not limited.
According to the acquisition of the geological property regularity of distribution in geological sample space described by above-mentioned steps 122 towards sample point
Process, by the way that geological property value of each unknown point in geological sample space is grouped and is counted, to obtain geology sample
The geological property regularity of distribution of this space towards unknown point.
Wherein, the grouping step-length that the selected geological property value to each unknown point in geology sample space is grouped, with
Used grouping step-length is consistent when obtaining the geological property regularity of distribution in geological sample space towards sample point, and according to
The grouping step-length is grouped the geological property value of each unknown point, and the number of packet of obtained unknown point is also and in step 122
The number of packet of divided sample point geological property value is consistent.
As a result, by being counted in the quantity of divided geological property value in each grouping, and each grouping of statistics respectively
The ratio of geological property value sum shared by the quantity of divided geological property value, and the geology category with unknown point is obtained by fitting
Property value be sample space probability statistical distribution function " F (p ') ".Probability statistical distribution function " F (p ') " is geological sample sky
Between the geological property regularity of distribution towards unknown point.
Illustratively, can also directly by the geological property value of each unknown point bring into probability statistical distribution function " f (p) " into
Row calculates, and can accordingly obtain probability statistical distribution function " F (p ') ".
Step 220, if the geological property regularity of distribution of the geological sample space towards unknown point with towards the ground of sample point
Deviation between matter property distribution rules is higher than preset value, updates to distance between any two sample point in geology sample space
The grouping step-length being grouped.
Wherein, the geological property regularity of distribution of the geological sample space towards unknown point and the geological property towards sample point point
Deviation between cloth rule, can pass through optimization objective function " min ∑ (F (p)-f (p))2" indicate, when calculated result is lower than pre-
If when value, indicating the coordinate attributes regularity of distribution and each sample point institute of the obtained geological sample space of step 122 towards sample point
It is adapted between the true distribution in geological sample space, the geological property value of each unknown point in obtained geological sample space
It is adapted with true geological property value, at this moment can ensure the essence of the geological property regularity of distribution corresponding to geological sample space
True property.
As a result, if gained deviation is higher than preset value, the coordinate attributes to geology sample space towards sample point are needed
The acquisition process of the regularity of distribution optimizes.By being grouped to distance between any two sample point in geology sample space
Grouping step-length " Δ h " optimize, and then realize optimization to the distribution statistics process of distance between any two sample, directly
True point of geology sample space where the coordinate attributes regularity of distribution and each sample point to geological sample space towards sample point
It is adapted between cloth, to obtain the coordinate attributes regularity of distribution of geology sample space where capable of accurately reflecting each sample point.
Illustratively, calculating is iterated to grouping step-length " Δ h " using gradient descent method, obtains optimized grouping step
It is long.It, can be by the grouping step-length to geological sample after every grouping step-length for carrying out a step gradient descent algorithm and obtaining optimization
The grouping step-length " Δ h " that distance is grouped between any two sample point in space is updated.
Calculation formula corresponding to steepest descent method indicates are as follows: Δ hk+1=Δ hk+akdk, wherein " ak" indicate step-length because
Son, the step factor are less than 1 setting value, and " k " indicates the step number of decline, " dk" indicate to decline the gradient of " k " step, the ladder
Degree includes positive direction gradient and negative direction gradient.
Step 230, geological property of the geological sample space towards unknown point is reacquired according to the grouping step-length updated
The regularity of distribution, until between the geological property regularity of distribution towards unknown point and the geological property regularity of distribution towards sample point
Deviation is lower than preset value.
Wherein, the grouping step-length " Δ h " that distance is grouped between any two sample point in geology sample space
After being updated, it can be re-execute the steps 121, step 122 and step 130 according to the grouping step-length of update, and according to institute
The geological property value of each unknown point is obtained, the geological property regularity of distribution of the geological sample space towards unknown point is calculated.
It executes repeatedly, up to the geological property regularity of distribution of the geological sample space towards unknown point and towards sample point
The geological property regularity of distribution between deviation be lower than preset value.
Fig. 5 is a kind of geological property acquisition device according to shown by an exemplary embodiment.As shown in figure 5, the device
Include:
Sample point attribute information obtains module 310, for obtaining corresponding geological sample of logging well according to well log interpretation curve
The attribute information of each sample point in space, the attribute information of the sample point include the coordinate where sample point in geology sample space
Value and geological property value；
Sample point property distribution rules obtain module 320 and obtain for being handled by the attribute information to each sample point
Obtain property distribution rules of the geological sample space towards sample point；
Unknown point geological property obtains module 330, for being advised according to geology sample space towards the property distribution of sample point
Rule calculates the geological property value of each unknown point in geological sample space；
Geological property distributed acquisition module 340, for using the geological property value of each sample point and each unknown point as geology
The geological property of sample space is distributed.
In a further exemplary embodiment, sample point property distribution rules acquisition module 320 includes:
Sample point distance computation unit calculates separately in geological sample space and appoints for the coordinate value according to each sample point
It anticipates the distance between two sample points；
Sample point spacing distribution characteristics acquiring unit, for by between any two sample point in geology sample space
The distribution statistics of distance obtain the distribution characteristics of distance between any two sample point in geological sample space.Geological sample is empty
Between between middle any two sample point distance distribution characteristics for reflecting coordinate attributes of the geological sample space towards sample point
The regularity of distribution.
In a further exemplary embodiment, above-mentioned sample point spacing distribution characteristics acquiring unit includes:
Sample point end value spacing obtains subelement, in the distance in geological sample space between any two sample point
In, maximum distance and minimum range are obtained respectively；
It is grouped step size computation subelement, for the difference according to acquired maximum distance and the minimum range, calculating pair
The grouping step-length that distance is grouped between any two sample point in geological sample space；
Sample point spacing is grouped subelement, is used for according to grouping step-length obtained to any two in geology sample space
The distance between sample point is grouped；
Sample point spacing counts subelement, for by statistics geology sample space between any two sample point away from
From the distribution situation in each grouping, the distribution characteristics of distance between any two sample point in geological sample space is obtained.
In a further exemplary embodiment, unknown point geological property acquisition module 330 includes:
Distance solves sample space acquiring unit, for being different from each unknown of sample point in geological sample space
Point, respectively from geological sample space several sample points of random screening as unknown point distance solve sample space；
Unknown point distance computation unit, for calculating in geological sample space unknown point and appointing apart from solving in sample space
The distance between meaning sample point；
Unknown point spacing distribution characteristics acquiring unit, for the spacing according to any two sample point in geological sample space
From distribution characteristics corresponding to acquisition process, obtain geological sample space in unknown point with apart from solve sample space in any sample
The distribution characteristics of distance between this point；
Unknown point geological property value acquiring unit, for by unknown point in geological sample space and apart from solution sample space
In between each sample point distance distribution characteristics, solved as unknown point relative distance in geological sample space each in sample space
The weight of sample point geological property value, by calculate distance solve sample space in each sample point geological property value weight and,
Obtain the geological property value of unknown point in geological sample space.
In a further exemplary embodiment, it further includes sample point geology category that sample point property distribution rules, which obtain module 320,
Property distribution characteristics acquiring unit, for by geology sample space each sample point geological property value carry out distribution statistics,
Obtain the geological property regularity of distribution of the geological sample space towards sample point.
In a further exemplary embodiment, geological property acquisition device further include:
Unknown point geological property value obtains module, for the geological property according to each unknown point in geological sample space
Value, statistics obtain the geological property regularity of distribution of the geological sample space towards unknown point；
Be grouped step-length update module, for the geological property regularity of distribution of the geological sample space towards unknown point with towards
In the case that deviation between the geological property regularity of distribution of sample point is higher than preset value, update to any in geology sample space
The grouping step-length that distance is grouped between two sample points；
Be grouped step-length update execution module, for according to the grouping step-length updated regain geological sample space towards
The geological property regularity of distribution of unknown point, up to the geological property regularity of distribution of the geological sample space towards unknown point and towards sample
Deviation between the geological property regularity of distribution of this point is lower than preset value.
It should be noted that method provided by device provided by above-described embodiment and above-described embodiment belongs to same structure
Think, the concrete mode that wherein modules execute operation is described in detail in embodiment of the method, no longer superfluous herein
It states.
In a kind of exemplary embodiment, a kind of electronic equipment, including processor and memory, wherein processor is used for
Execute the above method.
Fig. 6 is the structural schematic diagram of a kind of electronic equipment according to shown by an exemplary embodiment.It should be noted that
The electronic equipment is one and adapts to example of the invention, must not believe that there is provided to any of use scope of the invention
Limitation.
The hardware configuration of the electronic equipment can generate biggish difference due to the difference of configuration or performance, as shown in fig. 6,
Electronic equipment includes: power supply 410, interface 430, at least a memory 450 and an at least central processing unit 470.
Wherein, power supply 410 is used to provide operating voltage for each hardware device on electronic equipment.
Interface 430 includes an at least wired or wireless network interface 431, at least a string and translation interface 433, at least one defeated
Enter output interface 435 and at least USB interface 437 etc., is used for and external device communication.
The carrier that memory 450 is stored as resource, can be read-only memory, random access memory, disk or CD
Deng the resource stored thereon includes operating system 451, application program 453 or data 455 etc., and storage mode can be short
Temporary storage permanently stores.Wherein, operating system 451 is for managing and each hardware device in controlling electronic devices and answering
With program 453, to realize calculating and processing of the central processing unit 470 to mass data 455.Application program 453 is based on operation
The computer program that at least one particular job is completed on system 451, may include that an at least module (is not shown in Fig. 6
Out), each module can separately include the series of computation machine readable instruction to electronic equipment.Data 455, which can be, deposits
The key message etc. being stored in disk.
Central processing unit 470 may include the processor of one or more or more, and be set as through bus and memory
450 communications, for the mass data 455 in operation and processing memory 450.
Described in detail above, be applicable in electronic equipment of the invention will be by depositing in the reading storage of central processing unit 470
The form of the series of computation machine readable instruction of storage come complete above-mentioned digital cores hole and skeleton recognition methods.
It should be noted that the electronic equipment, which is one, adapts to example of the invention, it must not believe that there is provided right
Any restrictions of use scope of the invention.
In a kind of exemplary embodiment, a kind of computer readable storage medium is stored thereon with computer program, the meter
Calculation machine program realizes the above method when being executed by processor.
In addition, above-mentioned attached drawing is only the schematic theory of the processing according to included by the method for the application exemplary embodiment
It is bright, rather than limit purpose.It can be readily appreciated that the time that above-mentioned processing shown in the drawings did not indicated or limited these processing is suitable
Sequence.In addition, be also easy to understand, these processing, which can be, for example either synchronously or asynchronously to be executed in multiple modules.
Those skilled in the art after considering the specification and implementing the invention disclosed here, will readily occur to its of the application
His embodiment.This application is intended to cover any variations, uses, or adaptations of the application, these modifications, purposes or
Adaptive change follow the general principle of the application and including the undocumented common knowledge in the art of the application or
Conventional techniques.The description and examples are only to be considered as illustrative, and the true scope and spirit of the application are by claim
It points out.
Claims (10)
1. a kind of geological property acquisition methods, which is characterized in that described method includes following steps:
According to well log interpretation curve, the attribute information of each sample point in corresponding geological sample space of logging well, the sample are obtained
The attribute information of point includes the coordinate value and geological property value where sample point in geology sample space；
It is handled by the attribute information to each sample point, obtains the property distribution of the geological sample space towards sample point
Rule；
According to the property distribution rules of the geological sample space towards sample point, calculate each unknown in the geological sample space
The geological property value of point；
It is distributed the geological property value of the sample point and the unknown point as the geological property in the geological sample space.
2. the method according to claim 1, wherein the property distribution rules towards sample point include described
The coordinate attributes regularity of distribution of sample point, the attribute information by each sample point are handled, and obtain the geology sample
Property distribution rules of this space towards sample point, comprising:
According to the coordinate value of each sample point, the distance between any two sample point in the geological sample space is calculated separately；
By the distribution statistics to distance between any two sample point, distance between any two sample point is obtained
Distribution characteristics, the distribution characteristics of distance is for reflecting the geological sample space towards sample between any two sample point
The coordinate attributes regularity of distribution of this point.
3. according to the method described in claim 2, it is characterized in that, described by distance between any two sample point
Distribution statistics, obtain the distribution characteristics of distance between any two sample point, comprising:
In distance between any two sample point, maximum distance and minimum range are obtained respectively；
According to the difference of the maximum distance and the minimum range, calculating carries out distance between any two sample point
The grouping step-length of grouping；
The distance between described any two sample point is grouped according to the grouping step-length；
By counting distribution situation of the distance between the described any two sample point in each grouping, described any two are obtained
The distribution characteristics of distance between a sample point.
4. according to the method described in claim 3, it is characterized in that, it is described according to the geological sample space towards sample point
Property distribution rules calculate geological property value corresponding to each unknown point in the geological sample space, comprising:
To each unknown point for being different from the sample point, several sample points of random screening from the geological sample space respectively
Distance as the unknown point solves sample space；
It calculates the unknown point and the distance solves the distance between arbitrary sample point in sample space；
According to acquisition process corresponding to the distribution characteristics of distance between any two sample point, the unknown point and institute are obtained
State the distribution characteristics that distance solves distance between arbitrary sample point in sample space；
By the unknown point and apart from the distribution characteristics for solving distance between each sample point in sample space, as the unknown point
The relatively described distance solves the weight of each sample point geological property value in sample space, by calculating each sample point geological property value
Weight and, obtain the geological property value of the unknown point.
5. according to the described in any item methods of claim 2 to 4, which is characterized in that the property distribution rule towards sample point
Rule further includes the geological property regularity of distribution of the sample point, and the attribute information by each sample point is handled, obtained
Obtain the property distribution rules of the geological sample space towards sample point, further includes:
Distribution statistics are carried out by the geological property value to each sample point, obtain the ground of the geological sample space towards sample point
Matter property distribution rules.
6. according to the method described in claim 5, it is characterized in that, described according to each sample point in the geological sample space
The geological property regularity of distribution, calculate in the geological sample space after geological property value corresponding to each unknown point, the side
Method further include:
According to geological property value of each unknown point in the geological sample space, statistics obtain the geological sample space towards
The geological property regularity of distribution of unknown point；
If the geological property regularity of distribution of the geological sample space towards unknown point and the geology category towards sample point
Property the regularity of distribution between deviation be higher than preset value, update to the grouping that distance is grouped between any two sample point
Step-length；
The geological property regularity of distribution of the geological sample space towards unknown point is regained according to the grouping step-length updated,
Described in until between the geological property regularity of distribution and the geological property regularity of distribution towards sample point towards unknown point
Deviation is lower than the preset value.
7. a kind of geological property acquisition device, which is characterized in that described device includes:
Sample point attribute information obtains module, for according to well log interpretation curve, obtaining in corresponding geological sample space of logging well
The attribute information of each sample point, the attribute information of the sample point include coordinate value in geology sample space where sample point and
Geological property value；
Sample point property distribution rules obtain module, for being handled by the attribute information to each sample point, described in acquisition
Property distribution rules of the geological sample space towards sample point；
Unknown point geological property obtains module, for according to the property distribution rules of the geological sample space towards sample point,
Calculate the geological property value of each unknown point in the geological sample space；
Geological property distributed acquisition module, for using the geological property value of the sample point and the unknown point as the geology
The geological property of sample space is distributed.
8. device according to claim 7, which is characterized in that the sample point property distribution rules obtain module and include:
Sample point distance computation unit calculates separately in the geological sample space and appoints for the coordinate value according to each sample point
It anticipates the distance between two sample points；
Sample point spacing distribution characteristics acquiring unit, for being united by the distribution to distance between any two sample point
Meter, obtains the distribution characteristics of distance between any two sample point.
9. a kind of electronic equipment characterized by comprising
Processor；And
Memory, for storing the executable instruction of the processor；
Wherein, the processor be configured to via execute the executable instruction come perform claim require it is 1 to 6 described in any item
Method.
10. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program
Method as claimed in any one of claims 1 to 6 is realized when being executed by processor.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2018109338204 | 2018-08-16 | ||
| CN201810933820 | 2018-08-16 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109633781A true CN109633781A (en) | 2019-04-16 |
| CN109633781B CN109633781B (en) | 2021-03-02 |
Family
ID=66056697
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910004452.XA Active CN109633781B (en) | 2018-08-16 | 2019-01-03 | Geological property acquisition method and device, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109633781B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112016048A (en) * | 2020-07-28 | 2020-12-01 | 中交广州航道局有限公司 | Geological attribute prediction method based on multi-point geological statistics |
| CN113361263A (en) * | 2021-06-04 | 2021-09-07 | 中国人民解放军战略支援部队信息工程大学 | Character entity attribute alignment method and system based on attribute value distribution |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100228484A1 (en) * | 2007-12-06 | 2010-09-09 | Exxon-Mobile, Upstream Research Company | Volume of Investigation Based Image Processing |
| US20160047934A1 (en) * | 2013-04-02 | 2016-02-18 | Schlumberger Technology Corporation | Extended 1D Inversion Of Electromagnetic Measurements For Formation Evaluation |
| CN105893674A (en) * | 2016-03-31 | 2016-08-24 | 恒泰艾普石油天然气技术服务股份有限公司 | Method for performing geological attribute prediction with global covariance |
| CN106353811A (en) * | 2016-08-12 | 2017-01-25 | 北京天工新源科技有限公司 | Method for identifying fluid of multi-attribute fusion reservoirs on basis of well information attribute optimization |
| CN107703544A (en) * | 2017-09-27 | 2018-02-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Oil gas forecasting method is changed with offset distance based on the indication using prestack seismic amplitude of geostatistics |
| CN108073765A (en) * | 2017-12-11 | 2018-05-25 | 中海石油(中国)有限公司湛江分公司 | A kind of horizontal well routine well logging bed boundary identification is with border away from inversion method |
-
2019
- 2019-01-03 CN CN201910004452.XA patent/CN109633781B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100228484A1 (en) * | 2007-12-06 | 2010-09-09 | Exxon-Mobile, Upstream Research Company | Volume of Investigation Based Image Processing |
| US20160047934A1 (en) * | 2013-04-02 | 2016-02-18 | Schlumberger Technology Corporation | Extended 1D Inversion Of Electromagnetic Measurements For Formation Evaluation |
| CN105893674A (en) * | 2016-03-31 | 2016-08-24 | 恒泰艾普石油天然气技术服务股份有限公司 | Method for performing geological attribute prediction with global covariance |
| CN106353811A (en) * | 2016-08-12 | 2017-01-25 | 北京天工新源科技有限公司 | Method for identifying fluid of multi-attribute fusion reservoirs on basis of well information attribute optimization |
| CN107703544A (en) * | 2017-09-27 | 2018-02-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Oil gas forecasting method is changed with offset distance based on the indication using prestack seismic amplitude of geostatistics |
| CN108073765A (en) * | 2017-12-11 | 2018-05-25 | 中海石油(中国)有限公司湛江分公司 | A kind of horizontal well routine well logging bed boundary identification is with border away from inversion method |
Non-Patent Citations (2)
| Title |
|---|
| 王建 等: "基于测井数据的三维地质模型构建与可视化", 《测井技术》 * |
| 邵才瑞 等: "基于测井数据的地质曲面插值重构方法比较", 《测井技术》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112016048A (en) * | 2020-07-28 | 2020-12-01 | 中交广州航道局有限公司 | Geological attribute prediction method based on multi-point geological statistics |
| CN113361263A (en) * | 2021-06-04 | 2021-09-07 | 中国人民解放军战略支援部队信息工程大学 | Character entity attribute alignment method and system based on attribute value distribution |
| CN113361263B (en) * | 2021-06-04 | 2023-10-20 | 中国人民解放军战略支援部队信息工程大学 | Character entity attribute alignment method and system based on attribute value distribution |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109633781B (en) | 2021-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109800863B (en) | Logging phase identification method based on fuzzy theory and neural network | |
| Woo | Kernel estimation methods for seismic hazard area source modeling | |
| AU2013325186B2 (en) | Propagating fracture plane updates | |
| CN111079217B (en) | BIM-based geotechnical engineering comprehensive investigation information interpretation method and system | |
| CN103886216B (en) | A kind of multiple spot Geo-statistic Method based on geology Vector Message | |
| CN104854482A (en) | Rock facies prediction in non-cored wells from cored wells | |
| CN108645994A (en) | A kind of geology stochastic inversion methods and device based on Multiple-Point Geostatistics | |
| CN105044775B (en) | A kind of earthquake fluid inverting recognition methods and device | |
| CN108291439A (en) | The method that earth model is selected from multiple earth models | |
| CN105205239A (en) | Method and device for modeling reservoir physical property parameter | |
| CN109633781A (en) | Geological property acquisition methods and device, electronic equipment, storage medium | |
| RU2720115C1 (en) | Method of automated geological survey of wells and system for its implementation | |
| CN107316341A (en) | A kind of Multiple-Point Geostatistics facies modelization method | |
| CN114332391A (en) | Three-dimensional geologic body probability model modeling method | |
| CN103543478A (en) | Geologic morphological interpolation KM (Kriging and Multiple-point geostatistics) method | |
| CN102830430B (en) | A kind of horizon velocity modeling method | |
| AU2012385936B2 (en) | Systems and methods for estimating opportunity in a reservoir system | |
| US11899150B2 (en) | Velocity model for sediment-basement interface using seismic and potential fields data | |
| Yin et al. | A technical note on the application of interpolation method on modeling of 3D surface rock layer | |
| CN114153002A (en) | Three-dimensional geological modeling method and device for natural fracture of reservoir, electronic equipment and medium | |
| US20190265386A1 (en) | Meshless and Mesh-Based Technique for Modeling Subterranean Volumes | |
| Jia et al. | Accurate method for picking up the first arrival time of microseismic signals based on entropy theory | |
| CN108873059A (en) | Construct the method and device of loess thickness model | |
| CN118037980B (en) | Stratum determination method and device based on geospatial interpolation, medium and equipment | |
| CN117744454B (en) | Method and system for predicting equivalent scouring depth of marine structure |
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 |