A boundary recognition method, apparatus, device, and medium based on gravity data

Abstract

This application provides a boundary identification method, apparatus, device, and medium based on gravity data, relating to the field of geophysical exploration technology. The boundary identification method includes: acquiring first information; wherein the first information includes gravity data and coordinates of ground observation points; using the first information to obtain the positions of upper and lower interface feature points of a geological body in underground space; and obtaining boundary information of the geological body based on the positions of the upper and lower interface feature points, wherein the boundary information of the geological body includes the extension position of the geological body interface. The technical solution provided in this application can identify the boundary of a geological body by using the positions of the upper and lower interface feature points in underground space, achieving the identification of inclined boundaries, and is not limited by the assumption of the boundary's vertical angle, effectively improving the reliability of the geological body boundary identification results.

Description

Technical Field

[0001] This application relates to the field of geophysical exploration technology, and in particular to a boundary identification method, apparatus, equipment and medium based on gravity data. Background Technology

[0002] Boundary identification of subsurface geological bodies essentially involves using surface or near-surface observation data to invert the physical interfaces formed by density differences underground, providing a basis for geological interpretation, target area delineation, or risk assessment. Gravity exploration technology, due to its inherent sensitivity to density differences and its ability to collect data across multiple scenarios including surface, aerial, and marine environments, has become an important means of boundary identification.

[0003] Currently, boundary identification methods based on gravity data all assume vertical boundaries. However, in actual geological conditions, geological body boundaries are mostly non-vertical, leading to discrepancies between the identified boundary locations and the actual boundary locations under this assumption. Therefore, improving the reliability of boundary identification based on gravity data has become an urgent problem to be solved. Summary of the Invention

[0004] In view of this, the purpose of this application is to provide a boundary identification method, device, equipment and medium based on gravity data, which can identify the boundary of a geological body by the location of the upper interface feature points and the lower interface feature points of the geological body in underground space, realize the identification of inclined boundaries, and not be limited by the premise assumption of the vertical angle of the boundary, effectively improve the reliability of the geological body boundary identification results, and thus provide more solid evidence information for subsequent work such as geological interpretation.

[0005] This application mainly includes the following aspects:

[0006] In a first aspect, this application provides a boundary recognition method based on gravity data, the boundary recognition method comprising:

[0007] Obtain first information; wherein, the first information includes gravity data of the ground observation point and the coordinates of the ground observation point;

[0008] The location of the upper interface feature point and the lower interface feature point of the geological body in the underground space are obtained using the first information.

[0009] Based on the locations of the upper and lower interface feature points of the geological body, the boundary information of the geological body is obtained; wherein, the boundary information of the geological body includes the extension location of the interface of the geological body.

[0010] Furthermore, the step of obtaining the location of the upper interface feature point and the lower interface feature point of the geological body in the underground space using the first information includes:

[0011] Using the first information, gravity anomaly data corresponding to any depth in the underground space are obtained;

[0012] Obtain the scaling function corresponding to any depth in the underground space;

[0013] Based on the gravity anomaly data and the scaling function, the density distribution intensity at the first arbitrary point in the underground space is obtained;

[0014] In a first set consisting of the density distribution intensity at all first arbitrary points in the underground space, the positions of the upper interface feature points and the lower interface feature points of the geological body in the underground space are obtained through the extreme points in the first set.

[0015] Furthermore, the step of obtaining the boundary information of the geological body based on the positions of the upper interface feature points and the lower interface feature points of the geological body includes:

[0016] Based on the location of the upper interface feature point and the location of the lower interface feature point of the geological body, the axial constraint quantity at the first arbitrary point in the underground space is determined.

[0017] Using the axial constraints, a model weighting function is constructed at the first arbitrary point in the underground space;

[0018] Get the data weighting function;

[0019] Based on the model weighting function and the data weighting function, a physical property parameter model is obtained;

[0020] The linear information in the physical property parameter model is extracted using the finite difference method to obtain the boundary information of the geological body.

[0021] Furthermore, the step of obtaining gravity anomaly data corresponding to any depth in the underground space using the first information includes:

[0022] The gravity anomaly spectrum of the ground observation point is obtained by using the Fourier transform method on the gravity data of the ground observation point in the first information.

[0023] Based on the gravity anomaly spectrum of the ground observation points, determine the gravity anomaly spectrum at any depth in the underground space;

[0024] The gravity anomaly spectrum at any depth in the underground space is subjected to Fourier transform to obtain gravity anomaly data corresponding to any depth in the underground space.

[0025] Furthermore, the step of constructing the model weighting function at the first arbitrary point in the underground space using the axis constraints includes:

[0026] The axial constraint quantities at all first arbitrary points in the underground space are normalized to obtain a normalized weighted matrix;

[0027] Perform diagonal matrix operations on the normalized weighted matrix to obtain the first matrix;

