CN109633781A - Geological property acquisition methods and device, electronic equipment, storage medium - Google Patents

Geological property acquisition methods and device, electronic equipment, storage medium Download PDF

Info

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
Application number
CN201910004452.XA
Other languages
Chinese (zh)
Other versions
CN109633781B (en
Inventor
黄瑞
昝成
龙威
程浩然
李炜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Cleanergy Aike (shenzhen) Energy Technology Co Ltd
Shenzhen Research Institute Tsinghua University
Original Assignee
Cleanergy Aike (shenzhen) Energy Technology Co Ltd
Shenzhen Research Institute Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Cleanergy Aike (shenzhen) Energy Technology Co Ltd, Shenzhen Research Institute Tsinghua University filed Critical Cleanergy Aike (shenzhen) Energy Technology Co Ltd
Publication of CN109633781A publication Critical patent/CN109633781A/en
Application granted granted Critical
Publication of CN109633781B publication Critical patent/CN109633781B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V9/00Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (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

Geological property acquisition methods and device, electronic equipment, storage medium
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.
CN201910004452.XA 2018-08-16 2019-01-03 Geological property acquisition method and device, electronic equipment and storage medium Active CN109633781B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
王建 等: "基于测井数据的三维地质模型构建与可视化", 《测井技术》 *
邵才瑞 等: "基于测井数据的地质曲面插值重构方法比较", 《测井技术》 *

Cited By (3)

* Cited by examiner, † Cited by third party
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
CN107563574A (en) A kind of earth's magnetic field Forecasting Methodology based on long memory models Recognition with Recurrent Neural Network in short-term
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
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
CN118037980B (en) Stratum determination method and device based on geospatial interpolation, medium and equipment
CN110927793B (en) Reservoir prediction method and system based on sequential random fuzzy simulation
CN107316341A (en) A kind of Multiple-Point Geostatistics facies modelization method
CN106446514A (en) Fuzzy theory and neural network-based well-log facies recognition method
CN114332391A (en) Three-dimensional geologic body probability model modeling method
CN103543478A (en) Geologic morphological interpolation KM (Kriging and Multiple-point geostatistics) method
CN116757323A (en) Landslide displacement prediction method based on IVMD-IAO-BiLSTM
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
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
CN117759235B (en) Method for realizing mineral exploration based on limited drilling
CN114706126B (en) Quantitative carving method and system for seismic geologic body under geological awareness constraint

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