Method and device for determining the position of a fracture based on gravity data

Abstract

The application provides a method and device for determining a fracture position based on gravity data, and the method comprises the following steps: obtaining basic data of the gravity data, i.e., Bouguer gravity anomaly data; obtaining residual gravity anomaly gridded data from the Bouguer gravity anomaly gridded data, and performing pseudo-magnetic anomaly conversion on the residual gravity anomaly gridded data to obtain pseudo-magnetic anomaly gridded data; performing horizontal derivative calculation on the pseudo-magnetic anomaly gridded data to obtain pseudo-magnetic anomaly horizontal derivative gridded data; drawing a pseudo-magnetic anomaly horizontal derivative graph according to the pseudo-magnetic anomaly horizontal derivative gridded data; identifying maximum value positions on the pseudo-magnetic anomaly horizontal derivative graph, connecting the maximum value positions to obtain a fracture position. The application can reduce the horizontal gradient maximum value strip width of the fracture studied by using the gravity data by half, thereby improving the determination precision of the fracture position and improving the identification ability of the fracture with small scale and large depth.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of gravity exploration data processing, and particularly relates to a method and device for determining the position of a fracture based on gravity data. BACKGROUND

[0002] A fracture is a basic geological phenomenon in nature, and human activities and energy and mineral exploration and development are closely related to the fracture. Technical personnel have been exploring methods for improving the accuracy of fracture position interpretation. Due to the different densities of rocks of different strata and lithology, there is a gravity difference on both sides of the fracture, which makes it possible to study the fracture by using gravity data. Gravity method researchers have been exploring new methods to continuously improve the accuracy of determining the position of the fracture by using gravity data. Bouguer gravity anomaly is the basic data for gravity exploration. The fracture appears as a gravity gradient change zone (sometimes referred to as a gravity step zone) on the Bouguer gravity anomaly and the residual gravity anomaly obtained by removing the regional gravity field from the Bouguer gravity anomaly. People can understand the approximate position of the fracture by the position of the gravity gradient change zone of the Bouguer gravity anomaly or the residual gravity anomaly. However, due to the existence of a certain width of the gravity gradient change zone, it is impossible to directly determine the accurate position of the fracture from the zone. In order to obtain a more accurate fracture position, gravity researchers have been exploring processing and interpretation methods. In summary, there are two steps. The first step is to process the Bouguer gravity anomaly or the residual gravity anomaly to obtain a gravity anomaly that can highlight the position of the fracture. The traditional and representative method is to calculate the horizontal derivative to obtain a gravity horizontal derivative map (also known as a gravity horizontal gradient map). The second step is to draw a gravity anomaly map reflecting the position of the fracture. The position of the fracture is identified on the map. For the gravity horizontal derivative map, the position of the maximum value is identified on the map, and the line connecting the positions of the maximum values is interpreted as the position of the fracture. The above method is correct and effective, but there is a problem of unsatisfactory accuracy of the position of the fracture.

[0003] Because the gravity gradient change band is wide, the maximum of the horizontal derivative is obviously narrowed in a relatively narrow range, but the accuracy of determining the accurate position of the fault from the narrow band is still not ideal. In order to better solve the problem, researchers have taken various technical measures. The measures taken in the first step above are of several types: the first type is to filter the Bouguer gravity anomaly or residual gravity anomaly in a small sub-domain, wavelet transform and the like, so as to narrow the gravity gradient change band, and then calculate the horizontal derivative; the second type is to calculate the gravity horizontal total gradient, gravity total gradient and the like of the Bouguer gravity anomaly or residual gravity anomaly, instead of the gravity horizontal derivative, so as to narrow the range of the maximum value band; the third type is to further narrow the range of the maximum value band by calculating the vertical derivative, high-pass filtering and the like after calculating the horizontal derivative. The three types of methods all belong to the method of calculating the horizontal derivative of the Bouguer gravity anomaly or residual gravity anomaly, and mainly perform information enhancement before or after the derivative is calculated. Some researchers also combine two or more of the above methods to enhance the processing effect. The measures taken in the second step are of several types: the first type is to use the image method, that is, to simulate the horizontal gradient data as a three-dimensional relief-like terrain, and then to irradiate with a certain height angle, and the limit of light and shade change is the position of the maximum value, which is explained as the position of the fault; the second type is the value-seeking method, which finds the position of the maximum value through mathematical operation, and takes the line connecting the maximum values as the position of the fault. The two types of methods are both based on the gravity anomaly data obtained in the first step, and can improve the fault interpretation accuracy, but the accuracy is limited by the accuracy of the gravity anomaly data obtained in the first step.

[0004] How to narrow the width of the horizontal gradient maximum value band is a technical problem to be solved urgently to improve the accuracy of determining the fault position based on gravity data. SUMMARY

[0005] In view of the problems in the prior art, the present application provides a method and device for determining the position of a fault based on gravity data, which improves the accuracy of determining the position of a fault.

[0006] To solve the above technical problems, the present application provides the following technical solutions:

[0007] In a first aspect, the present application provides a method for determining the position of a fault based on gravity data, comprising:

[0008] obtaining Bouguer gravity anomaly data as the basic data of gravity data;

[0009] obtaining residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data, performing pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data to obtain pseudo-magnetic anomaly gridding data;

[0010] The pseudo-magnetic force anomaly gridded data is subjected to horizontal derivative calculation to obtain pseudo-magnetic force anomaly horizontal derivative gridded data.

[0011] A pseudo-magnetic force anomaly horizontal derivative map is drawn according to the pseudo-magnetic force anomaly horizontal derivative gridded data.

[0012] A maximum value position is identified on the pseudo-magnetic force anomaly horizontal derivative map, and each maximum value position is connected to obtain a fracture position.

[0013] After the Bouguer gravity anomaly data of the basic data of the gravity data is obtained, the method further comprises:

[0014] The Bouguer gravity anomaly data is subjected to gridded interpolation processing to obtain Bouguer gravity anomaly gridded data.

[0015] If a regional background gravity field exists in the Bouguer gravity anomaly gridded data, then residual gravity anomaly gridded data is obtained from the Bouguer gravity anomaly gridded data.

[0016] The residual gravity anomaly gridded data is obtained from the Bouguer gravity anomaly gridded data, and comprises:

[0017] The residual gravity anomaly gridded data is obtained from the Bouguer gravity anomaly gridded data by using a potential field continuation method or a wave number domain filtering method.

[0018] The pseudo-magnetic force anomaly conversion of the residual gravity anomaly gridded data comprises:

[0019] The pseudo-magnetic force anomaly conversion of the residual gravity anomaly gridded data is performed by using a Poisson formula of magnetic potential and gravitational potential of a same uniform geological body.

[0020] In a second aspect, the present application provides a device for determining a fracture position based on gravity data, comprising:

[0021] An acquisition unit is configured to acquire Bouguer gravity anomaly data of basic data of gravity data.

[0022] An extraction unit is configured to obtain residual gravity anomaly gridded data from the Bouguer gravity anomaly gridded data, and to perform pseudo-magnetic force anomaly conversion on the residual gravity anomaly gridded data to obtain pseudo-magnetic force anomaly gridded data.

[0023] A calculation unit is configured to perform horizontal derivative calculation on the pseudo-magnetic force anomaly gridded data to obtain pseudo-magnetic force anomaly horizontal derivative gridded data.

[0024] A drawing unit is configured to draw a pseudo-magnetic force anomaly horizontal derivative map according to the pseudo-magnetic force anomaly horizontal derivative gridded data.

[0025] A determining unit is configured to identify maximum value positions on the horizontal derivative map of the pseudo magnetic force anomaly, connect each maximum value position, and obtain the fracture position.

[0026] Further comprising:

[0027] An interpolation unit is configured to perform grid interpolation processing on the Bouguer gravity anomaly data to obtain Bouguer gravity anomaly grid data.

[0028] A gravity field unit is configured to determine a regional background gravity field in the Bouguer gravity anomaly grid data, and obtain residual gravity anomaly grid data from the Bouguer gravity anomaly grid data.

[0029] The extraction unit comprises:

[0030] A first extraction subunit is configured to obtain residual gravity anomaly grid data from the Bouguer gravity anomaly grid data by using a potential field continuation method or a wave number domain filtering method.

[0031] The extraction unit comprises:

[0032] A second extraction subunit is configured to convert the residual gravity anomaly grid data into pseudo magnetic force anomaly by using a Poisson formula of the magnetic potential and the gravitational potential of the same uniform geological body.

[0033] In a third aspect, the present application provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps of the method for determining a fracture position based on gravity data.

[0034] In a fourth aspect, the present application provides a computer readable storage medium having a computer program stored thereon, wherein the computer program is executable on a processor to implement the steps of the method for determining a fracture position based on gravity data.

[0035] According to the above technical solution, the present application provides a method and device for determining a fracture position based on gravity data, which obtains Bouguer gravity anomaly data as the basic data of gravity data, obtains residual gravity anomaly grid data from the Bouguer gravity anomaly grid data, converts the residual gravity anomaly grid data into pseudo magnetic force anomaly to obtain pseudo magnetic force anomaly grid data, performs horizontal derivative calculation on the pseudo magnetic force anomaly grid data to obtain pseudo magnetic force anomaly horizontal derivative grid data, draws a pseudo magnetic force anomaly horizontal derivative map according to the pseudo magnetic force anomaly horizontal derivative grid data, identifies maximum value positions on the pseudo magnetic force anomaly horizontal derivative map, connects each maximum value position, and obtains the fracture position. Attached Figure Description

[0036] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0037] Figure 1 This is a schematic diagram of the first process of the method for determining the fracture location based on gravity data in an embodiment of the present invention.

[0038] Figure 2 This is a schematic diagram of the second process of the method for determining the fracture location based on gravity data in an embodiment of the present invention.

[0039] Figure 3 This is a schematic diagram of the step model in the method for determining the fracture location based on gravity data in an embodiment of the present invention.

[0040] Figure 4 This is a schematic diagram of the gravity anomaly curve of the step model in the method for determining the fracture location based on gravity data in an embodiment of the present invention.