[0028] Based on the first matrix, the model weighting function at the first arbitrary point in the underground space is obtained.

[0029] Furthermore, the data weighting function is obtained through the following steps:

[0030] Based on the coordinates of the first arbitrary point in the underground space and the coordinates of the ground observation point, a first forward operator is constructed corresponding to the first arbitrary point of the second arbitrary point in the ground observation point.

[0031] Based on the first forward operator, the second forward operator corresponding to the second arbitrary point is obtained;

[0032] The data weighting function is obtained through the second forward operator.

[0033] Secondly, this application also provides a boundary recognition device based on gravity data, the boundary recognition device comprising:

[0034] An acquisition module is used to acquire first information; wherein, the first information includes gravity data of a ground observation point and the coordinates of the ground observation point;

[0035] The processing module is used to obtain the location of the upper interface feature point and the lower interface feature point of the geological body in the underground space using the first information;

[0036] The identification module is used to obtain the boundary information of the geological body based on the location of the upper interface feature points and the location of the lower interface feature points of the geological body; wherein, the boundary information of the geological body includes the extension position of the interface of the geological body.

[0037] Furthermore, the processing module is specifically used for:

[0038] Using the first information, gravity anomaly data corresponding to any depth in the underground space are obtained;

[0039] Obtain the scaling function corresponding to any depth in the underground space;

[0040] Based on the gravity anomaly data and the scaling function, the density distribution intensity at the first arbitrary point in the underground space is obtained;

[0041] In a first set consisting of the density distribution intensity at all first arbitrary points in the underground space, the positions of the upper interface feature points and the lower interface feature points of the geological body in the underground space are obtained through the extreme points in the first set.

[0042] Thirdly, this application also provides an electronic device, including: a processor, a memory, and a bus, wherein the memory stores machine-readable instructions executable by the processor, and when the electronic device is running, the processor communicates with the memory via the bus, and when the machine-readable instructions are executed by the processor, the steps of the boundary recognition method based on gravity data as described above are performed.

[0043] Fourthly, this application also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, performs the steps of the boundary recognition method based on gravity data as described above.

[0044] This application provides a boundary identification method, apparatus, device, and medium based on gravity data. The boundary identification method includes: acquiring first information; wherein the first information includes gravity data of a ground observation point and the coordinates of the ground observation point; using the first information to obtain the positions of upper interface feature points and lower interface feature points of a geological body in underground space; and obtaining boundary information of the geological body based on the positions of the upper interface feature points and the lower interface feature points of the geological body; wherein the boundary information of the geological body includes the extension position of the interface of the geological body.

[0045] Thus, the technical solution provided in this application can identify the boundary of a geological body by the location of the upper interface feature points and the lower interface feature points in the underground space, thereby achieving the identification of inclined boundaries. It is not limited by the premise assumption of the vertical angle of the boundary, effectively improving the reliability of the geological body boundary identification results, and thus providing more solid evidence for subsequent work such as geological interpretation.

[0046] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, preferred embodiments are described below in detail with reference to the accompanying drawings. Attached Figure Description

[0047] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0048] Figure 1 A flowchart of a boundary recognition method based on gravity data provided in an embodiment of this application is shown;

[0049] Figure 2 A schematic diagram of an underground space grid structure provided in an embodiment of this application is shown;

[0050] Figure 3 This illustration shows one of the schematic diagrams of boundary information provided in an embodiment of this application;

[0051] Figure 4 This illustrates a second schematic diagram of boundary information provided in an embodiment of this application;

[0052] Figure 5 This illustration shows a third schematic diagram of boundary information provided in an embodiment of this application;

[0053] Figure 6 This illustration shows a fourth schematic diagram of boundary information provided in an embodiment of this application;

[0054] Figure 7 This illustration shows a structural schematic diagram of a boundary recognition device based on gravity data provided in an embodiment of this application;

[0055] Figure 8 A schematic diagram of the structure of an electronic device provided in an embodiment of this application is shown. Detailed Implementation

[0056] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. It should be understood that the drawings in this application are for illustrative and descriptive purposes only and are not intended to limit the scope of protection of this application. Furthermore, it should be understood that the schematic drawings are not drawn to scale. The flowcharts used in this application illustrate operations implemented according to some embodiments of this application. It should be understood that the operations in the flowcharts may not be implemented in sequence, and steps without logical contextual relationships may be reversed or implemented simultaneously. In addition, those skilled in the art, guided by the content of this application, may add one or more other operations to the flowcharts, or remove one or more operations from the flowcharts.

[0057] Furthermore, the described embodiments are merely some, not all, of the embodiments of this application. The components of the embodiments of this application described and illustrated herein can typically be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of this application provided in the accompanying drawings is not intended to limit the scope of the claimed application, but merely to illustrate selected embodiments of the application. All other embodiments obtained by those skilled in the art based on the embodiments of this application without inventive effort are within the scope of protection of this application.

