CN115113281A - Method, device, equipment and storage medium for determining marker well and marker layer thereof - Google Patents
Method, device, equipment and storage medium for determining marker well and marker layer thereof Download PDFInfo
- Publication number
- CN115113281A CN115113281A CN202210540165.2A CN202210540165A CN115113281A CN 115113281 A CN115113281 A CN 115113281A CN 202210540165 A CN202210540165 A CN 202210540165A CN 115113281 A CN115113281 A CN 115113281A
- Authority
- CN
- China
- Prior art keywords
- well
- similarity
- marker
- stratum
- wells
- 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
- 239000003550 marker Substances 0.000 title claims abstract description 108
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 34
- 239000011159 matrix material Substances 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 230000000877 morphologic effect Effects 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 5
- 238000012217 deletion Methods 0.000 claims description 4
- 230000037430 deletion Effects 0.000 claims description 4
- 238000005755 formation reaction Methods 0.000 abstract description 29
- 230000000694 effects Effects 0.000 abstract description 2
- 230000000875 corresponding effect Effects 0.000 description 84
- 238000010586 diagram Methods 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
Abstract
The method comprises the steps of obtaining well tracks, well logging curves and formation data of all wells, and obtaining first similarity of the same well logging curves of the formations corresponding to every two wells according to the well tracks, the formation data and the well logging curves; selecting a marker well from each well by adopting a greedy search idea according to the first similarity; and selecting the marker layer of the marker well according to the first similarity. This application has the effect that accurate efficient confirmed sign well and mark layer thereof.
Description
Technical Field
The present application relates to the technical field of formation data processing, and in particular, to a marker well, a method, an apparatus, a device, and a storage medium for determining a marker layer of the marker well.
Background
The determination of the marker well and the marker layer thereof is used as a core technology for fine oil and gas reservoir description, is the basis for developing geological research of the oil and gas reservoir, and is also the key premise for subsequently making scientific and efficient oil and gas reservoir development schemes. In the face of a large number of scattered oil and gas reservoirs, and thousands of wells in each block of the oil and gas reservoir, the conventional method is to manually observe the data of each well in the block to realize the selection of the marker wells and the marker layers thereof.
The traditional manual method seriously depends on experience recognition of service personnel on a specific block, consumes much time and energy of the service personnel, has the problems of strong subjectivity, long time consumption, inconsistent results and the like, and brings huge working pressure to oil reservoir evaluation personnel, oil reservoir research personnel and well logging interpreters. Therefore, an intelligent and automatic method is urgently needed to replace the traditional manual method, so that the workload of personnel is reduced, and the efficiency of determining the marking well and the marking layer thereof is improved.
Disclosure of Invention
In order to accurately and efficiently determine a marker well and a marker layer thereof, the application provides a method, a device, equipment and a storage medium for determining the marker well and the marker layer thereof.
In a first aspect, the present application provides a marker well and a method for determining a marker layer thereof, which adopts the following technical scheme:
a method of determining a marker well and its marker formation, comprising:
acquiring a well track, a logging curve and formation data of each well, and acquiring first similarity of the same logging curve of the formation corresponding to each two wells according to the well track, the formation data and the logging curve;
selecting a marker well from each well by adopting a greedy search idea according to the first similarity;
and selecting the marker layer of the marker well according to the first similarity.
By adopting the technical scheme and adopting a greedy algorithm, the first similarity is processed, and the marker well can be intelligently and automatically selected; in addition, the marking layer of the marking well can be automatically selected through the first similarity, and the problems of high working strength, low efficiency, poor repeatability and the like in the manual determination method in the prior art are solved.
Preferably, the first similarity includes a mean similarity, a variance similarity and a morphology similarity; the method for acquiring the well track, the logging curve and the stratum data of each well, and acquiring the first similarity of the same logging curve of the stratum corresponding to every two wells according to the well track, the stratum data and the logging curve comprises the following steps:
preprocessing the logging curve, acquiring depths corresponding to all stratums of a target well according to the stratum data, acquiring observation points corresponding to all stratums in a well track according to the depths corresponding to all the stratums of the target well, and acquiring a curve value corresponding to each observation point on the preprocessed logging curve;
obtaining the mean similarity and variance similarity of the same logging curve of the stratum corresponding to each two wells according to the curve values;
and adopting a DTW algorithm to obtain the form similarity of the same logging curves of the stratum corresponding to every two wells.
By adopting the technical scheme, the logging curve is preprocessed, and the preprocessed logging curve data is more complete and accurate; matching the formation data, the well track and the logging curve according to the depth, so as to facilitate the subsequent calculation of the first similarity; the first similarity is obtained according to three different similarities, namely mean similarity, variance similarity and form similarity, so that the obtained first similarity is more representative, the similarity between wells can be more accurately represented, and the foundation technology construction is made for subsequently selecting the marker well and the marker layer thereof.
Preferably, the preprocessing the well log comprises: longitudinal splicing, invalid value deletion, curve attribution, structuralization, logarithm taking processing and normalization processing.
Preferably, the obtaining of the morphological similarity of the same kind of logging curves of the stratum corresponding to each two wells by using the DTW algorithm includes:
acquiring an observation point distance matrix of the stratum corresponding to each two wells by adopting a DTW algorithm;
carrying out boundary constraint on the observation point distance matrix, backtracking a non-boundary constraint part of the observation point distance matrix, and obtaining a shortest path;
and according to the shortest path, acquiring the form similarity of the same logging curves of the stratums corresponding to the two wells.
By adopting the technical scheme, the shortest path is always in the diagonal part of the matrix according to the experiment, so that the boundary constraint is added, the backtracking is only carried out on the non-boundary constraint part, and the boundary constraint part does not participate in the backtracking calculation, thereby reducing the search space and reducing the calculation amount.
Preferably, the selecting a marker well from the wells according to the first similarity by using a greedy search concept includes:
acquiring a second similarity of the logging curve of each well according to the first similarity;
setting a first threshold value according to a sigma principle, and selecting a selection curve according to the second similarity and the first threshold value;
acquiring a third similarity of the two wells according to the similarity of the selection curves corresponding to the two known wells;
and setting a second threshold according to a sigma principle, and selecting the marker well by using a greedy thought according to the third similarity and the second threshold.
By adopting the technical scheme, a marking well selection method of greedy search is adopted, and wells corresponding to well logging curves with good curve forms are screened as marking wells.
Preferably, the selecting the marker layer of the marker well according to the first similarity includes:
and acquiring the stratum similarity of the stratum of the marker well according to the first similarity, and selecting the marker layer of the marker well according to the stratum similarity.
In a second aspect, the present application provides a marker well and a device for determining a marker layer thereof, which adopts the following technical solutions:
a marking well and a device for determining a marking layer thereof comprise,
the acquisition module is used for acquiring well tracks, logging curves and stratum data of all wells, and acquiring first similarity of the same logging curves of the stratum corresponding to every two wells according to the well tracks, the stratum data and the logging curves;
the selecting module is used for selecting a marker well from each well by adopting a greedy search idea according to the first similarity; and the number of the first and second groups,
and the selecting module is used for selecting the marking layer of the marking well according to the first similarity.
In a third aspect, the present application provides a marker well and a device for determining a marker bed thereof, which adopts the following technical scheme:
a marker well and its marker layer determination apparatus comprising a memory and a processor, the memory having stored thereon a computer program that can be loaded by the processor and that performs the method of marker well and its marker layer determination of any one of the first aspect.
In a fourth aspect, the present application provides a computer-readable storage medium, which adopts the following technical solutions:
a computer readable storage medium storing a computer program capable of being loaded by a processor and executing the method of determining a marker well and its marker bed according to any of the first aspect.
Drawings
Fig. 1 is a schematic flow chart of the determination of a marker well and its marker layer provided in the embodiments of the present application.
Fig. 2 is a schematic diagram of continuous partial data provided by an embodiment of the present application.
Fig. 3 is a schematic diagram of a first similarity provided in an embodiment of the present application.
Fig. 4 is a schematic diagram of a third mean value, a fourth mean value and a fifth mean value provided in an embodiment of the present application.
Fig. 5 is a schematic diagram of a third similarity provided in the embodiment of the present application.
Fig. 6 is a schematic diagram of a similarity matrix provided in an embodiment of the present application.
Fig. 7 is a block diagram of a marker well and a marker layer determination device thereof according to an embodiment of the present application.
Fig. 8 is a schematic structural diagram of a marker well and a marker layer determination device thereof provided in an embodiment of the present application.
Fig. 9 is a schematic diagram of an observation point distance matrix provided in an embodiment of the present application.
Fig. 10 and fig. 11 are schematic diagrams of a process of finding a shortest path according to an embodiment of the present application.
FIG. 12 is a schematic diagram of a boundary constraint section provided in an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application.
The present embodiment provides a method for determining a marker well and a marker layer thereof, as shown in fig. 1, the main flow of the method is described as follows (steps S101 to S103):
step S101: and acquiring the well track, the logging curve and the stratum data of each well, and acquiring the first similarity of the logging curves of the same kind of stratum corresponding to every two wells according to the well track, the stratum data and the logging curve.
Where each well may be all wells in a block of the reservoir.
The method for acquiring the well track of the well comprises the following specific steps:
and acquiring the wellhead data of each well and the XY offset in the well deviation data, and calculating to obtain the XY coordinates corresponding to the current vertical depth according to the XY offset of the vertical depth in the well deviation data relative to the wellhead data.
The current vertical depth is Z, and the corresponding coordinates are (X, Y). Wherein X ═ X Well head +X Offset of ；Y＝Y Well head +Y Offset of (ii) a Wherein the wellhead data comprises X Well head And Y Well head (ii) a The XY offset comprises X Offset of And Y Offset of 。
The well track comprises a plurality of observation points and coordinates of the observation points, and the observation points are connected in sequence according to the coordinates of the observation points to obtain the well track of the well. Wherein the coordinates of the observation point are (X, Y, Z).
The method comprises the steps of obtaining a logging curve and preprocessing the logging curve, wherein the preprocessing comprises longitudinal splicing, invalid value deletion, curve attribution, structuring, log taking processing (namely log processing) and normalization processing. The preprocessing sequence can be longitudinal splicing, invalid value deletion, curve attribution, structuralization, logarithm processing and normalization processing in sequence.
Wherein, vertically splice: acquiring logging curves of different series of wells, and if the logging curves have data files corresponding to different depth sections, longitudinally splicing the logging curves according to the depth, namely splicing the data files of the different depth sections according to the depth value sequence; if the well log is a data file of a complete depth segment, no longitudinal concatenation is required.
Deleting invalid values: and deleting invalid values in each log according to the continuity and the size of the numerical values. Such as-999.25, -9999, and null, etc.
Attribution of curves: and attributing the logging curves to corresponding wells, and associating the logging curves belonging to the same well. For example, 8 logging curves (in the present application, each of the multiple logging curves of one well is a unique type; the type may refer to a name; hereinafter, the type of the logging curve may refer to a name) are obtained, and after analysis, 3 logging curves are curves of a well a, the 3 logging curves are classified into the well a, and the other 5 logging curves are curves of a well b, the 5 logging curves are classified into the well b; correlating 3 logging curves of the well a one by one according to depths, namely, corresponding all values on the 3 logging curves one by one according to the depths; similarly, the 5 well logs of well b are correlated one by one according to depth.
Structuring: and acquiring the stratum data of the well, and correspondingly matching the stratum data, the well track and the logging curve of the well. Specifically, the depth measurement in the well log corresponds to the vertical depth in the well trace, each observation point of the well trace can find a corresponding well log value on the corresponding well log, and the depth value in the formation data (the depth value in the formation data is also the depth measurement) corresponds to the depth measurement in the well log and the vertical depth in the well trace respectively.
And (3) logarithm processing: the resistance is a part of the data in the logging curve data, and the resistance is subjected to logarithmic processing according to the artificial stratum contrast process.
Normalization treatment: in order to ensure that the dimensions of the data of the logging curves are consistent when the similarity is calculated subsequently and the comparison is convenient, and in order to prevent the influence of non-uniformity of the dimensions on the subsequent target well marking stratum dividing effect, the logging curves are subjected to normalization processing.
Referring to fig. 2, in all data of the log, continuous data of the formation part of the well corresponding to the log are intercepted: acquiring the top depth of a first stratum and the bottom depth of a last stratum, and solving the average thickness of each stratum; acquiring an upper limit and a lower limit according to the top depth of the first stratum, the bottom depth of the last stratum and the average thickness, wherein the upper limit is the difference of the top depth of the first stratum minus the average thickness, and the lower limit is the difference of the bottom depth of the last stratum plus the average thickness; first data corresponding to the upper limit on the logging curve and second data corresponding to the lower limit on the logging curve are obtained, and the continuous partial data are all data (including the first data and the second data) between the first data and the second data. Respectively solving a maximum value and a minimum value in continuous partial data, and normalizing the logging curve data according to the maximum value and the minimum value, wherein a specific formula is as follows:
wherein x is part For successive partial data, x i The corresponding logging curve value of the depth value on the logging curve is obtained; x is the number of j Is normalized x i 。
By the formula, the normalization processing of the whole logging curve is completed, and the influence of the head and tail non-stratum parts on final data is reduced.
The first similarity comprises mean similarity, variance similarity and form similarity.
Mean similarity: and acquiring the depth corresponding to each stratum of the well according to the stratum data, and acquiring a curve value corresponding to each depth value on the preprocessed logging curve according to the depth corresponding to the stratum (namely the curve values corresponding to all observation points corresponding to the stratum on the preprocessed logging curve).
And respectively acquiring the mean value similarity of the same logging curves of the stratums corresponding to every two wells (the corresponding stratums refer to two stratums with the same name, and the same logging curves refer to logging curves with the same name) according to the corresponding curve values of the observation points on the logging curves. Specifically, the formation corresponding to two wells refers to a formation in which the two wells have the same formation name. Referring to fig. 3, the first column and the second column respectively refer to the names of two wells, the third column refers to the names of the stratums corresponding to the two wells, and the fourth column refers to the names of the log curves of the same kind of the two wells. Respectively defining the stratums corresponding to the two wells as a first corresponding stratum and a second corresponding stratum, calculating the mean value of curve values of all observation points of the first corresponding stratum on the corresponding logging curve, and defining the mean value as a first mean value; similarly, calculating the mean value of curve values of all the observation points of the second corresponding stratum on the corresponding logging curve, and defining the mean value as a second mean value; the mean value similarity of the same kind of logging curves of the stratums corresponding to the two wells is the difference between the first mean value and the second mean value. The smaller the value of the mean similarity is, the higher the similarity of the formations corresponding to the two wells is, and the logging curve segments of the formations corresponding to the two wells are considered to be within the same fluctuation threshold. It should be noted that the two well logs used for calculating the first average and the second average are the same kind of well log.
Variance similarity: the method for calculating the variance similarity is consistent with the method for calculating the mean similarity in principle, the variance of the curve values of all the observation points of the first corresponding stratum on the corresponding logging curve is calculated, and the variance is defined as a first variance; similarly, calculating the variance of the curve values of all the observation points of the second corresponding stratum on the corresponding logging curve, and defining the variance as a second variance; the variance similarity of the same kind of well logs of the corresponding stratum of the two wells is the difference of the first variance and the second variance. The smaller the value of the variance similarity is, the higher the similarity of the formations corresponding to the two wells is, and the fluctuation degrees of the logging curve sections of the formations corresponding to the two wells tend to be the same. It should be noted that the two logs used to calculate the first variance and the second variance are the same kind of log.
Morphological similarity: and adopting a DTW algorithm to obtain the form similarity of the same logging curves of the stratum corresponding to every two wells. The method comprises the following specific steps:
(1) and calculating the distance between the observation points of the stratums corresponding to the two wells to obtain an observation point distance matrix. The formula for solving the distance between the observation points of the stratums corresponding to the two wells is as follows:
d(i，j)＝|one i -two j |；
wherein d (i, j) is the distance between the observation points of the stratum corresponding to the two wells; one (E) i Representing a curve value corresponding to an observation point i of a well stratum (namely a first corresponding stratum) on a logging curve of the well; two (two) j And representing the curve value corresponding to the observation point j of the formation corresponding to the other well (namely, the second corresponding formation) on the logging curve of the well.
The observation point distance matrix is shown in fig. 9.
(2) The distance of each possible path is calculated in the observation point distance matrix, noting that the selected path must start from the lower left square of the observation point distance matrix and end at the upper right square, and that each square cannot be traversed repeatedly. When the distance from any square in ((i-1, j-1), (i-1, j) or (i, j-1)) to the next square (i, j) is d (i, j) if the distance is horizontal or vertical, and is 2d (i, j) if the distance is diagonal.
The calculation formula of the distance g (i, j) between the observation points of the stratums corresponding to the two wells is as follows:
(3) and (5) backtracking the observation point distance matrix by adopting a backtracking method to obtain the shortest path.
The process of finding the shortest path is described with reference to fig. 10 and 11.
(4) And according to the shortest path, acquiring the form similarity of the same logging curves of the stratums corresponding to the two wells (the logging curve of one stratum is a section of curve on the whole logging curve of the well).
Redefining the calculation mode of the shortest path, so that under the condition that the lengths of the same type of well logs of the stratum of one well and the stratum corresponding to the other well are closer (the closer the lengths are, the closer the number of observation points corresponding to the well log of one stratum and the number of observation points corresponding to the well log of the other stratum are), the smaller the acquired form similarity is, and the well logs of the stratum of one well and the stratum corresponding to the other well are more similar.
The formula for solving the morphological similarity best _ dis of the same type of logging curves of the stratums corresponding to the two wells is as follows:
best_dis＝g(m，n)/min(m，n)；
wherein m and n are the lengths of the logging curves of the same type of the well stratum and the stratum corresponding to the other well respectively, and g (m and n) is the distance value of the shortest path obtained in the step (3).
Since the measurement mode of the morphological similarity is measured according to the distance between the observation points of the two wells, the smaller the value of the morphological similarity is, the higher the similarity of the stratum corresponding to the two wells is.
Optionally, an improvement is performed on the basis of the DTW algorithm, specifically, in the step (4), the boundary constraint is performed on the observation point distance matrix, and a backtracking method is adopted to backtrack the non-boundary constraint part of the observation point distance matrix to obtain the shortest path. Selecting a right triangle at the upper left corner of the observation point distance matrix as a shadow triangle, and taking N percent of the abscissa from left to right 1 Taking N percent of the ordinate from top to bottom as a rectangular length of the shadow triangle 2 Another length of right angle as a shadow triangle(ii) a Similarly, a right triangle is also selected from the lower right corner of the observation point distance matrix as a shadow triangle, and N percent of the abscissa is taken from right to left 3 Taking N percent of ordinate from bottom to top as a length of a right angle of the shadow triangle 4 As another right angle length of the shaded triangle. In this example, N 1 、N 2 、 N 3 And N 4 Are all twenty.
The area of the two hatched triangles is the boundary constraint part, see fig. 12.
Step S102: and selecting the marker wells from the wells according to the first similarity by adopting a greedy search idea.
(1) And acquiring a second similarity of the log curve of each well according to the first similarity:
calculating the mean value of the mean value similarity of the same well logging curves of all the stratums, and defining the mean value as a third mean value; calculating the mean value of the variance similarity of the same well logging curves of all the stratums of the well, and defining the mean value as a fourth mean value; and calculating the mean value of the morphological similarity of the same well log of all the stratums of the well, and defining the mean value as a fifth mean value. Referring to fig. 4, the first column of fig. 4 is a fifth mean value, which is the mean value of the morphological similarity of CBL _ cementing quality logs (log names) of all the formations of the two wells in fig. 1, and the second column of fig. 4 is a third mean value, which is the mean value of the mean similarity of CBL _ cementing quality logs (log names) of all the formations of the two wells in fig. 1; the third column of FIG. 4 is the fourth mean value, which is the mean value of the variance similarity of CBL _ Wellmass logs (log names) of all the formations of the two wells in FIG. 1.
And respectively carrying out weighted summation on the third average value, the fourth average value and the fifth average value of the logging curves of the same kind of the two wells to obtain the second similarity of the logging curves. Referring to FIG. 4, the result of weighted summation of the third, fourth and fifth means between two wells is a second similarity.
(2) Setting a first threshold value according to a sigma principle, and selecting a selection curve according to the second similarity and the first threshold value:
the first threshold is a difference between a mean of the second similarities of all the logs and a variance of the second similarities of all the logs. Judging whether the second similarity is smaller than a first threshold value; if so, taking the logging curve corresponding to the second similarity as the selection curve.
(3) According to the similarity of the selection curves corresponding to the two known wells, acquiring the third similarity of the two wells:
wherein, the similarity of the selected curve is the corresponding second similarity. Referring to fig. 5, in all the same well logging curves of two wells, the average value of the similarity of the corresponding selected curves is selected, and the average value is the final similarity of the two wells, namely the third similarity.
(4) Setting a second threshold value according to a sigma principle, selecting a marker well according to a greedy search idea and according to a third similarity and the second threshold value:
referring to FIG. 6, a similarity matrix between wells is created based on the third similarity.
Wherein the second threshold is a difference between a mean of all the third similarities and a variance of all the third similarities (all the third similarities refer to all the values in the similarity matrix).
The method comprises the following steps: and selecting the minimum third similarity in the similarity matrix.
Step two: judging whether the selected minimum third similarity is smaller than a second threshold value or not; if yes, executing the third step; if not, executing the step six.
Step three: and selecting the two wells corresponding to the minimum third similarity as the marking wells.
Step four: judging whether the number of the marker wells is larger than a first preset value and smaller than a second preset value; if so, stopping selecting the marker wells; if not, executing the fifth step; the first preset value is P percent of the number of the reference wells, the second preset value is Q percent of the number of the reference wells, Q is greater than P, P is ten and Q is fifty in the embodiment.
Step five: and (4) taking the row and the column where the minimum third similarity exists (the row and the column where the minimum third similarity belongs in the similarity matrix) as the current row and column, removing the minimum third similarity in the current row and column, continuously selecting the minimum third similarity in the remaining third similarities of the current row and column, and returning to the step two.
Step six: and (4) removing the row and the column where the historical minimum third similarity is located (the historical minimum third similarity is all the minimum third similarities related to the first step to the fifth step), selecting the minimum third similarity from the rest rows and columns, and returning to the second step.
It is noted that the third similarity at the same position in the similarity matrix is not selected repeatedly.
After the marked wells are selected, the remaining wells except the marked wells are non-marked wells in all the wells.
Step S103: and selecting the marker layer of the marker well according to the first similarity.
The method comprises the steps that a marked well is provided with a plurality of stratums, and the mean similarity, the variance similarity and the form similarity of all logging curve sections of each stratum (all logging curve sections of a stratum are curve sections corresponding to all logging curves of one stratum respectively) are obtained respectively; carrying out weighted summation on the mean similarity, the variance similarity and the form similarity belonging to the same logging curve segment to obtain a sum value; and calculating the average value of the sum values of all the logging curve segments of the stratum, wherein the average value is the stratum similarity of the stratum.
And selecting a marking layer of the marking well according to the stratum similarity. For example, in all the stratums of the marker well, the stratum similarity is sorted from small to large, and the stratums corresponding to the first n stratum similarities are selected as the marker layers of the marker well in the sorting result; for another example, a first preset threshold is set, whether the stratum similarity is smaller than the first preset threshold is judged, and if yes, the stratum corresponding to the stratum similarity is used as a marker layer of the marker well.
Optionally, after dividing the marker well into a plurality of strata according to the stratum data, obtaining the thickness of the strata before obtaining the observation points of the strata corresponding to the well trajectory, and judging whether the thicknesses of a plurality of continuous strata are all smaller than a thickness threshold; if so, it indicates that the thickness of the multiple continuous strata is relatively thin, and the characteristics of the relatively thin strata are not sufficiently obvious, so the multiple continuous strata with relatively thin thickness can also be called small strata, and the multiple continuous small strata are combined into one stratum, so that the characteristics of the strata are sufficiently obvious.
Optionally, if the marker well is deep, combining a plurality of continuous small layers into one stratum; if the marker well is shallow, there is no need to merge multiple successive sub-layers into one formation, since the shallower marker well itself has no several formations to combine.
Correspondingly, the method for selecting the marker layer of the marker well can also comprise the following steps: and calculating the similarity of the same logging curves of the two marker wells corresponding to the small layers, wherein the similarity also comprises mean similarity, variance similarity and form similarity, and calculating the stratum similarity of the small layers (the calculation method of the stratum similarity of the small layers is consistent with the calculation method of the stratum similarity of the medium stratum in principle, and is not repeated here). And selecting a plurality of small layers with smaller stratum similarity from the stratum similarities of all the small layers of the mark layer to be recombined into a new mark layer. For example, in all the small layers of the mark layer, the stratum similarities of the small layers are sorted from small to large, the small layers corresponding to the first m stratum similarities are selected in the sorting result, and the selected small layers are recombined into a new mark layer; for another example, a second preset threshold is set, and the small layers corresponding to the stratum similarity of the small layers smaller than the second preset threshold are recombined into a new mark layer.
In order to better implement the method, the embodiment of the present application further provides a device for determining a marker well and a marker layer thereof, which may be specifically integrated in a device for determining a marker well and a marker layer thereof, such as a terminal or a server, where the terminal may include, but is not limited to, a mobile phone, a tablet computer, or a desktop computer.
Fig. 7 is a block diagram of a structure of a marker well and a marker layer determination apparatus thereof according to an embodiment of the present application, and as shown in fig. 7, the apparatus mainly includes:
the acquisition module 201 is configured to acquire a well trajectory, a logging curve and formation data of each well, and acquire a first similarity of the same logging curve of the formation corresponding to each two wells according to the well trajectory, the formation data and the logging curve;
a selecting module 202, configured to select a marker well from each well according to the first similarity by using a greedy search idea; and the number of the first and second groups,
and the selecting module 203 is used for selecting the marking layer of the marking well according to the first similarity.
Various changes and specific examples in the method provided in the above embodiments are also applicable to the marking well and the device for determining the marking layer thereof in this embodiment, and through the foregoing detailed description of the marking well and the method for determining the marking layer thereof, those skilled in the art can clearly know the method for implementing the marking well and the device for determining the marking layer thereof in this embodiment, and for the sake of brevity of the description, details are not described here.
In order to better execute the program of the method, the embodiment of the present application further provides a device for determining a marker well and a marker layer thereof, as shown in fig. 8, the device 300 for determining a marker well and a marker layer thereof includes a memory 301 and a processor 302.
The marker well and its marker layer determination device 300 may be implemented in various forms including devices such as cell phones, tablets, palm top computers, laptop and desktop computers.
The memory 301 may be used to store, among other things, instructions, programs, code sets, or instruction sets. The memory 301 may include a program storage area and a data storage area, wherein the program storage area may store instructions for implementing an operating system, instructions for at least one function (for example, obtaining a first similarity, a second similarity, a third similarity, and the like), and instructions for implementing the method for determining the marker wells and the marker layers thereof provided in the foregoing embodiments, and the like; the storage data area may store data and the like involved in the method for determining the marker wells and the marker layers thereof provided in the above embodiments.
An embodiment of the present application provides a computer-readable storage medium, including: a U-disk, a removable hard disk, a Read Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, an optical disk, or other various media capable of storing program codes. The computer readable storage medium stores a computer program that can be loaded by a processor and executed to perform the method of determining a marker well and its marker bed of the above embodiments.
The specific embodiments are merely illustrative and not restrictive, and various modifications that do not materially contribute to the embodiments may be made by those skilled in the art after reading this specification as required, but are protected by patent laws within the scope of the claims of this application.
Claims (9)
1. A method of determining a marker well and its marker formation, comprising:
acquiring a well track, a logging curve and formation data of each well, and acquiring first similarity of the same logging curve of the formation corresponding to each two wells according to the well track, the formation data and the logging curve;
selecting a marker well from each well by adopting a greedy search idea according to the first similarity;
and selecting the marker layer of the marker well according to the first similarity.
2. The method of claim 1, wherein the first similarity comprises a mean similarity, a variance similarity, and a morphology similarity; the method for acquiring the well track, the logging curve and the stratum data of each well and acquiring the first similarity of the logging curves of the same kind of stratum corresponding to every two wells according to the well track, the stratum data and the logging curve comprises the following steps:
preprocessing the logging curve, acquiring depths corresponding to all stratums of a target well according to the stratum data, acquiring observation points corresponding to all stratums in a well track according to the depths corresponding to all the stratums of the target well, and acquiring a curve value corresponding to each observation point on the preprocessed logging curve;
obtaining the mean similarity and variance similarity of the same logging curve of the stratum corresponding to each two wells according to the curve values;
and adopting a DTW algorithm to obtain the form similarity of the same logging curves of the stratum corresponding to every two wells.
3. The method of claim 2, wherein the pre-processing the well log comprises: longitudinal splicing, invalid value deletion, curve attribution, structurization, logarithm processing and normalization processing.
4. The method of claim 2, wherein the obtaining the morphological similarity of the same log curves of the formation corresponding to each two wells by using the DTW algorithm comprises:
acquiring an observation point distance matrix of the stratum corresponding to each two wells by adopting a DTW algorithm;
carrying out boundary constraint on the observation point distance matrix, backtracking a non-boundary constraint part of the observation point distance matrix, and obtaining a shortest path;
and according to the shortest path, acquiring the form similarity of the same logging curves of the stratums corresponding to the two wells.
5. The method of claim 1, wherein selecting a marker well from the wells based on the first similarity using a greedy search concept comprises:
acquiring a second similarity of the logging curve of each well according to the first similarity;
setting a first threshold according to a sigma principle, and selecting a selection curve according to the second similarity and the first threshold;
acquiring a third similarity of the two wells according to the similarity of the selection curves corresponding to the two known wells;
and setting a second threshold value according to a sigma principle, and selecting the marker well according to the third similarity and the second threshold value.
6. The method of claim 1, wherein selecting a marker layer of the marker well based on the first similarity comprises:
and acquiring the stratum similarity of the stratum of the marker well according to the first similarity, and selecting the marker layer of the marker well according to the stratum similarity.
7. A device for determining a marker well and a marker layer thereof is characterized by comprising,
the acquisition module is used for acquiring well tracks, logging curves and stratum data of all wells, and acquiring first similarity of the same logging curves of the stratum corresponding to every two wells according to the well tracks, the stratum data and the logging curves;
the selecting module is used for selecting a marker well from each well by adopting a greedy search idea according to the first similarity; and the number of the first and second groups,
and the selecting module is used for selecting the marking layer of the marking well according to the first similarity.
8. A marker well and marker layer determination apparatus therefor, comprising a memory and a processor, the memory having stored thereon a computer program which can be loaded by the processor and which performs the method of any one of claims 1 to 6.
9. A computer-readable storage medium, in which a computer program is stored which can be loaded by a processor and which executes the method of any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210540165.2A CN115113281B (en) | 2022-05-18 | 2022-05-18 | Method, device, equipment and storage medium for determining marker well and marker layer thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210540165.2A CN115113281B (en) | 2022-05-18 | 2022-05-18 | Method, device, equipment and storage medium for determining marker well and marker layer thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115113281A true CN115113281A (en) | 2022-09-27 |
| CN115113281B CN115113281B (en) | 2023-01-03 |
Family
ID=83327210
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210540165.2A Active CN115113281B (en) | 2022-05-18 | 2022-05-18 | Method, device, equipment and storage medium for determining marker well and marker layer thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115113281B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116201535A (en) * | 2023-02-06 | 2023-06-02 | 北京月新时代科技股份有限公司 | Automatic dividing method, device and equipment for oil and gas reservoir target well sign stratum |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110852144A (en) * | 2019-08-23 | 2020-02-28 | 中国石油大学(北京) | DTW-based intelligent stratum comparison method and system |
| CN111472765A (en) * | 2020-04-24 | 2020-07-31 | 北京月新时代科技股份有限公司 | Stratigraphic division method and device for target well |
| CN112459767A (en) * | 2020-10-30 | 2021-03-09 | 中国石油天然气集团有限公司 | Standard layer selection method and device |
| CN113236124A (en) * | 2021-05-07 | 2021-08-10 | 四川页岩气勘探开发有限责任公司 | Deep shale gas horizontal well geological guiding method |
| CN113420936A (en) * | 2021-07-01 | 2021-09-21 | 中国海洋石油集团有限公司 | Stratum comparison method based on dynamic time programming optimization algorithm |
| US20220035065A1 (en) * | 2018-09-28 | 2022-02-03 | Schlumberger Technology Corporation | Elastic adaptive downhole acquisition system |
-
2022
- 2022-05-18 CN CN202210540165.2A patent/CN115113281B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20220035065A1 (en) * | 2018-09-28 | 2022-02-03 | Schlumberger Technology Corporation | Elastic adaptive downhole acquisition system |
| CN110852144A (en) * | 2019-08-23 | 2020-02-28 | 中国石油大学(北京) | DTW-based intelligent stratum comparison method and system |
| CN111472765A (en) * | 2020-04-24 | 2020-07-31 | 北京月新时代科技股份有限公司 | Stratigraphic division method and device for target well |
| CN112459767A (en) * | 2020-10-30 | 2021-03-09 | 中国石油天然气集团有限公司 | Standard layer selection method and device |
| CN113236124A (en) * | 2021-05-07 | 2021-08-10 | 四川页岩气勘探开发有限责任公司 | Deep shale gas horizontal well geological guiding method |
| CN113420936A (en) * | 2021-07-01 | 2021-09-21 | 中国海洋石油集团有限公司 | Stratum comparison method based on dynamic time programming optimization algorithm |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116201535A (en) * | 2023-02-06 | 2023-06-02 | 北京月新时代科技股份有限公司 | Automatic dividing method, device and equipment for oil and gas reservoir target well sign stratum |
| CN116201535B (en) * | 2023-02-06 | 2024-02-09 | 北京月新时代科技股份有限公司 | Automatic dividing method, device and equipment for oil and gas reservoir target well sign stratum |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115113281B (en) | 2023-01-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2019201880B2 (en) | System and method for automatically correlating geologic tops | |
| US11704748B2 (en) | System and method for automatically correlating geologic tops | |
| WO2017114443A1 (en) | Method and device for determining dissolution degree of karst reservoir | |
| CN108897066B (en) | Carbonate rock crack density quantitative prediction method and device | |
| CN115113299B (en) | Intelligent marking stratum positioning and dividing method, device, equipment and storage medium | |
| CN110457501B (en) | Rock soil core sample information editing and recording system based on image recognition | |
| CN115113281B (en) | Method, device, equipment and storage medium for determining marker well and marker layer thereof | |
| CN109298464B (en) | Tight sandstone reservoir diagenetic facies logging identification method and device | |
| CN113792936A (en) | Intelligent lithology while drilling identification method, system, equipment and storage medium | |
| CN113838005A (en) | Intelligent rock fracture identification and three-dimensional reconstruction method and system based on dimension conversion | |
| CN106650446A (en) | Identification method and system of malicious program behavior, based on system call | |
| CN112800590B (en) | Grid coarsening method for machine learning-assisted two-phase flow oil reservoir random modeling | |
| CN110232208A (en) | The Reservoir Parameter Models construction method and equipment of advantage seepage flow section for identification | |
| CN110046414A (en) | The Reservoir Parameter Models construction method and equipment of advantage seepage flow section for identification | |
| CN116201535B (en) | Automatic dividing method, device and equipment for oil and gas reservoir target well sign stratum | |
| CN112081582A (en) | Prediction method, system and device for dominant channel in water-drive oil reservoir development | |
| CN115977615A (en) | Fracturing effect evaluation method, system, medium and electronic device | |
| CN110552693A (en) | layer interface identification method of induction logging curve based on deep neural network | |
| CN114612394B (en) | Concrete internal damage prediction method based on Yolo network and DIC technology | |
| CN112612935B (en) | Complete logging data acquisition method based on self-reasoning model | |
| CN112990320A (en) | Lithology classification method and device, electronic equipment and storage medium | |
| RU2760105C1 (en) | System, machine-readable media and method for core analysis by images | |
| CN109254330B (en) | Identification method and device for fractured segments of igneous rock reservoir | |
| CN113971351B (en) | Method and device for determining porosity of crack | |
| CN117911587A (en) | Digital mapping method based on geological survey |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240122 Address after: 101112, 1st Floor, Building 13, Courtyard 11, Xinyu 3rd Street, Tongzhou District, Beijing 101-566 Patentee after: Beijing Jincheng Qiuyuan Technology Co.,Ltd. Country or region after: China Address before: 100089 Courtyard 2, Dinghui North Lane, Haidian District, Beijing Patentee before: Jiang Qingqing Country or region before: China Effective date of registration: 20240122 Address after: 100089 Courtyard 2, Dinghui North Lane, Haidian District, Beijing Patentee after: Jiang Qingqing Country or region after: China Address before: Room 4001, 4th Floor, Building 27, No. 25, North Third Ring West Road, Haidian District, Beijing 100089 Patentee before: BEIJING YUEXIN TIMES TECHNOLOGY CO.,LTD. Country or region before: China |