[0041] Figure 5 This is a schematic diagram of the pseudo-magnetic anomaly curve of the step model in the method for determining the fracture location based on gravity data in an embodiment of the present invention.

[0042] Figure 6 This is a schematic diagram of the horizontal derivative curve and half-maximum width of the pseudo-magnetic anomaly of the step model in the method for determining the fracture location based on gravity data in this embodiment of the invention.

[0043] Figure 7 This is a schematic diagram of the horizontal derivative curve of gravity anomaly and the half-maximum width of the step model in the method for determining the fracture location based on gravity data in this embodiment of the invention.

[0044] Figure 8 This is a schematic diagram of the Bouguer gravity anomaly map in the method for determining the fracture location based on gravity data in an embodiment of the present invention.

[0045] Figure 9 This is a schematic diagram of a pseudo-magnetic anomaly map in the method for determining the fracture location based on gravity data in an embodiment of the present invention.

[0046] Figure 10 This is a schematic diagram of the total gradient of pseudo-magnetic anomalies in the method for determining the fracture location based on gravity data in an embodiment of the present invention.

[0047] Figure 11It is a schematic diagram of a Bouguer gravity anomaly horizontal total gradient map for a method for determining a fracture position based on gravity data in an embodiment of the present application.

[0048] Figure 12 It is a structural schematic diagram for determining a fracture position based on gravity data in an embodiment of the present application.

[0049] Figure 13 It is a structural schematic diagram of an electronic device in an embodiment of the present application. DETAILED DESCRIPTION

[0050] To make the objectives, technical solutions and advantages of embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are some but not all of the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by a person of ordinary skill in the art without creative work fall within the protection scope of the present application.

[0051] The present application provides an embodiment of a method for determining a fracture position based on gravity data, referring to Figure 1 The method for determining a fracture position based on gravity data specifically includes the following contents:

[0052] S101: Obtain Bouguer gravity anomaly data of basic data of gravity data;

[0053] In this step, gravity exploration is deployed in a study area, gravity data collection is performed, point coordinates and gravity values are obtained, and various correction calculations such as normal field correction, elevation correction, intermediate layer correction and terrain correction are performed to obtain Bouguer gravity anomaly data.

[0054] S102: Obtain residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data, perform pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data, and obtain pseudo-magnetic anomaly gridding data;

[0055] In this step, residual gravity anomaly gridding data is obtained from the Bouguer gravity anomaly gridding data by using a potential field continuation method or a wave number domain filtering method.

[0056] The Bouguer gravity anomaly data includes the Bouguer gravity anomaly gridding data.

[0057] The Poisson formula of the magnetic potential and the gravitational potential of the same uniform geological body is used to perform pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data.

[0058] In this step, the Poisson formula of the magnetic potential and the gravitational potential is:

[0059]

[0060] In the formula, U is a magnetic potential of the geologic body, V is a gravity potential of the geologic body, J is a magnetization intensity, G is a gravity constant, and δ is a density.

[0061] S103: horizontal derivative calculation is performed on the pseudo-magnetic force anomaly gridded data to obtain pseudo-magnetic force anomaly horizontal derivative gridded data;

[0062] In this step, the horizontal derivative can be a horizontal direction derivative, a horizontal total gradient or a total gradient.

[0063] S104: a pseudo-magnetic force anomaly horizontal derivative map is drawn according to the pseudo-magnetic force anomaly horizontal derivative gridded data;

[0064] In this step, the map can be a plane contour map or a three-dimensional image map.

[0065] S105: maximum value positions are identified on the pseudo-magnetic force anomaly horizontal derivative map, and each maximum value position is connected to obtain a fracture position.

[0066] In this step, maximum value positions are identified on the pseudo-magnetic force anomaly horizontal derivative map, and each related maximum value position is connected to obtain a fracture plane position, so that the fracture position is accurately determined based on the gravity data.

[0067] As can be seen from the above description, the present application provides a method for determining a fracture position based on gravity data, which comprises the following steps: obtaining basic data of the gravity data, i.e., Bouguer gravity anomaly data; obtaining residual gravity anomaly gridded data from the Bouguer gravity anomaly gridded data, performing pseudo-magnetic force anomaly conversion on the residual gravity anomaly gridded data to obtain pseudo-magnetic force anomaly gridded data, performing horizontal derivative calculation on the pseudo-magnetic force anomaly gridded data to obtain pseudo-magnetic force anomaly horizontal derivative gridded data, drawing a pseudo-magnetic force anomaly horizontal derivative map according to the pseudo-magnetic force anomaly horizontal derivative gridded data, and identifying maximum value positions on the pseudo-magnetic force anomaly horizontal derivative map and connecting each maximum value position to obtain a fracture position. The method can reduce the horizontal gradient maximum value strip width by half when the fracture is studied based on the gravity data, thereby improving the determination accuracy of the fracture position and improving the identification ability of the fracture with small size and large depth.

[0068] In an embodiment of the present application, referring to Figure 2 , the method for determining a fracture position based on gravity data comprises the following steps after step S101:

[0069] S106: gridded interpolation processing is performed on the Bouguer gravity anomaly data to obtain Bouguer gravity anomaly gridded data;

[0070] In this step, the grid spacing is the geodetic distance, and the gridding range is the data range. The purpose of this step is to prepare the data in the required format for subsequent processing.

[0071] S107: Determine whether there is a regional background gravity field in the Bouguer gravity anomaly gridded data, and obtain the residual gravity anomaly gridded data from the Bouguer gravity anomaly gridded data.

[0072] In this step, in combination with known historical geological data and historical exploration results, it is determined whether there is a regional background gravity field in the Bouguer gravity anomaly. If there is, residual gravity anomaly extraction processing is performed on the Bouguer gravity anomaly gridded data to obtain residual gravity anomaly gridded data; if there is no obvious regional background field, residual gravity anomaly extraction processing can not be performed.

[0073] The extraction of residual gravity anomaly can adopt various effective methods such as potential field continuation method and wave number domain filtering method, and the obtained residual gravity anomaly should suppress the regional background field to the greatest extent.

[0074] As can be seen from the above description, the technical scheme provided by the embodiment adopts the method for determining the fracture position based on gravity data, so that the fracture position determined by the gravity data is more accurate.

[0075] The embodiment of the application provides a specific example of the method for determining the fracture position based on gravity data, and specifically includes the following contents:

[0076] The following respectively shows the implementation steps and method effects of determining the fracture position by using model gravity profile data and measured gravity plane data.

[0077] First, the implementation steps and effects of determining the fracture position by using model gravity profile data are shown.

[0078] The model selects a vertical step two-dimensional model, the step top surface is buried at a depth of 500 meters, the step bottom surface is buried at a depth of 1000 meters, the step (fracture) position is at the middle 0 meter, and the density difference on both sides is 1.0 g / cm 3 , see Figure 3 . The gravity calculation profile vertically passes through the step. Through gravity forward calculation, the gravity anomaly curve of the above-mentioned model is obtained, see Figure 4 , which corresponds to step one of implementing the application, obtaining gravity anomaly data.

[0079] For gravity profile research, data gridding is not required, so step two is not implemented.

[0080] When setting the model, no regional background field is superimposed, so step three is not implemented.

[0081] Step four is implemented to obtain the residual gravity anomaly of the model gravity anomaly Figure 4) to the model pseudo-magnetic anomaly curve ( Figure 5 ) In the conversion, the magnetic field strength is set to 50000 nT, the magnetic inclination is set to 90°, which is equivalent to vertical magnetization, and the magnetic susceptibility is set to 1000 x 10 -5 SI. Comparison Figure 5 and Figure 4 It can be seen that in the middle of the profile, near the position of the step (fracture) of 0 meters, the width of the gradient change of the pseudo-magnetic anomaly ( Figure 7 ) is much narrower than that of the gravity anomaly gradient change zone ( Figure 4 ).

[0082] Step five is implemented, the horizontal derivative of the pseudo-magnetic anomaly is calculated, and the pseudo-magnetic anomaly horizontal derivative data is obtained.

[0083] Step six is implemented, the pseudo-magnetic anomaly horizontal derivative curve ( Figure 6 ) is plotted with the pseudo-magnetic anomaly horizontal derivative data. According to the index of the half-extreme anomaly width often used in gravity and magnetic force research, the half-extreme anomaly width of the model pseudo-magnetic anomaly horizontal derivative is about 700 meters.

[0084] Step seven is implemented, the position of the maximum value of the pseudo-magnetic anomaly horizontal derivative is interpreted, and the position of the maximum value is the position of the fracture. From the curve, the maximum value is located at 0 meters, indicating that the position of the fracture is near 0 meters, which is consistent with the position of the step (fracture) in the step model. Thus, the fracture position is identified from the gravity data of the step model.

[0085] In order to compare the effect with the existing method of determining the fracture position from the gravity anomaly horizontal derivative, the horizontal derivative of the model gravity anomaly ( Figure 4 ) is directly calculated, and the gravity anomaly horizontal derivative curve is shown in Figure 7 The maximum value of the horizontal derivative of the gravity anomaly is located at 0 meters, which is consistent with the position reflected by the horizontal derivative of the pseudo-magnetic anomaly, indicating that the method of the present application is correct; the half-extreme anomaly width of the gravity anomaly horizontal derivative is about 1400 meters, which is twice the half-extreme anomaly width of the pseudo-magnetic anomaly horizontal derivative, indicating that the method of the present application is more accurate. It can be seen that the reflection of the fracture position by the pseudo-magnetic anomaly horizontal derivative is obviously better than that by the gravity anomaly horizontal derivative, that is, the method of the present application is superior to the existing method of directly using the gravity anomaly horizontal derivative to determine the fracture position.

[0086] Then the implementation steps and method effects of determining the fracture position from the measured gravity plane data are shown. The study area is selected at the basin boundary of a certain basin, and the fracture is obliquely through the study area.

[0087] In the first step, a 500m × 500m gravity measurement network was deployed in the study area. Through gravity data acquisition, the coordinates and gravity values ​​of each measurement point were obtained. Various correction calculations were performed, including normal field correction, elevation correction, intermediate layer correction, and terrain correction, to obtain Bouguer gravity anomaly data.