[0058] In order to enable those skilled in the art to use the content of this application, and in combination with the specific application scenario of "boundary recognition based on gravity data", the following implementation method is given. For those skilled in the art, the general principles defined herein can be applied to other embodiments and application scenarios without departing from the spirit and scope of this application.

[0059] The methods, apparatus, electronic devices, or computer-readable storage media described in this application can be applied to any scenario requiring boundary recognition based on gravity data. This application does not limit specific application scenarios. Any solution using a boundary recognition method, apparatus, electronic device, or storage medium based on gravity data provided in this application is within the protection scope of this application.

[0060] It is worth noting that boundary identification methods based on gravity data have been widely developed. Compared with other geophysical methods, their advantages are: (1) Gravity exploration methods have advantages such as economy and strong applicability. The collected data is relatively stable and less affected by the external environment, and has high reliability. Therefore, it is one of the most mature geophysical exploration methods; (2) One of the most typical advantages of gravity data is that it has high lateral resolution, which means that it has a good identification effect on the horizontal projection position of geological body interfaces, contact surfaces, faults, etc. Therefore, boundary identification methods are also one of the most important research directions of gravity data processing methods. This branch of technology has been widely studied, developed and applied, and has achieved good application results.

[0061] However, the main problem with current boundary identification / enhancement technologies based on gravity data is that they all assume vertical boundaries. In actual geological conditions, geological body boundaries (including contact surfaces, faults, and lithological interfaces) are mostly non-vertical. Under this assumption, the calculated boundary positions will deviate from the actual positions, limiting the reliability of the inversion calculation results in depicting the actual geological body. This affects both the reliability of the inversion results and their wide applicability. Therefore, how to fully leverage the lateral resolution advantage of gravity data to further effectively depict inclined boundary information, thereby broadening the applicability of gravity data-based boundary identification and improving the reliability of boundary identification results, has always been a pressing issue to be addressed in boundary identification / enhancement technologies based on gravity data.

[0062] Based on this, this application proposes a boundary identification method, apparatus, device, and medium based on gravity data. The boundary identification method includes: acquiring first information; wherein the first information includes gravity data of a ground observation point and the coordinates of the ground observation point; using the first information to obtain the positions of the upper interface feature points and the lower interface feature points of a geological body in the underground space; and obtaining the boundary information of the geological body based on the positions of the upper interface feature points and the lower interface feature points of the geological body; wherein the boundary information of the geological body includes the extension position of the interface of the geological body.

[0063] Thus, the technical solution provided in this application can identify the boundary of a geological body by the location of the upper interface feature points and the lower interface feature points in the underground space, thereby achieving the identification of inclined boundaries. It is not limited by the premise assumption of the vertical angle of the boundary, effectively improving the reliability of the geological body boundary identification results, and thus providing more solid evidence for subsequent work such as geological interpretation.

[0064] To facilitate understanding of this application, the technical solutions provided in this application will be described in detail below with reference to specific embodiments.

[0065] Please see Figure 1 , Figure 1 A flowchart of a boundary recognition method based on gravity data provided in an embodiment of this application is shown, as follows: Figure 1 As shown, the boundary recognition method includes:

[0066] S101, Obtain first information.

[0067] For example, the first information can be in the form of a file or in the form of data. Users can upload the first information in the form of a file or input the first information in the form of data.

[0068] The first piece of information includes gravity data from ground observation points. Coordinates of ground observation points and the number of ground observation points .

[0069] S102. Using the first information, obtain the location of the upper interface feature point and the lower interface feature point of the geological body in the underground space.

[0070] For example, a grid coordinate system can be established with the ground observation surface as the X-axis and the direction perpendicular to the ground observation surface as the Z-axis to construct the grid structure of the underground space. The grid structure of the underground space can be used to determine the location of the upper interface feature points and the lower interface feature points of the geological body in the underground space.

[0071] For example, before constructing the grid structure of the underground space, users can upload files containing calculation parameters or input calculation parameters, including underground space partitioning parameters (such as the grid size in the X direction). Z-axis mesh size Number of mesh cells in the X direction and the number of mesh elements in the Z direction ), the ratio of the two orders and And the coefficients of the inversion process (e.g., maximum number of iterations) and maximum allowable fitting error wait).

[0072] For example, please refer to Figure 2 , Figure 2 This illustration shows a schematic diagram of an underground space grid structure provided in an embodiment of this application. A computer can read the underground space partitioning parameters from the first information and calculation parameters to construct the underground space grid structure, such as... Figure 2 As shown, the underground space is divided into multiple uniformly sized grid cells, each with a size of [missing information]. The rectangle is of size [size missing], and the triangular points on the X-axis represent the observation points. Indicates the observation point coordinates The x-coordinate on the observation surface is Gravity observations (gravity data) at the location, with dots representing the subdivision points of the underground space grid structure. Figure 2 The “…” in the text indicates omitted subdivision points. Indicates the partition point coordinates Indicates the depth of underground space.

