CN113240771A - Automatic generation method and device of stratum isopachous map - Google Patents
Automatic generation method and device of stratum isopachous map Download PDFInfo
- Publication number
- CN113240771A CN113240771A CN202110582465.2A CN202110582465A CN113240771A CN 113240771 A CN113240771 A CN 113240771A CN 202110582465 A CN202110582465 A CN 202110582465A CN 113240771 A CN113240771 A CN 113240771A
- Authority
- CN
- China
- Prior art keywords
- target layer
- thickness
- dem
- hole
- stratum
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
Abstract
The invention discloses a method and a device for automatically generating a stratum contour map, wherein the method comprises the following steps: acquiring drilling information of a target layer hole and a non-target layer hole in a mapping area, and carrying out thickness value T on the target layer with the target layer holeEye shieldGenerating a thickness DEM through interpolation; calculating the average thickness of the target layer with the target layer holes; taking the absolute value of the thickness value of the target layer hole on the thickness DEM corresponding to the point position information of the target layer hole; taking the arithmetic mean of the absolute values of the average thickness of the target layer with the target layer holes and the thickness without the target holes, and multiplying the arithmetic mean by-1 to obtain the final negative value thickness of the target layer without the target layer holes; taking the thickness of a target layer, the negative value thickness of a hole without the target layer, drilling point location information corresponding to the thickness and the negative value thickness of the hole without the target layer and a geological map boundary line as discrete data, and interpolating the discrete data; carrying out contour tracing, namely generating a target stratum contour map; the invention provides an automatic generation method of a stratum isophote map with higher precision.
Description
Technical Field
The invention belongs to the field of field engineering geological analysis and evaluation, and particularly relates to a method and a device for automatically generating a formation isopachrome.
Background
The evaluation of the site engineering geological conditions is an important component of site earthquake risk evaluation and regional earthquake safety evaluation. The safety of the ground building is directly influenced by the condition of the site. Therefore, it is very important to study the geological conditions of the field project, especially in the distribution areas of important buildings or major projects. The stratum isobologram is one of geological condition distribution maps of site engineering, and refers to isobolograms of various stratums
In the prior art, the influence of formation drilling and a geological map on the formation isophotic line is not accurately considered when the isophotic map of a certain stratum is generated, when the formation is not generated at a certain drilling point is neglected, the computer forces the thickness of the stratum of the drilling hole to be 0, and interpolation is still carried out at the place without the stratum, so that the stratum still exists near the drilling hole without the stratum as the interpolation result, and the finally generated isophotic map cannot truly express the change of the thickness of the stratum; or manual delineation is time-consuming and labor-consuming, and when the stratum boundary is divided, the stratum is selected between a drill hole with a certain stratum and a drill hole without the stratum, so that the fluctuation of the stratum close to the pinch-out position is abnormal, and the precision of the generated isopachrome is reduced.
Accurate division of stratigraphic boundaries in the mapping process is a key to solving the above problems. And the location of the formation pinch-out is the formation boundary. According to geological rules, the thickness of the stratum has a direct influence on the pinch-out position, and generally, the thicker the stratum is, the farther the pinch-out distance from the hole is.
In the prior art, when the stratum boundary is divided, the influence of the thickness of the stratum is not considered when the pinch-out position of the stratum is judged, and the restriction of a geological map on the pinch-out position cannot be considered, so that the accuracy of the divided stratum boundary is obviously reduced. The geological map data directly indicate the pinch-out position of the stratum exposed out of the earth surface, and the stratum pinch-out position is a stratum boundary (namely the position of the stratum with the thickness of 0); the layered boundary line of the exposed ground surface is objective data obtained through field actual investigation, and the accuracy is high.
Disclosure of Invention
In order to overcome the technical defects, the invention provides a method and a device for automatically generating a formation contour map, and aims to provide a method for automatically generating a formation contour map with higher precision.
In order to solve the problems, the invention is realized according to the following technical scheme:
in a first aspect, the present invention provides a method for automatically generating a formation contour map, where the method includes:
acquiring drilling information of a target layer hole and a non-target layer hole in a mapping area, wherein the drilling information comprises point location information and layering information; the target layer is a stratum of which the current stratum boundary needing to be divided is in the mapping area;
calculating the thickness value T of the target layer with the target layer holesEye shieldAnd calculating the average thickness of the target layer, denoted as "T1”；
Perforating the T with the target layerEye shieldAnd corresponding point location information is used as point location data, the point location data is used as discrete data, interpolation is carried out on the point location data to generate a thickness DEM, and the thickness DEM is recorded as DEM _ TEye shield；
Slave DEM _ TEye shieldTaking the absolute value of the thickness value corresponding to the point location information of the hole without the target layer, and recording the absolute value as TMesh 2(ii) a Get T1And TMesh 2Multiplying the average value by a negative coefficient n to obtain the final target layer negative value thickness without target layer holes, and recording the thickness as TEye-friendly food(ii) a All the holes without the target layer are processed by the step, and the corresponding T is calculatedEye-friendly food；
By TEye-friendly food、TEye shieldThe corresponding point location information and the thickness value are used as point location data, the point location data and the geological map boundary line are used as discrete data, interpolation is carried out on the discrete data, and a thickness DEM is constructed and recorded as DEM _ T0Wherein the geological map boundary line assignment thickness is 0;
thickness DEM _ T0The part of the target layer with thickness larger than 0 is obtained, and the thickness DEM of the target layer is recorded as DEM _ TEyes of a user", wherein DEM _ TEyes of a userThe position with the middle thickness of 0 is the formation boundary of the target layer; for "DEM _ TEyes of a userAnd carrying out contour tracing, namely generating a target stratum contour map.
Further, before obtaining the drilling information of the target layer hole and the non-target layer hole in the imaging area, all the drilling information in the imaging area needs to be obtained, and the drilling is divided into three types including the target layer hole, the non-target layer hole and the unknown hole; therefore, the drilling information of the target layer hole and the non-target layer hole in the imaging area is obtained.
Further, the method adopted by the interpolation is a thin plate spline function method.
Further, the negative coefficient n is-1.
In a second aspect, the present invention provides an apparatus for automatically generating a contour map of a formation, including:
the information acquisition module is used for acquiring drilling information with a target layer hole and without the target layer hole in the mapping area, wherein the drilling information comprises point location information and layering information; the target layer is a stratum of which the current stratum boundary needing to be divided is in the mapping area;
a first calculating module for calculating the thickness value T of the target layer with the target layer holeEye shieldAnd calculating the average thickness of the target layer, denoted as "T1”；
A first interpolation module for calculating the thickness value T of the target layer with the target layer holeEye shieldAnd point location information is used as point location data, the point location data is used as discrete data, interpolation is carried out on the discrete data, a thickness DEM is constructed and recorded as DEM _ TEye shield；
A second calculation module for calculating the second time from DEM _ TEye shieldTaking the absolute value of the thickness value corresponding to the point location information of the hole without the target layer, and recording the absolute value as TMesh 2(ii) a Get T1And TMesh 2Multiplying the average value by a negative coefficient n to obtain the final target layer negative value thickness without target layer holes, and recording the thickness as TEye-friendly food(ii) a All the holes without the target layer are processed by the step, and the corresponding T is calculatedEye-friendly food；
A second interpolation module for interpolating with TEye-friendly food、TEye shieldThe corresponding point location information and the thickness value are used as point location data, the point location data and the geological map boundary line are used as discrete data, interpolation is carried out on the discrete data, and a thickness DEM is constructed and recorded as DEM _ T0Wherein the geological map boundary line assignment thickness is 0;
a tracking module for taking the thickness DEM _ T0The part of the target layer with thickness larger than 0 is obtained, and the thickness DEM of the target layer is recorded as DEM _ TEyes of a user"; wherein, DEM _ TEyes of a userThe position with the middle thickness of 0 is the stratum boundary; for "DEM _ TEyes of a userAnd carrying out contour tracing, namely generating a target stratum contour map.
In a third aspect, the present invention provides a terminal device, including a processor, a memory, and a program stored in the memory and configured to be executed by the processor, wherein the processor implements the method for automatically generating a stratigraphic isophote map according to the first aspect when executing the program.
Compared with the prior art, the invention has the beneficial effects that:
the invention discloses a method and a device for automatically generating a formation contour map; according to the automatic generation method of the formation isopiestic chart, drill holes are classified according to drill hole point positions and drill hole layering information, and the formation isopiechart is obtained by combining a thickness interpolation mode. The invention not only considers whether the drilling point has a target stratum, but also fully utilizes the influence of stratum thickness information on the stratum pinch-off distance, and accurately divides the stratum boundary after fully considering the constraint action of the geological map data. The method has a perfect automatic mapping process, and realizes the upgrading and conversion from the qualitative information to the formation boundary generated according to the quantitative information; the method can quickly and automatically generate the formation isopachrome without using an artificial drawing method, thereby improving the drawing efficiency and quality and reducing the construction cost.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a flow chart of a method for automatically generating a contour map of a formation according to this embodiment 1;
FIG. 2 is a schematic diagram of a pinch-out position obtained by a thickness interpolation method according to the "one-half pinch-out" rule in this embodiment 1;
FIG. 3 is a schematic diagram of pinch-out positions obtained by the prior art interpolation method without considering the formation thickness according to the "one-half pinch-out" rule in the case of multiple boreholes in this embodiment 1;
FIG. 4 is a schematic diagram comparing the effect of the "one-half pinch-out" rule of the present invention on pinch-out position obtained by the interpolation method according to the prior art without considering the thickness when multiple holes are drilled in the embodiment 1;
fig. 5 is a diagram for analyzing the thickness and gradient of soft soil in a research area generated by the method for generating a formation contour map in example 1;
FIG. 6 is a graph of soft soil thickness gradient analysis generated by the prior art in example 1;
FIG. 7 is a soft soil isobologram generated by the formation isobologram generation method in example 1
Fig. 8 is a soft soil iso-contour map generated by the prior art in example 1.
Fig. 9 is a flowchart of an automatic generation method of a soft soil contour map in embodiment 1;
fig. 10 is a program interface of the automatic soft soil contour map generation method in embodiment 1.
Detailed Description
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. It is to be understood that the embodiments described are only a few embodiments of the present application and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The flow diagrams depicted in the figures are merely illustrative and do not necessarily include all of the elements and operations/steps, nor do they necessarily have to be performed in the order depicted. For example, some operations/steps may be decomposed, combined or partially combined, so that the actual execution sequence may be changed according to the actual situation.
Example 1
In this embodiment, n strata are counted from 1 to n according to the stratum new and old sequence, wherein the stratum No. 1 is the newest, the stratum No. n is the oldest, the target stratum is represented by X, and the target stratum is a stratum of which the current stratum boundary to be divided is in the mapping area. The embodiment provides a method for generating an equal-thickness line graph of a certain stratum (target stratum X) for dividing a stratum boundary of a certain area.
Referring to fig. 1, the method for automatically generating a formation contour map according to the embodiment includes the following steps:
s10, obtaining drilling information of all holes in the mapping area, wherein the drilling information comprises point location information and layering information, the point location information is a space coordinate of the drilling hole, and the layering information comprises stratum age, lithology, burial depth and the like.
Specifically, according to the stratum information of the drill holes, all the drill holes are divided into: 3 types of target stratum holes, target stratum holes and unknown holes are formed; the target stratum hole represents that the target stratum X appears in the layered information table of the drill hole; the non-target stratum hole indicates that the target stratum X does not appear in the drilling hierarchical table, and the stratum deepest in the current drilling is older than the stratum age of the target stratum X, so that the target stratum X does not exist at the drilling point according to the standard stratum sequence (under the condition of no special construction activity, a new stratum is always above the old stratum); an unknown hole means that the hole has no target stratum X in the current hole depth, and the stratum at the deepest part of the current borehole is newer than the stratum of the target stratum X, so that whether the X stratum exists at the position of the hole can not be determined, and if drilling is continued, the target stratum X may be met, or the target stratum X may not be met, so that the hole is called as an unknown hole. In the subsequent steps, the unknown hole does not participate in any calculations.
S20, calculating the thickness value T of the target layer with the target layer holesEye shieldAnd calculating the average thickness of the target layer, denoted as "T1”。
S30, taking the thickness value T of the target layer with the target layer holesEye shieldAnd point location information is used as point location data, the point location data is used as discrete data, interpolation is carried out on the discrete data, a thickness DEM is constructed and recorded as DEM _ TEye shield。
In particular, the thickness DEM is constructed according to a predetermined interpolation method, which may be, but is not limited to, the thin-plate spline method (Wahba, 1990). The interpolation method adopts an iterative finite difference interpolation technology, has the fast calculation efficiency of a local interpolation method through optimization, and does not sacrifice the surface continuity of a global interpolation method.
S40 slave DEM _ TEye shieldTaking the absolute value of the thickness value corresponding to the point location information of the hole without the target layer, and recording the absolute value as TMesh 2(ii) a Get T1And TMesh 2Multiplying the average value by a negative coefficient n to obtain the final target layer negative value thickness without target layer holes, and recording the thickness as TEye-friendly food(ii) a All the holes without the target layer are processed by the step, and the corresponding T is calculatedEye-friendly food；
Wherein, T1And TMesh 2Are all positive numbers; multiplied by a negative factor n. The range of the negative coefficient n is (- ∞, 0). According to the principle of one-half pinch-out idea, the negative coefficient of n is defaulted to be-1', so that the final target layer negative value thickness without target layer holes is obtained.
In particular, TMesh 2Focusing on expressing the influence of local borehole information on the X-free layer thickness value, T1The two modes have respective advantages and disadvantages by focusing on expressing the influence of the drilling information in the whole imaging area on the thickness value of the X-layer-free layer. In this embodiment, the arithmetic mean is taken again for the negative thickness of the two methods, thereby balancing the influence of the local information and the global information on the value of the negative thickness.
S50, using TEye-friendly food、TEye shieldThe corresponding point location information and the thickness value are used as point location data, the point location data and the geological map boundary line are used as discrete data, interpolation is carried out on the discrete data, and a thickness DEM is constructed and recorded as DEM _ T0And the geological map boundary line is assigned with a thickness of 0.
S60, taking thickness DEM _ T0The part of the target layer with thickness larger than 0 is obtained, and the thickness DEM of the target layer is recorded as DEM _ TEyes of a user", wherein DEM _ TEyes of a userThe position with the middle thickness of 0 is the formation boundary of the target layer; for "DEM _ TEyes of a userBy contour tracing, i.e. generatingA target formation isophote map.
Based on the rule of 'one-half pinch-out', by comparing the stratum thickness-considered interpolation method of the invention with the stratum tip positions divided by the thickness-not-considered interpolation method of the prior art, it can be seen that:
when only two drill holes exist in the imaging range, the thickness of a target layer of the first drill hole is 3, and when a target layer of the second drill hole is absent, the stratum pinch-off position divided by the existing method is the middle of the two drill holes, namely the naive 'one-half pinch-off' rule; referring to fig. 2, when calculating the pinch-off position by the thickness interpolation method based on the rule of 'one-half pinch-off' of the method of the present invention, a 'negative thickness', namely '-3 m', is firstly taken for the second hole, and then the pinch-off position is simulated by the thickness interpolation method, wherein the pinch-off position is also the middle of the two holes.
However, when a plurality of drill holes exist, the thickness of each drill hole is not uniform, and compared with the conventional forced half pinch-out mode without considering the thickness of the stratum and the effect of calculating the pinch-out position by the thickness interpolation mode of the method disclosed by the invention: referring to fig. 3, the difference between the thicknesses of the left and right holes is obvious, but under the existing 'one-half pinch-out' rule without considering the thickness, the pinch-out distances of the left and right holes are the same and are both one-half pitch. The simplest geological rule that the larger the stratum thickness is, the farther the pinch-off distance from the hole is not met. This can result in the formation heave varying in abnormally large amplitude (left-hand formation with a large slope, but right-hand formation with a small slope), and the map rationality is significantly questionable. Referring to fig. 4, the method for dividing the stratum pinch-off position provided by the invention considers the thickness of the stratum, obviously improves the consistency and uniformity of the fluctuation of the stratum, completely accords with the geological rule that the larger the thickness of the stratum is, the farther the pinch-off distance from the hole is, and has more reasonable stratum boundary and stratum thickness variation trend.
And when the thickness coefficient of the hole without the target layer is minus 1, the mapping area integrally ensures that the stratum coverage area occupies half of the whole mapping range and conforms to the 'half pinch-out' rule.
Comparing the automatic mapping method based on accurate division of the stratum boundary with the ordinary forced one-half pinch-out mapping method, the difference between the two methods is as follows: the forced half pinch-out is qualitative division, namely, the formation boundary is divided only based on the qualitative information of existence and nonexistence of the target formation, but the thickness-based formation boundary of the method not only considers the existence of the formation, but also fully utilizes the formation thickness information, and realizes the purposes of quantitatively, accurately and reasonably dividing the formation boundary and generating a formation isopiegram based on the formation thickness.
As a more specific embodiment, 5km somewhere in Guangzhou2For a research area, 8 borehole stratum information is utilized to automatically generate a soft soil isobologram, and the method is specifically described as an embodiment, and a method flow chart refers to fig. 9.
1. And acquiring point location information and layering information of the eight drill holes.
2. According to the drill hole layering information, dividing all drill holes into: soft soil holes, no soft soil holes and unknown holes. The specific classification method comprises the following steps: inquiring from the ground surface downwards, wherein a soft soil hole indicates that soft soil appears in the stratum information of the drill hole;
the fact that no soft soil hole exists means that no soft soil exists in the surface of the stratum of the drilled hole, and the stratum at the deepest part of the drilled hole is older than the stratum of the soft soil, so that the new stratum is always above the old stratum under normal conditions, and therefore the drilled hole is determined to have no soft soil;
the unknown hole means that the hole has no soft soil when viewed according to the current hole depth, the stratum at the deepest part of the current drill hole is newer than the soft soil stratum, if drilling is continued, the soft soil may occur or may not occur, and therefore whether the hole has the soft soil or not cannot be determined, and the hole is called as the unknown hole. In the subsequent steps of the method of the invention, the unknown pore does not participate in any calculations.
3. Calculating the thickness of the soft soil in each hole and recording as TSoft soilAnd calculating the average thickness of the soft soil with the soft soil holes and recording as T1”。
4. With point location information of soft soil hole and corresponding soft soil thickness TSoft soilAs point location data, all point location data are used as discrete data, and the discrete data are processedThe data is interpolated by a discretized thin-plate spline method of iterative finite difference interpolation to generate a thickness DEM which is recorded as DEM _ TSoft soil”。
5. Point location slave DEM _ T without soft soil holeSoft soilThe absolute value of the corresponding thickness value is taken and recorded as TSoft soil 2。
6. Taking the thickness T1And TSoft soil 2Average thickness (T) of1And TSoft soil 2Positive numbers) and then multiplied by a negative coefficient n (n is-1), so as to obtain the final soft soil negative value thickness without soft soil holes, which is recorded as TSoft soil. All the holes without soft soil are processed by the step, and the corresponding T is calculatedSoft soil；
7. Will TSoft soil、TSoft soilThe corresponding point location information and the thickness value are used as point location data, the point location data and a geological map boundary line (the assigned thickness is 0) are used as discrete data, interpolation is carried out on the discrete data to generate a thickness DEM, and the thickness DEM is recorded as DEM _ T0。
8. Thickness DEM _ T0The part of the soft soil DEM which is larger than 0 is obtained, and is marked as DEM _ TSoft soil”。
9. For "DEM _ TSoft soil' carrying out contour tracing, namely generating a soft soil contour map.
And programming and integrating the steps to form a functional module for automatically generating the soft soil isobologram. And selecting the drill holes participating in the drawing during drawing, clicking a 'field seismic engineering geological condition' → 'isophotic line' → 'soft soil' in a program interface (figure 10), entering a 'soft soil isophotic line' interface, clicking a 'create …' menu at the lower right corner in the figure, and adjusting all entries to generate a soft soil isophotic graph.
Fig. 5 shows an analysis diagram of the soft soil thickness and gradient in the research area generated by the formation contour map generation method in this embodiment 1, where an average gradient of 0.89 degrees and a standard deviation of the gradient of 0.38 can be found, and fig. 6 shows an analysis diagram of the soft soil thickness and gradient generated in the prior art, where an average gradient of 0.82 degrees and a standard deviation of the gradient of 0.517 can be found. Compared with the two methods, the automatic generation method of the formation isophotic graph in the embodiment fully considers the influence of the formation thickness, the gradient difference is obviously smaller than that of the existing method, the consistency and uniformity of the formation fluctuation are obviously improved, and the geological rule that the larger the formation thickness is, the farther the pinch-off distance from the hole is better met.
Fig. 7 shows a soft soil contour map generated by using the formation contour map generation method in this embodiment 1, fig. 8 shows a soft soil contour map generated by the prior art, the soft soil contour map generated by comparing the two methods is compared, the contour map generated by dividing the formation boundary based on the quantitative formation thickness fully considers the influence of the formation thickness, and the divided soft soil boundary is highly matched with the geological map, which indicates that the contour map generated by the method of the present invention is more accurate.
In this embodiment, the dividing of the formation boundary in consideration of the formation thickness refers to quantitatively dividing a boundary of a certain formation by using the formation thickness according to the information of the borehole formation, and the accuracy of the dividing is a key link for automatically generating a formation isopieogram.
Example 2
An embodiment of the present invention provides an automatic generation apparatus for a formation isophote map, including:
the information acquisition module is used for acquiring drilling information with a target layer hole and without the target layer hole in the mapping area, wherein the drilling information comprises point location information and layering information; the target layer is a stratum of which the current stratum boundary needing to be divided is in the mapping area;
a first calculating module for calculating the thickness value T of the target layer with the target layer holeEye shieldAnd calculating the average thickness of the target layer, denoted as "T1”；
A first interpolation module for using the point location information of the target layer hole and the corresponding thickness value TEye shieldAnd taking the point location data as discrete data as point location data, interpolating the point location data to generate a thickness DEM, and recording the thickness DEM as DEM _ TEye shield；
A second calculation module for calculating the second time from DEM _ TEye shieldTaking the absolute value of the thickness value corresponding to the point location information of the hole without the target layer, and recording the absolute value as TMesh 2(ii) a Get T1And TMesh 2The average value of the target layer is multiplied by a negative coefficient n to obtain the final target layer negative value without target layer holesThickness, noted as TEye-friendly food(ii) a All the holes without the target layer are processed by the step, and the corresponding T is calculatedEye-friendly food；
A second interpolation module for interpolating with TEye-friendly food、TEye shieldThe corresponding point location information and the thickness value are used as point location data, the point location data and the geological map boundary line are used as discrete data, interpolation is carried out on the discrete data, and a thickness DEM is constructed and recorded as DEM _ T0Wherein the geological map boundary line assignment thickness is 0;
a tracking module for taking the thickness DEM _ T0The part of the target layer with thickness larger than 0 is obtained, and the thickness DEM of the target layer is recorded as DEM _ TEyes of a user", wherein DEM _ TEyes of a userThe position with the middle thickness of 0 is the stratum boundary; for "DEM _ TEyes of a userAnd carrying out contour tracing, namely generating a target stratum contour map.
Example 3
The embodiment of the invention provides a terminal device, which comprises a processor and a memory, wherein the memory is used for storing a computer program; the processor is configured to execute the computer program and implement the method for automatically generating a formation contour map provided in embodiment 1 of the present invention when the computer program is executed.
It should be understood that the Processor may be a Central Processing Unit (CPU), and the Processor may be other general purpose processors, Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field Programmable Gate Arrays (FPGAs) or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components, etc. Wherein a general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
While the invention has been described with reference to specific embodiments, the invention is not limited thereto, and various equivalent modifications and substitutions can be easily made by those skilled in the art within the technical scope of the invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (6)
1. A method for automatically generating a formation isophote map, the method comprising:
acquiring drilling information of a target layer hole and a non-target layer hole in a mapping area, wherein the drilling information comprises point location information and layering information; the target layer is a stratum of which the current stratum boundary needing to be divided is in the mapping area;
calculating the thickness value T of the target layer with the target layer holesEye shieldAnd calculating the average thickness of the target layer, denoted as "T1”；
Perforating the T with the target layerEye shieldAnd corresponding point location information is used as point location data, the point location data is used as discrete data, interpolation is carried out on the point location data to generate a thickness DEM, and the thickness DEM is recorded as DEM _ TEye shield；
Slave DEM _ TEye shieldTaking the absolute value of the thickness value corresponding to the point location information of the hole without the target layer, and recording the absolute value as TMesh 2(ii) a Get T1And TMesh 2Multiplying the average value by a negative coefficient n to obtain the final target layer negative value thickness without target layer holes, and recording the thickness as TEye-friendly food(ii) a All the holes without the target layer are processed by the step, and the corresponding T is calculatedEye-friendly food；
By TEye-friendly food、TEye shieldThe corresponding point location information and the thickness value are used as point location data, the point location data and the geological map boundary line are used as discrete data, interpolation is carried out on the discrete data, and a thickness DEM is constructed and recorded as DEM _ T0Wherein the geological map boundary line assignment thickness is 0;
thickness DEM _ T0The part of the target layer with thickness larger than 0 is obtained, and the thickness DEM of the target layer is recorded as DEM _ TEyes of a user", wherein DEM _ TEyes of a userThe position with the middle thickness of 0 is the formation boundary of the target layer; for "DEM _ TEyes of a userAnd carrying out contour tracing, namely generating a target stratum contour map.
2. The method according to claim 1, wherein before obtaining the drilling information of the target layer hole and the non-target layer hole in the mapping area, all the drilling information in the mapping area needs to be obtained, and the drilling is divided into three types, including the target layer hole, the non-target layer hole and the unknown hole; therefore, the drilling information of the target layer hole and the non-target layer hole in the imaging area is obtained.
3. The method of claim 1, wherein the interpolation is performed by thin-plate spline.
4. The method of claim 1, wherein the negative coefficient n is-1.
5. An apparatus for automatically generating a contour map of a formation, comprising:
the information acquisition module is used for acquiring drilling information with a target layer hole and without the target layer hole in the mapping area, wherein the drilling information comprises point location information and layering information; the target layer is a stratum of which the current stratum boundary needing to be divided is in the mapping area;
a first calculating module for calculating the thickness value T of the target layer with the target layer holeEye shieldAnd calculating the average thickness of the target layer, denoted as "T1”；
A first interpolation module for calculating the thickness value T of the target layer with the target layer holeEye shieldAnd point location information is used as point location data, the point location data is used as discrete data, interpolation is carried out on the discrete data, a thickness DEM is constructed and recorded as DEM _ TEye shield；
A second calculation module for calculating the second time from DEM _ TEye shieldTaking the absolute value of the thickness value corresponding to the point location information of the hole without the target layer, and recording the absolute value as TMesh 2(ii) a Get T1And TMesh 2Multiplying the average value by a negative coefficient n to obtain the final target layer negative value thickness without target layer holes, and recording the thickness as TEye-friendly food(ii) a All the holes without the target layer are processed by the step, and the corresponding T is calculatedEye-friendly food；
A second interpolation module for interpolating with TEye-friendly food、TEye shieldThe corresponding point location information and the thickness value are used as point location data, the point location data and the geological map boundary line are used as discrete data, interpolation is carried out on the discrete data, and a thickness DEM is constructed and recorded as DEM _ T0Wherein the geological map boundary line assignment thickness is 0;
a tracking module for taking the thickness DEM _ T0The part of the target layer with thickness larger than 0 is obtained, and the thickness DEM of the target layer is recorded as DEM _ TEyes of a user"; wherein, DEM _ TEyes of a userThe position with the middle thickness of 0 is the stratum boundary; for "DEM _ TEyes of a userAnd carrying out contour tracing, namely generating a target stratum contour map.
6. A terminal device comprising a processor, a memory, and a program stored in the memory and configured to be executed by the processor, the processor implementing the method of automatically generating a stratigraphic isophote map according to any of claims 1 to 5 when executing the program.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110582465.2A CN113240771B (en) | 2021-05-27 | 2021-05-27 | Automatic generation method and device of stratum isopachous map |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110582465.2A CN113240771B (en) | 2021-05-27 | 2021-05-27 | Automatic generation method and device of stratum isopachous map |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113240771A true CN113240771A (en) | 2021-08-10 |
| CN113240771B CN113240771B (en) | 2021-12-07 |
Family
ID=77139084
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110582465.2A Active CN113240771B (en) | 2021-05-27 | 2021-05-27 | Automatic generation method and device of stratum isopachous map |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113240771B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040181341A1 (en) * | 2003-03-10 | 2004-09-16 | Neff Dennis B. | Contour-based attribute seismic mapping |
| CN102760307A (en) * | 2011-04-29 | 2012-10-31 | 黄桂芝 | Design method of three-dimensional simulation system for mining |
| CN108986213A (en) * | 2018-07-02 | 2018-12-11 | 南京师范大学 | A kind of three dimensional contour line method based on stacking technology |
| CN109979011A (en) * | 2019-03-22 | 2019-07-05 | 李潇 | Plains region three-dimensional geological model building method based on multi-source heterogeneous data |
| CN111472763A (en) * | 2019-12-30 | 2020-07-31 | 北京国双科技有限公司 | Stratum thickness prediction method and device |
| CN111950046A (en) * | 2020-07-03 | 2020-11-17 | 成都理工大学 | Drilling data model construction method based on BIM |
| CN112381937A (en) * | 2020-11-03 | 2021-02-19 | 广东省地质调查院 | Multi-source geological data coupling modeling method based on drilling and complex geological profile |
| CN112529981A (en) * | 2020-11-24 | 2021-03-19 | 中国煤炭地质总局地球物理勘探研究院 | Depth domain coal seam floor contour line mapping method |
| CN112837413A (en) * | 2021-02-08 | 2021-05-25 | 南京师范大学 | Geological drilling-oriented virtual stratum inference method and device |
-
2021
- 2021-05-27 CN CN202110582465.2A patent/CN113240771B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040181341A1 (en) * | 2003-03-10 | 2004-09-16 | Neff Dennis B. | Contour-based attribute seismic mapping |
| CN102760307A (en) * | 2011-04-29 | 2012-10-31 | 黄桂芝 | Design method of three-dimensional simulation system for mining |
| CN108986213A (en) * | 2018-07-02 | 2018-12-11 | 南京师范大学 | A kind of three dimensional contour line method based on stacking technology |
| CN109979011A (en) * | 2019-03-22 | 2019-07-05 | 李潇 | Plains region three-dimensional geological model building method based on multi-source heterogeneous data |
| CN111472763A (en) * | 2019-12-30 | 2020-07-31 | 北京国双科技有限公司 | Stratum thickness prediction method and device |
| CN111950046A (en) * | 2020-07-03 | 2020-11-17 | 成都理工大学 | Drilling data model construction method based on BIM |
| CN112381937A (en) * | 2020-11-03 | 2021-02-19 | 广东省地质调查院 | Multi-source geological data coupling modeling method based on drilling and complex geological profile |
| CN112529981A (en) * | 2020-11-24 | 2021-03-19 | 中国煤炭地质总局地球物理勘探研究院 | Depth domain coal seam floor contour line mapping method |
| CN112837413A (en) * | 2021-02-08 | 2021-05-25 | 南京师范大学 | Geological drilling-oriented virtual stratum inference method and device |
Non-Patent Citations (6)
| Title |
|---|
| 孟凡刚等: "基于钻孔勘探数据的煤层厚度分布与构造预判", 《煤炭科学技术》 * |
| 曾梦秋: "基于钻孔数据的三维地层模型建立及工程应用研究", 《中国优秀硕士学位论文全文数据库（电子期刊）工程科技I辑》 * |
| 李程等: "基于GOCAD的钻孔数据建模方法与实例研究", 《江西科学》 * |
| 王金艳等: "郯庐断裂带江苏段新生界三维地质构造建模", 《地震学报》 * |
| 郑先昌等: "基于GIS与稀疏钻孔数据的工程地质等值线绘制方法", 《工程勘察》 * |
| 陈小芳等: "基于支持向量机的地震灾害自动成图方法", 《灾害学》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113240771B (en) | 2021-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2657715C (en) | Providing a simplified subterranean model | |
| US20180321127A1 (en) | Method for Providing a Numerical Model of a Sample of Rock | |
| NO328482B1 (en) | Method, product and system for automated, non-artificial, three-dimensional fault surface-based horizon modeling of the subsurface | |
| US20180119523A1 (en) | Oilfield Reservoir Saturation and Permeability Modeling | |
| US11066905B2 (en) | Oilfield reservoir saturation and permeability modeling | |
| CN113240771B (en) | Automatic generation method and device of stratum isopachous map | |
| US20220206175A1 (en) | Training a machine learning system using hard and soft constraints | |
| CN110765665B (en) | Dynamic modeling method and system for geography | |
| CN113269879B (en) | Automatic construction method and device for three-dimensional model of geologic body | |
| CN109188522B (en) | Velocity field construction method and device | |
| Rodriguez | Integration of topographic and bathymetric digital elevation model using ArcGIS interpolation methods: a case study of the Klamath River Estuary | |
| CN113269852B (en) | Stratum partition map and stratum isobath mapping method and device | |
| US11209572B2 (en) | Meshless and mesh-based technique for modeling subterranean volumes | |
| CN112965142B (en) | Method, device, equipment and storage medium for generating karst ancient political view | |
| Omeragic et al. | Workflow to automatically update geological models during well placement with high angle and horizontal well log interpretation results | |
| Vaes et al. | APWP-online. org: a global reference database and open-source tools for calculating apparent polar wander paths and relative paleomagnetic displacements | |
| CN113219534B (en) | Pre-stack depth migration speed quality control method and device, medium and electronic equipment | |
| CN113495303B (en) | Method and device for recovering denudation quantity | |
| CN113219533B (en) | Pre-stack time migration velocity modeling method and device, medium and electronic equipment | |
| CN112346117B (en) | Reservoir characteristic prediction method and device based on seismic attribute fusion | |
| Araújo et al. | Geostatistical simulation with heterotopic soft data without the LMC | |
| Zhao et al. | A Method of Simulating Regional Groundwater level Distribution Based on Digital Elevation Model | |
| Anisimov et al. | Method of Multivariate Seismic Tomography Based on Stochastic Analysis of Input Data Reliability | |
| CN117828538A (en) | Multi-source information comprehensive analysis method and system based on weight distribution | |
| CN117930327A (en) | Excitation well depth extraction method and device and electronic equipment |
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 |