[0088] Step two involves performing gridded interpolation on the Bouguer gravity anomaly data to obtain gridded Bouguer gravity anomaly data. The grid spacing is set to the gravity point distance, i.e., 500 meters. The gridded area is the data range; blank areas in the figure represent regions where no gravity measurement points were deployed. The existing Kriging interpolation method is used for gridding. A schematic diagram of the Bouguer gravity anomaly map can be found here. Figure 8 .

[0089] Step three involves collecting and analyzing regional geological data. It is known that this area is located in a stable platform zone, and theoretically, there is no obvious regional background field. Through the collection and analysis of large-scale gravity data, it is known that there is no obvious regional background field in this area. Therefore, this embodiment does not need to extract residual gravity anomalies. The Bouguer gravity anomaly in this area is equivalent to the residual gravity anomaly. Therefore, subsequent processing is based on the Bouguer gravity anomaly data.

[0090] Step four involves using Poisson's formula for the magnetic and gravitational potentials of the same homogeneous geological body. The Bouguer gravity anomaly gridded data was transformed into pseudo-magnetic anomaly gridded data. The transformation parameters were selected according to the actual parameters of the survey area. The pseudo-magnetic anomaly map is shown below. Figure 9 .

[0091] Step five involves calculating the horizontal derivative of the pseudo-magnetic anomaly gridded data to obtain the pseudo-magnetic anomaly horizontal derivative gridded data. The horizontal gradient calculation method used here is to obtain the total horizontal gradient, the same method used later for gravity anomaly processing to facilitate comparison.

[0092] Step six involves using plotting software to create a plane contour map of the pseudo-magnetic anomaly's horizontal derivative using gridded data. (See attached image.) Figure 10 Drawing plane contour maps using gridded data is a common method, and it can be done with various existing drawing software. The graphic form here is the most common plane contour map.

[0093] Step seven involves implementing the overall gradient map at the pseudo-magnetic anomaly level. Figure 10 The location of the maximum value of the horizontal derivative of the pseudo-magnetic anomaly is identified on the graph. Connecting these maximum values ​​yields the planar location of the fracture, thus enabling the determination of the fracture location based on gravity data. This method of determining the fracture location using the maximum value of the horizontal gradient of the pseudo-magnetic anomaly is the same as the method of determining the fracture location using the maximum value of the horizontal gradient of the gravity anomaly; therefore, the fracture location is not plotted in the figure.