[0073] In some embodiments, the step of obtaining the location of the upper interface feature point and the lower interface feature point of a geological body in an underground space using the first information includes:

[0074] S1021. Using the first information, obtain gravity anomaly data corresponding to any depth in the underground space.

[0075] In some embodiments, the step of obtaining gravity anomaly data corresponding to any depth in the underground space using the first information includes:

[0076] 1) Using the Fourier transform method, the gravity anomaly spectrum of the ground observation points is obtained from the gravity data of the ground observation points in the first information.

[0077] For example, using gravity data The Fourier transform method is used to obtain the results for any observation point. Gravity anomaly spectrum The formula is as follows:

[0078]

[0079] in, For wave number, , This is the symbol for the Fourier transform.

[0080] 2) Determine the gravity anomaly spectrum at any depth in the underground space based on the gravity anomaly spectrum of ground observation points.

[0081] For example, to obtain the corresponding underground space at any depth The spectrum of gravity anomalies The formula is as follows:

[0082]

[0083] in, For underground space Wavenumber at each partition point along the X-axis, .

[0084] 3) The Fourier transform method is used to obtain the gravity anomaly data corresponding to any depth in the underground space by applying the gravity anomaly spectrum at any depth.

[0085] For example, the inverse Fourier transform can be used to obtain the depth of any object. Corresponding gravity anomaly data The formula is as follows:

[0086]

[0087] in, This is the symbol for the inverse Fourier transform.

[0088] In this way, a set of gravity anomaly data at any point in the underground space can be obtained. The details are as follows:

[0089]

[0090] S1022. Obtain the scale function corresponding to any depth in the underground space.

[0091] For example, the computer reads the ratio of two orders from the calculation parameters. and Determine any depth of underground space Corresponding scaling function The specific formula is as follows:

[0092]

[0093] in:

[0094]

[0095]

[0096] here, Indicates the gradient sign. express Pointed order ratio, express Pointed Order ratio.

[0097] S1023. Based on gravity anomaly data and scale function, obtain the density distribution intensity at the first arbitrary point in the underground space.

[0098] For example, the product of gravity anomaly data and the corresponding scaling function can be used to determine the density distribution intensity at the corresponding subdivision point (the first arbitrary point) in the underground space. The specific formula is as follows:

[0099]

[0100] S1024. In the first set consisting of the density distribution intensity at all first arbitrary points in the underground space, the positions of the upper interface feature points and the lower interface feature points of the geological body in the underground space are obtained through the extreme points in the first set.

[0101] For example, the first set here can be the set of density distribution intensities obtained by normalizing the density distribution intensities at all first arbitrary points in the underground space; firstly, the set of density distribution intensities formed by arbitrary dividing points (i.e., the first arbitrary points) in the underground space. It can be represented as follows:

[0102]

[0103] Then, the density distribution intensity at the first arbitrary point in the underground space can be normalized to obtain the normalized density distribution intensity. The specific formula is as follows:

[0104]

[0105] in, It is the first set The maximum value in, It is the first set The minimum value in.

[0106] Finally, we can obtain the set of normalized density distribution intensities at all first arbitrary points in the underground space (i.e., the first set). , means as follows:

[0107]

[0108] For example, determine the first set Maximum point in Minimum point Construct the characteristic formula of the inclined surface using the maximum and minimum points. The specific formula is as follows:

[0109]

[0110] in, It is a maximum point x-coordinate It is a maximum point The ordinate, It is a local minimum point x-coordinate It is a local minimum point The ordinate.

[0111] For example, the maximum point can be determined as the location of the upper interface feature point of the geological body in the underground space, and the minimum point can be determined as the location of the lower interface feature point of the geological body in the underground space.

[0112] S103. Based on the location of the upper interface feature points and the lower interface feature points of the geological body, the boundary information of the geological body is obtained.

[0113] Among them, the boundary information of the geological body includes the extension position of the geological body interface. The extension position of the geological body interface mainly includes the interface dip, which refers to the direction of the interface tilt, and the direction of the vertical strike line when projected downward along the tilt surface onto the horizontal plane.

[0114] In some embodiments, the step of obtaining the boundary information of a geological body based on the locations of feature points at its upper and lower interfaces includes:

[0115] S1031. Based on the location of the upper interface feature point and the lower interface feature point of the geological body, determine the axial constraint quantity at the first arbitrary point in the underground space.

[0116] For example, when the x-coordinate of the maximum point is equal to the x-coordinate of the minimum point (or the difference between the x-coordinates of the maximum and minimum points is within a preset range), the boundary surface of the geological body is determined to be perpendicular to the observation surface. When the x-coordinates of the maximum and minimum points are not equal (or the difference between the x-coordinates of the maximum and minimum points is not within a preset range), the boundary surface of the geological body is determined to be inclined to the observation surface, that is, when... At that time, the tilt angle of the feature surface ;when At that time, the tilt angle of the feature surface The formula is as follows:

[0117]

[0118] make:

[0119]

[0120]

[0121]

[0122] in, , , as well as All of these represent intermediate parameters. Represents the x-coordinate of the first arbitrary point. Represents the ordinate of the first arbitrary point.

[0123] In this way, the axial constraint at the first arbitrary point in the underground space can be obtained. The specific formula is as follows:

[0124]

[0125] S1032. Construct the model weighting function at the first arbitrary point in the underground space using axis constraints.

[0126] In some embodiments, the step of constructing a model weighting function at a first arbitrary point in the underground space using axis constraints includes:

[0127] First, normalize the axial constraints at all first arbitrary points in the underground space to obtain the normalized weighted matrix.

[0128] For example, the formula for normalizing the axial constraint at all first arbitrary points in the underground space is shown below:

[0129]

[0130] in, It is used to avoid tiny positive numbers with a denominator of 0, and can be taken as 0.001.

[0131] Perform diagonal matrix operations on the normalized weighted matrix to obtain the first matrix.

[0132] For example, the specific formula for performing diagonal matrix operations on the normalized weighted matrix is ​​as follows:

[0133]

[0134] in, This represents diagonal matrix operations.

[0135] III. Based on the first matrix, obtain the model weighting function at the first arbitrary point in the underground space.

[0136] For example, determining the model weighting function at the first arbitrary point in the underground space. The formula is expressed as follows:

[0137]

[0138] in, Represents the y-coordinate of the first arbitrary point. This represents the ordinate of the observation point when the observation surface is horizontal. It is a constant.

[0139] Therefore, the model weighting function composed of all the partition points of the underground space can be obtained. The expression is as follows:

[0140]

[0141] S1033, Obtain the data weighting function.

[0142] In some embodiments, the data weighting function is obtained through the following steps:

[0143] ① Based on the coordinates of the first arbitrary point in the underground space and the coordinates of the ground observation point, construct the first forward operator corresponding to the first arbitrary point of the second arbitrary point in the ground observation point.

[0144] For example, constructing a second arbitrary point among ground observation points. Corresponding to the first arbitrary point First orthogonal operator To obtain the second arbitrary point The second forward modeling operator corresponding to the underground space grid structure ,in, The coordinates are , The coordinates are The first forward operator The expression is:

[0145]

[0146] in, Represents the gravitational constant, and its values ​​can be... ; express The size of the directional mesh. express The size of the directional mesh.

[0147] ② Based on the first forward operator, obtain the second forward operator corresponding to the second arbitrary point.

[0148] For example, the second forward operator The expression is:

[0149]

[0150] ③ Obtain the data weighting function through the second forward operator.

[0151] For example, according to the second forward operator Generate data weighting function The specific formula is as follows:

[0152]

[0153] in, This indicates the transpose operation.

[0154] S1034. Based on the model weighting function and the data weighting function, the physical property parameter model is obtained.

[0155] For example, when obtaining the physical property parameter model, firstly, according to the model weighting function... Determine the weighted initial model parameters The specific formula is as follows:

[0156]

[0157] in, This represents the density value of each grid cell in the initial underground space, which can be set to 0 by default.

[0158] Then, according to the data weighting function Second forward operator and model weighting function Determine the weighted forward operators The specific formula is as follows:

[0159]

[0160] Secondly, it can be achieved through weighted forward operators. The specific formulas for determining physical property parameters are as follows:

[0161]

[0162]

[0163]

[0164] in, , as well as These are physical property parameters; This represents the regularization parameter, which can generally be selected through calculation or experience. Here, we can take the value 2 for calculation. Represents the identity matrix; This represents gravity data.

[0165] Then, based on the physical property parameters , and initial model parameters Determine the objective function For initial model parameters derivative The expression is as follows:

[0166]

[0167] Then, the corresponding initial search step size can be determined based on the derivative. The formula is as follows:

[0168]

[0169] in, Represents process quantities during the operation. This indicates the transpose operation.

[0170] Then, based on the initial model parameters Initial search step size and process quantities The conjugate gradient algorithm is used to iteratively calculate and update the physical property parameter model of the objective function. Specifically, the steps are as follows: when the number of iterations... or data fitting function At that time, the weighted physical property parameter model is updated in the weighted parameter domain, as shown in the following formula:

[0171]

[0172] in, , Indicates the first The weighted physical property parameter model after the next iteration Indicates the first Next search step size, Indicates the first Secondary weighted property parameter model Search direction.