[0094] To compare its effectiveness with existing methods for determining fracture locations using the horizontal derivative of gravity anomalies, conventional methods were used to analyze measured gravity anomalies (…). Figure 8 The total horizontal gradient was directly calculated, and a total horizontal gradient map of gravity anomalies was plotted. Figure 11 The method of interpreting fracture locations using total horizontal gradient data of gravity anomalies involves determining the locations of maxima on the total horizontal gradient map and interpreting them as fracture locations. Therefore, when the maximum value band is narrow, the determined maximum value location is more accurate, and the interpreted fracture location has higher precision. Conversely, when the maximum value band is wide, the final interpreted fracture location will have lower precision. (Comparison with pseudo-magnetic anomaly total horizontal gradient map...) Figure 10 ) and total gradient map of gravity anomaly level ( Figure 11 As can be seen, the maximum band of the horizontal total gradient of the pseudo-magnetic anomaly is significantly narrower than that of the horizontal total gradient of the gravity anomaly. Therefore, the fracture location explained using the horizontal total gradient of the pseudo-magnetic anomaly is necessarily more accurate than that explained using the horizontal total gradient of the gravity anomaly. Furthermore, it can be seen that in the horizontal total gradient map of the pseudo-magnetic anomaly (… Figure 10 On the map, there is a noticeable break and relative displacement in the middle of the maximum value band (reflecting the fracture location), while in the total gradient map of the gravity anomaly level ( Figure 11 This phenomenon is difficult to observe on the surface, which further reflects that the present invention has a significantly better ability to distinguish fractures than existing methods that directly use the total gradient of gravity anomalies to determine the fracture location.

[0095] As can be seen from the above description, the method for determining fracture location based on gravity data provided in the embodiments of the present invention uses this technology to reduce the width of the horizontal gradient maximum stripe of fractures studied using gravity data by half, thereby improving the accuracy of fracture location determination and improving the ability to identify fractures that are small in scale but large in depth.

[0096] This invention provides a specific implementation of a device for determining fracture location based on gravity data, capable of realizing all the functions of the method for determining fracture location based on gravity data. See [link to specific implementation details]. Figure 12 The device for determining the fracture location based on gravity data specifically includes the following components:

[0097] Acquisition unit 10 is used to acquire the basic data of gravity data, Bouguer gravity anomaly data;

[0098] Extraction unit 20 is used to obtain residual gravity anomaly gridded data from the Bouguer gravity anomaly gridded data, and to perform pseudo-magnetic anomaly transformation on the residual gravity anomaly gridded data to obtain pseudo-magnetic anomaly gridded data.

[0099] The computing unit 30 is used for performing horizontal derivative calculation on the pseudo-magnetic anomaly gridded data to obtain pseudo-magnetic anomaly horizontal derivative gridded data.

[0100] The drawing unit 40 is used for drawing a pseudo-magnetic anomaly horizontal derivative graph according to the pseudo-magnetic anomaly horizontal derivative gridded data.

[0101] The determining unit 50 is used for identifying maximum value positions on the pseudo-magnetic anomaly horizontal derivative graph, connecting the maximum value positions to obtain a fracture position.

[0102] Further comprising:

[0103] The interpolation unit is used for performing gridded interpolation processing on the Bouguer gravity anomaly data to obtain Bouguer gravity anomaly gridded data.

[0104] The gravity field unit is used for determining that a regional background gravity field exists in the Bouguer gravity anomaly gridded data, and then obtaining residual gravity anomaly gridded data from the Bouguer gravity anomaly gridded data.

[0105] The extraction unit comprises:

[0106] The first extraction subunit is used for obtaining the residual gravity anomaly gridded data from the Bouguer gravity anomaly gridded data by using a potential field continuation method or a wave number domain filtering method.

[0107] The extraction unit comprises:

[0108] The second extraction subunit is used for performing pseudo-magnetic anomaly conversion on the residual gravity anomaly gridded data by using a Poisson formula of magnetic potential and gravitational potential of the same uniform geological body.

[0109] The embodiment of the device for determining a fracture position based on gravity data provided by the application can be used to execute the processing flow of the embodiment of the method for determining a fracture position based on gravity data in the above embodiment, and the function thereof will not be repeated here, and the detailed description of the above method embodiment can be referred to.

[0110] As known from the above description, the device for determining a fracture position based on gravity data provided by the embodiment of the application reduces the horizontal gradient maximum value strip width of the fracture by half by using gravity data to study the fracture, thereby improving the determination accuracy of the fracture position and improving the identification ability of the fracture with small scale and large depth.

[0111] The embodiment of the electronic device for implementing all or part of the contents of the method for determining a fracture position based on gravity data provided by the application specifically comprises the following contents:

[0112] A processor, a memory, a communications interface, and a bus; wherein the processor, the memory, the communications interface complete the communication among each other through the bus; the communications interface is used to realize the information transmission between the related devices; the electronic device can be a desktop computer, a tablet computer, a mobile terminal, and the like, and the embodiments are not limited thereto. In the embodiments, the electronic device can be implemented by referring to the embodiments of the method for determining the fracture position based on the gravity data and the embodiments of the device for determining the fracture position based on the gravity data, and the contents are incorporated herein, and the repeated parts will not be described herein.

[0113] Figure 13 A schematic block diagram of the system configuration of the electronic device 9600 of the embodiments of the present application is shown in FIG. 9. As shown in FIG. 9, the electronic device 9600 can include a central processor 9100 and a memory 9140; the memory 9140 is coupled to the central processor 9100. It is worth noting that the structure shown in FIG. 9 is exemplary; other types of structures can also be used to supplement or replace the structure to realize the telecommunication function or other functions. Figure 13 Figure 13 The structure shown in FIG. 9 is exemplary; other types of structures can also be used to supplement or replace the structure to realize the telecommunication function or other functions.

[0114] In an embodiment, the function of determining the fracture position based on the gravity data can be integrated into the central processor 9100. Wherein, the central processor 9100 can be configured to control as follows:

[0115] Obtaining the base data of the gravity data, the Bouguer gravity anomaly data; obtaining the residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data, performing the pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data to obtain the pseudo-magnetic anomaly gridding data; performing the horizontal derivative calculation on the pseudo-magnetic anomaly gridding data to obtain the pseudo-magnetic anomaly horizontal derivative gridding data; drawing the pseudo-magnetic anomaly horizontal derivative graph according to the pseudo-magnetic anomaly horizontal derivative gridding data; identifying the maximum value positions on the pseudo-magnetic anomaly horizontal derivative graph, connecting each maximum value position to obtain the fracture position.

[0116] ​As can be seen from the above description, the electronic device provided by the embodiment of the application obtains the Bouguer gravity anomaly data as the basic data of the gravity data; obtains the residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data, performs pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data to obtain the pseudo-magnetic anomaly gridding data; performs horizontal derivative calculation on the pseudo-magnetic anomaly gridding data to obtain the pseudo-magnetic anomaly horizontal derivative gridding data; draws the pseudo-magnetic anomaly horizontal derivative graph according to the pseudo-magnetic anomaly horizontal derivative gridding data; and identifies the maximum value positions on the pseudo-magnetic anomaly horizontal derivative graph, connects the various maximum value positions to obtain the fracture position, so that the width of the horizontal gradient maximum value strip of the fracture studied by using the gravity data is reduced by half, thereby improving the determination accuracy of the fracture position and improving the identification ability of the fracture with small scale and large depth.

[0117] In another embodiment, the device for determining the fracture position based on the gravity data can be configured separately from the central processor 9100, for example, the device for determining the fracture position based on the gravity data can be configured as a chip connected with the central processor 9100, and the function of determining the fracture position based on the gravity data is realized through the control of the central processor.

[0118] As shown in Figure 13 , the electronic device 9600 can further include a communication module 9110, an input unit 9120, an audio processor 9130, a display 9160, and a power supply 9170. It should be noted that the electronic device 9600 does not necessarily include all the components shown in Figure 13 ; in addition, the electronic device 9600 can further include components not shown in Figure 13 , which can be referred to the prior art.

[0119] As shown in Figure 13 , the central processor 9100, sometimes also referred to as a controller or operation control, can include a microprocessor or other processor device and / or logic device, which receives input and controls the operation of various components of the electronic device 9600.

[0120] The memory 9140, for example, can be one or more of a buffer, a flash memory, a hard drive, a removable media, a volatile memory, a non-volatile memory, or other suitable device. The information related to the failure can be stored, and in addition, programs for executing the information can be stored. The central processor 9100 can execute the programs stored in the memory 9140 to realize information storage or processing, etc.

[0121] The input unit 9120 provides input to the central processing unit 9100. The input unit 9120 is, for example, a key or a touch input device. The power supply 9170 is for supplying power to the electronic device 9600. The display 9160 is for displaying display objects such as images and characters. The display is, for example, an LCD display, but is not limited thereto.

[0122] The memory 9140 can be a solid-state memory such as a read only memory (ROM), a random access memory (RAM), a SIM card, and the like. It can also be a memory that retains information even when power is off, can be selectively erased, and is provided with more data, an example of which is sometimes referred to as an EPROM or the like. The memory 9140 can also be some other type of device. The memory 9140 includes a buffer memory 9141 (sometimes referred to as a buffer). The memory 9140 can include an application / function storage section 9142 for storing application programs and function programs or a flow for executing the operation of the electronic device 9600 by the central processing unit 9100.

[0123] The memory 9140 can also include a data storage section 9143 for storing data such as contacts, digital data, pictures, sounds, and / or any other data used by the electronic device. A driver storage section 9144 of the memory 9140 can include various drivers of the electronic device for a communication function and / or for executing other functions of the electronic device such as a messaging application, an address book application, and the like.

[0124] The communication module 9110 is a transmitter / receiver 9110 that transmits and receives signals via an antenna 9111. The communication module (transmitter / receiver) 9110 is coupled to the central processing unit 9100 to provide input signals and receive output signals, which can be the same as in the case of a conventional mobile communication terminal.

[0125] Based on different communication technologies, a plurality of communication modules 9110 can be provided in the same electronic device, such as a cellular network module, a Bluetooth module, and / or a wireless local area network module, and the like. The communication module (transmitter / receiver) 9110 is also coupled to a speaker 9131 and a microphone 9132 via an audio processor 9130 to provide audio output via the speaker 9131 and receive audio input from the microphone 9132, thereby implementing a conventional telecommunication function. The audio processor 9130 can include any suitable buffer, decoder, amplifier, and the like. In addition, the audio processor 9130 is also coupled to the central processing unit 9100, thereby enabling recording on the local device through the microphone 9132 and enabling playing of a sound stored on the local device through the speaker 9131.

[0126] The embodiment of the present application also provides a computer readable storage medium capable of realizing all steps of the method for determining the fracture position based on gravity data in the above embodiment, and the computer program is stored on the computer readable storage medium and realizes all steps of the method for determining the fracture position based on gravity data in the above embodiment when the processor executes the computer program, for example, the processor realizes the following steps when executing the computer program:

[0127] The basic data of the gravity data is Bouguer gravity anomaly data, the residual gravity anomaly gridding data is obtained from the Bouguer gravity anomaly gridding data, the pseudo-magnetic anomaly conversion is performed on the residual gravity anomaly gridding data to obtain the pseudo-magnetic anomaly gridding data, the horizontal derivative calculation is performed on the pseudo-magnetic anomaly gridding data to obtain the pseudo-magnetic anomaly horizontal derivative gridding data, the pseudo-magnetic anomaly horizontal derivative graph is drawn according to the pseudo-magnetic anomaly horizontal derivative gridding data, the maximum value positions are identified on the pseudo-magnetic anomaly horizontal derivative graph, and the fracture position is obtained by connecting the maximum value positions.

[0128] As can be seen from the above description, the computer readable storage medium provided by the embodiment of the present application can obtain the Bouguer gravity anomaly data as the basic data of the gravity data, obtain the residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data, perform the pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data to obtain the pseudo-magnetic anomaly gridding data, perform the horizontal derivative calculation on the pseudo-magnetic anomaly gridding data to obtain the pseudo-magnetic anomaly horizontal derivative gridding data, draw the pseudo-magnetic anomaly horizontal derivative graph according to the pseudo-magnetic anomaly horizontal derivative gridding data, identify the maximum value positions on the pseudo-magnetic anomaly horizontal derivative graph, and obtain the fracture position by connecting the maximum value positions, so that the horizontal gradient maximum value strip width of the fracture studied by using the gravity data is reduced by half, the determination accuracy of the fracture position is improved, and the identification ability of the fracture with small scale and large depth is improved.

[0129] Although the present application provides the method operation steps as described in the embodiments or flowcharts, more or less operation steps can be included based on conventional or non-inventive labor. The step sequence listed in the embodiments is only one of the many step execution sequences, and does not represent the only execution sequence. When the device or client product is executed in practice, the method sequence shown in the embodiments or the drawings can be executed in sequence or in parallel (for example, in a parallel processor or multi-thread processing environment).

[0130] Those skilled in the art will appreciate that embodiments of the present specification can be readily used as a method, apparatus (system) or computer program product. Accordingly, embodiments of the present specification can take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROMs, optical storage devices, etc.) embodying computer readable program code thereon for use by or in connection with an instruction execution system. Program Code

[0131] The present invention is described in reference to the flowchart and / or block diagrams of the method, apparatus (system) and computer program product according to embodiments of the present invention. It will be understood that each block of the flowchart and / or block diagrams, and combinations of blocks in the flowchart and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks.

[0132] These computer program instructions can also be stored in a computer- readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instructions which implement the function specified in the flowchart and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks.

[0133] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks.

[0134] Each of the embodiments in the specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other. Each of the embodiments focuses on the difference from other embodiments. In particular, the system embodiments are described simply because they are basically similar to the method embodiments, and the relevant parts can be referred to the part of the method embodiments. In this document, the relationship terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between the entities or operations. It should be noted that the embodiments and features in the embodiments can be combined with each other without conflict. The present application is not limited to any single aspect or embodiment and is not limited to any single aspect or embodiment described, but can include combinations of aspects and / or embodiments.

[0135] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, and not to limit them. Although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that the technical solutions recorded in the foregoing embodiments can still be modified, or some or all of the technical features can be replaced by equivalents. The modification or replacement does not change the essence of the corresponding technical solution, which should be covered in the scope of the claims and the specification of the present application.

Claims

1. A method for determining the location of a fault based on gravity data, characterized in that, The method comprises the following steps: obtaining Bouguer gravity anomaly data as basic data of gravity data; obtaining residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data, and performing pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data to obtain pseudo-magnetic anomaly gridding data; performing horizontal derivative calculation on the pseudo-magnetic anomaly gridding data to obtain pseudo-magnetic anomaly horizontal derivative gridding data; drawing a pseudo-magnetic anomaly horizontal derivative graph according to the pseudo-magnetic anomaly horizontal derivative gridding data; identifying maximum value positions on the pseudo-magnetic anomaly horizontal derivative graph, and connecting the maximum value positions to obtain a fracture position; the pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data comprises: performing pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data by using Poisson formula of magnetic potential and gravitational potential of the same uniform geological body.

2. The method of claim 1, wherein, After the step of obtaining the Bouguer gravity anomaly data as the basic data of the gravity data, the method further comprises the following steps: performing gridding interpolation processing on the Bouguer gravity anomaly data to obtain Bouguer gravity anomaly gridding data; determining that there is a regional background gravity field in the Bouguer gravity anomaly gridding data, and then obtaining the residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data.

3. The method of claim 1, wherein, the step of obtaining the residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data comprises: obtaining the residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data by using a potential field continuation method or a wave number domain filtering method.

4. An apparatus for determining the location of a fracture based on gravity data, comprising: The method comprises the following steps: an acquisition unit is configured to obtain Bouguer gravity anomaly data as basic data of gravity data; an extraction unit is configured to obtain residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data, and perform pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data to obtain pseudo-magnetic anomaly gridding data; a calculation unit is configured to perform horizontal derivative calculation on the pseudo-magnetic anomaly gridding data to obtain pseudo-magnetic anomaly horizontal derivative gridding data; a drawing unit is configured to draw a pseudo-magnetic anomaly horizontal derivative graph according to the pseudo-magnetic anomaly horizontal derivative gridding data; a determination unit is configured to identify maximum value positions on the pseudo-magnetic anomaly horizontal derivative graph, and connect the maximum value positions to obtain a fracture position. the extraction unit comprises: a second extraction subunit is configured to perform pseudo-magnetic anomaly conversion on the residual gravity anomaly gridding data by using Poisson formula of magnetic potential and gravitational potential of the same uniform geological body.

5. The apparatus for determining a fracture location based on gravity data according to claim 4, wherein, The method further comprises the following steps: an interpolation unit is configured to perform gridding interpolation processing on the Bouguer gravity anomaly data to obtain Bouguer gravity anomaly gridding data; a gravity field unit is configured to determine that there is a regional background gravity field in the Bouguer gravity anomaly gridding data, and then obtain the residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data.

6. The apparatus for determining a fracture location based on gravity data according to claim 4, wherein, the extraction unit comprises: a first extraction subunit is configured to obtain the residual gravity anomaly gridding data from the Bouguer gravity anomaly gridding data by using a potential field continuation method or a wave number domain filtering method.

7. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor implements the steps of the method for determining a fracture position based on gravity data according to any one of claims 1 to 3 when executing the program.

8. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program, which is executed by a processor, realizes the steps of the method for determining the fracture location based on gravity data according to any one of claims 1 to 3.

Images