[0173] Then, the first Weighted property parameter model after the second iteration The conversion is performed, and the expression is as follows:

[0174]

[0175] in, Represents the physical property parameter model, This represents the prior density parameter.

[0176] Then, based on the physical property parameter model Determine the data fitting function The formula is as follows:

[0177]

[0178] Then, based on the physical property parameters , and the Weighted property parameter model after the second iteration Determine the objective function For the first Weighted property parameter model after the second iteration first derivative The expression is as follows:

[0179]

[0180] In the loop, the first Secondary weighted property parameter model Search direction The expression is:

[0181]

[0182] in, Describe the objective function For the first Weighted property parameter model after the second iteration The first derivative, Describe the objective function For the first Weighted property parameter model after the second iteration The first derivative; This represents intermediate parameters.

[0183] Secondly, the step size corresponding to the search direction The expression is:

[0184]

[0185] Finally, the loop ends and the physical property parameter model after gain is output. ,in, Indicates the first Secondary property parameter model Search direction.

[0186] S1035. The linear information in the extractive parameter model is extracted using the finite difference method to obtain the boundary information of the geological body.

[0187] For example, the linear features of the gravity data boundary are extracted and obtained using the finite difference method. The expression is as follows:

[0188]

[0189]

[0190]

[0191] in, This indicates that linear features in the X direction enhance linear information. This indicates that the linear features in the Z direction enhance the linear information. underground space Orientation mesh size, Indicates underground space Orientation mesh size.

[0192] For example, through linear features It can obtain the boundary information of geological bodies, including the extension location of the geological body boundary and maps, including the model forward curve, model location and boundary identification results.

[0193] For example, Table 1 shows the theoretical model parameters for the experiment. It can be seen that the burial depth of the upper and lower bases, the tilt angle, and the corner points of the upper interface of model bodies A, B, C, and D are... Axial projection position and residual density.

[0194] Table 1 - Theoretical Model Parameter Table

[0195]

[0196] Please refer to Figures 3 to 6 , Figures 3 to 6 For the graphics corresponding to Table 1, from Figure 3 As can be seen, the upper base of model A is buried at a depth of 5m, the lower base at a depth of 20m, and the inclination angle is 45°. Figure 4 It can be seen from this that the upper base of model B has a burial depth of 5m, the lower base has a burial depth of 20m, and the inclination angle is 60°. Figure 5 It can be seen from this that the upper base of model C has a burial depth of 5m, the lower base has a burial depth of 20m, and the tilt angle is 90°. Figure 6 It can be seen that the upper base of model D is buried at a depth of 5m, the lower base at a depth of 20m, and the tilt angle is 120°.

[0197] This application provides a boundary identification method based on gravity data. The boundary identification method includes: acquiring first information; wherein the first information includes gravity data of a ground observation point and the coordinates of the ground observation point; using the first information to obtain the positions of the upper interface feature points and the lower interface feature points of a geological body in underground space; and obtaining the boundary information of the geological body based on the positions of the upper interface feature points and the lower interface feature points of the geological body; wherein the boundary information of the geological body includes the extension position of the interface of the geological body.

[0198] Thus, the technical solution provided in this application can identify the boundary of a geological body by the location of the upper interface feature points and the lower interface feature points in the underground space, thereby achieving the identification of inclined boundaries. It is not limited by the premise assumption of the vertical angle of the boundary, effectively improving the reliability of the geological body boundary identification results, and thus providing more solid evidence for subsequent work such as geological interpretation.

[0199] Based on the same application concept, this application also provides a boundary recognition device based on gravity data, which corresponds to the boundary recognition method based on gravity data provided in the above embodiment. Since the principle of the device in this application is similar to the boundary recognition method based on gravity data in the above embodiment, the implementation of the device can refer to the implementation of the method, and the repeated parts will not be described again.

[0200] This application provides an embodiment of a boundary recognition device based on gravity data. Please refer to [link to relevant documentation]. Figure 7 , Figure 7 This is a schematic diagram of a boundary recognition device based on gravity data provided in an embodiment of this application, as shown below. Figure 7 As shown, the boundary identification device 710 includes:

[0201] The acquisition module 711 is used to acquire first information; wherein, the first information includes gravity data of the ground observation point and the coordinates of the ground observation point;

[0202] Processing module 712 is used to obtain the location of the upper interface feature point and the lower interface feature point of the geological body in the underground space using the first information;

[0203] The identification module 713 is used to obtain the boundary information of the geological body based on the location of the upper interface feature points and the location of the lower interface feature points of the geological body; wherein, the boundary information of the geological body includes the extension position of the interface of the geological body.

[0204] Optionally, the processing module 712 is specifically used for:

[0205] Using the first information, gravity anomaly data corresponding to any depth in the underground space are obtained;

[0206] Obtain the scaling function corresponding to any depth in the underground space;

[0207] Based on the gravity anomaly data and the scaling function, the density distribution intensity at the first arbitrary point in the underground space is obtained;

[0208] In a first set consisting of the density distribution intensity at all first arbitrary points in the underground space, the positions of the upper interface feature points and the lower interface feature points of the geological body in the underground space are obtained through the extreme points in the first set.

[0209] Optionally, the identification module 713 is specifically used for:

[0210] Based on the location of the upper interface feature point and the location of the lower interface feature point of the geological body, the axial constraint quantity at the first arbitrary point in the underground space is determined.

[0211] Using the axial constraints, a model weighting function is constructed at the first arbitrary point in the underground space;

[0212] Get the data weighting function;

[0213] Based on the model weighting function and the data weighting function, a physical property parameter model is obtained;

[0214] The linear information in the physical property parameter model is extracted using the finite difference method to obtain the boundary information of the geological body.

[0215] Optionally, when the processing module 712 is used to obtain gravity anomaly data corresponding to any depth in the underground space using the first information, the processing module 712 is specifically used for:

[0216] The gravity anomaly spectrum of the ground observation point is obtained by using the Fourier transform method on the gravity data of the ground observation point in the first information.

[0217] Based on the gravity anomaly spectrum of the ground observation points, determine the gravity anomaly spectrum at any depth in the underground space;

[0218] The gravity anomaly spectrum at any depth in the underground space is subjected to Fourier transform to obtain gravity anomaly data corresponding to any depth in the underground space.

[0219] Optionally, when the identification module 713 is used to construct the model weighting function at the first arbitrary point in the underground space using the axis constraint, the identification module 713 is specifically used for:

[0220] The axial constraint quantities at all first arbitrary points in the underground space are normalized to obtain a normalized weighted matrix;

[0221] Perform diagonal matrix operations on the normalized weighted matrix to obtain the first matrix;

[0222] Based on the first matrix, the model weighting function at the first arbitrary point in the underground space is obtained.

[0223] Optionally, when the identification module 713 is used to obtain the data weighting function, the identification module 713 is further used to:

[0224] Based on the coordinates of the first arbitrary point in the underground space and the coordinates of the ground observation point, a first forward operator is constructed corresponding to the first arbitrary point of the second arbitrary point in the ground observation point.

[0225] Based on the first forward operator, the second forward operator corresponding to the second arbitrary point is obtained;

[0226] The data weighting function is obtained through the second forward operator.

[0227] This application provides a boundary recognition device based on gravity data. The boundary recognition device includes: an acquisition module for acquiring first information, wherein the first information includes gravity data of a ground observation point and the coordinates of the ground observation point; a processing module for using the first information to obtain the positions of the upper interface feature points and the lower interface feature points of a geological body in underground space; and an identification module for obtaining the boundary information of the geological body based on the positions of the upper interface feature points and the lower interface feature points of the geological body, wherein the boundary information of the geological body includes the extension position of the interface of the geological body.

[0228] Thus, the technical solution provided in this application can identify the boundary of a geological body by the location of the upper interface feature points and the lower interface feature points in the underground space, thereby achieving the identification of inclined boundaries. It is not limited by the premise assumption of the vertical angle of the boundary, effectively improving the reliability of the geological body boundary identification results, and thus providing more solid evidence for subsequent work such as geological interpretation.

[0229] Please see Figure 8 , Figure 8 This is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. Figure 8 As shown, the electronic device 800 includes a processor 810, a memory 820, and a bus 830.

[0230] The memory 820 stores machine-readable instructions executable by the processor 810. When the electronic device 800 is running, the processor 810 and the memory 820 communicate via the bus 830. When the machine-readable instructions are executed by the processor 810, they can perform the operations described above. Figure 1 The steps of the boundary recognition method based on gravity data in the method embodiment shown are described in detail in the method embodiment, and will not be repeated here.

[0231] This application also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, can perform the above-described actions. Figure 1 The steps of the boundary recognition method based on gravity data in the method embodiment shown are described in detail in the method embodiment, and will not be repeated here.

[0232] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0233] In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods can be implemented in other ways. The apparatus embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. Furthermore, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the shown or discussed mutual couplings, direct couplings, or communication connections may be through some communication interfaces; indirect couplings or communication connections between devices or units may be electrical, mechanical, or other forms.

[0234] In addition, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0235] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a processor-executable, non-volatile, computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0236] Finally, it should be noted that the above-described embodiments are merely specific implementations of this application, used to illustrate the technical solutions of this application, and not to limit them. The scope of protection of this application is not limited thereto. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments, or make equivalent substitutions for some of the technical features, within the scope of the technology disclosed in this application. Such modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application, and should all be covered within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A boundary recognition method based on gravity data, characterized in that, The boundary recognition method includes: Obtain first information; wherein, the first information includes gravity data of the ground observation point and the coordinates of the ground observation point; The location of the upper interface feature point and the lower interface feature point of the geological body in the underground space are obtained using the first information. Based on the locations of the upper and lower interface feature points of the geological body, the boundary information of the geological body is obtained; wherein, the boundary information of the geological body includes the extension location of the interface of the geological body. The step of obtaining the location of the upper interface feature point and the lower interface feature point of the geological body in the underground space using the first information includes: Using the first information, gravity anomaly data corresponding to any depth in the underground space are obtained; Obtain the scaling function corresponding to any depth in the underground space; Based on the gravity anomaly data and the scaling function, the density distribution intensity at the first arbitrary point in the underground space is obtained; In a first set consisting of density distribution intensities at all first arbitrary points in the underground space, the positions of the upper interface feature points and the lower interface feature points of the geological body in the underground space are obtained through the extreme points in the first set. The step of obtaining the boundary information of the geological body based on the positions of the upper and lower interface feature points of the geological body includes: Based on the location of the upper interface feature point and the location of the lower interface feature point of the geological body, the axial constraint quantity at the first arbitrary point in the underground space is determined. Using the axial constraints, a model weighting function is constructed at the first arbitrary point in the underground space; Obtain the data weighting function; Based on the model weighting function and the data weighting function, a physical property parameter model is obtained; The linear information in the physical property parameter model is extracted using the finite difference method to obtain the boundary information of the geological body.

2. The boundary recognition method according to claim 1, characterized in that, The step of obtaining gravity anomaly data corresponding to any depth in the underground space using the first information includes: The gravity anomaly spectrum of the ground observation point is obtained by using the Fourier transform method on the gravity data of the ground observation point in the first information. Based on the gravity anomaly spectrum of the ground observation points, determine the gravity anomaly spectrum at any depth in the underground space; The gravity anomaly spectrum at any depth in the underground space is subjected to Fourier transform to obtain gravity anomaly data corresponding to any depth in the underground space.

3. The boundary recognition method according to claim 1, characterized in that, The step of constructing the model weighting function at the first arbitrary point in the underground space using the axis constraints includes: The axial constraint quantities at all first arbitrary points in the underground space are normalized to obtain a normalized weighted matrix; Perform diagonal matrix operations on the normalized weighted matrix to obtain the first matrix; Based on the first matrix, the model weighting function at the first arbitrary point in the underground space is obtained.

4. The boundary recognition method according to claim 1, characterized in that, Obtain the data weighting function through the following steps: Based on the coordinates of the first arbitrary point in the underground space and the coordinates of the ground observation point, a first forward operator is constructed corresponding to the first arbitrary point of the second arbitrary point in the ground observation point. Based on the first forward operator, the second forward operator corresponding to the second arbitrary point is obtained; The data weighting function is obtained through the second forward operator.

5. A boundary recognition device based on gravity data, characterized in that, The boundary identification device includes: An acquisition module is used to acquire first information; wherein, the first information includes gravity data of a ground observation point and the coordinates of the ground observation point; The processing module is used to obtain the location of the upper interface feature point and the lower interface feature point of the geological body in the underground space using the first information; The identification module is used to obtain the boundary information of the geological body based on the location of the upper interface feature points and the location of the lower interface feature points; wherein, the boundary information of the geological body includes the extension position of the interface of the geological body; The processing module is specifically used for: Using the first information, gravity anomaly data corresponding to any depth in the underground space are obtained; Obtain the scaling function corresponding to any depth in the underground space; Based on the gravity anomaly data and the scaling function, the density distribution intensity at the first arbitrary point in the underground space is obtained; In a first set consisting of density distribution intensities at all first arbitrary points in the underground space, the positions of the upper interface feature points and the lower interface feature points of the geological body in the underground space are obtained through the extreme points in the first set. The identification module is specifically used for: Based on the location of the upper interface feature point and the location of the lower interface feature point of the geological body, the axial constraint quantity at the first arbitrary point in the underground space is determined. Using the axial constraints, a model weighting function is constructed at the first arbitrary point in the underground space; Obtain the data weighting function; Based on the model weighting function and the data weighting function, a physical property parameter model is obtained; The linear information in the physical property parameter model is extracted using the finite difference method to obtain the boundary information of the geological body.

6. An electronic device, characterized in that, include: The device includes a processor, a memory, and a bus. The memory stores machine-readable instructions executable by the processor. When the electronic device is running, the processor communicates with the memory via the bus. The machine-readable instructions are executed by the processor to perform the steps of the boundary recognition method based on gravity data as described in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, performs the steps of the boundary recognition method based on gravity data as described in any one of claims 1 to 4.

Images