Method, device, equipment and storage medium for determining crust shortening

Abstract

The application discloses a method, device and equipment for determining crust shortening amount and a storage medium, and belongs to the technical field of geophysical exploration. The method comprises the following steps: acquiring gravity data and geological data of a detection area comprising a compression structure area and a gravity background area; determining a first gravity anomaly integral amount of the detection area based on the gravity data of the detection area; determining a gravity contrast model corresponding to a specified area based on the geological data of the detection area, and further calculating a second gravity anomaly integral amount of the detection area; and calculating the crust shortening amount of the detection area based on the first gravity anomaly integral amount and the second gravity anomaly integral amount. According to the method, the gravity data is combined with the deformation of the compression structure based on the gravity anomaly of the crust compression structure, so that the crust shortening amount can be obtained simply, quickly, accurately and reliably, and reliable prior information can be provided for compression structure research and geophysical structure modeling.

Description

Technical Field

[0001] This application relates to the field of geophysical exploration technology, and in particular to a method, apparatus, device, and storage medium for determining crustal shortening. Background Technology

[0002] When sedimentary strata in nature are subjected to horizontal compressive stress, they undergo deformations such as folding, uplift, or fracturing. This causes geological bodies within the sedimentary strata to move out of their original positions, resulting in a reduction in the width of the sedimentary strata. The amount of this reduction in width is known as crustal shortening. In the field of geophysical exploration, crustal shortening is an important parameter for stratigraphic research and oil and gas exploration.

[0003] Currently, the amount of crustal shortening is often estimated using the envelope π-value method or the equilibrium profile method. When using this method, it is usually necessary to first identify marker strata in the sedimentary strata, then trace these marker strata using methods such as artificial drilling or establishing stratigraphic models, and finally, express the amount of crustal shortening in the sedimentary strata as the amount of shortening of the marker strata.

[0004] However, in areas with complex stratigraphic structures, it is often difficult to identify and trace marker strata, which makes it impossible for the aforementioned envelope π-value method or equilibrium profile method to obtain the amount of crustal shortening in that area. Summary of the Invention

[0005] This application provides a method, apparatus, device, and storage medium for determining crustal shortening, which solves the problem in related technologies where it is impossible to accurately obtain the crustal shortening in regions with complex crustal structures. The technical solution is as follows:

[0006] Firstly, a method for determining crustal shortening is provided, the method comprising:

[0007] The coordinates, elevation, and gravity values ​​of multiple measuring points in the detection area are obtained. The detection area includes a compression structure area and a gravity background area. The multiple measuring points are linearly distributed at equal intervals and perpendicularly pass through the main body of the compression structure area. The distance between two adjacent measuring points is the first point distance.

[0008] Based on the coordinate elevation and gravity value corresponding to the multiple measuring points, and the first point distance, the first gravity anomaly integral quantity of the detection area is determined;

[0009] The location information of multiple calculation points in the detection area is obtained, along with the horizontal coordinates, depth, width, and thickness of the specified area, the average density of the strata in the specified area, and the average density of the strata in the tectonic background area. The specified area is located within the extrusion tectonic area. The multiple calculation points are linearly distributed at equal intervals and perpendicularly pass through the main body of the specified area. At least three of the calculation points are located on the specified area, and the distance between two adjacent calculation points is the second point distance. The tectonic background area is a region located outside the extrusion tectonic area with a depth equal to the depth of the specified area and a thickness equal to the thickness of the specified area.

[0010] Calculate the density difference between the average density of the strata in the specified area and the average density of the strata in the tectonic background area;

[0011] The horizontal coordinates, depth, width, thickness, and extension direction of the specified region, as well as the density difference, are determined as a gravity comparison model;

[0012] Based on the gravity forward model formula of the horizontally infinitely extending square column model, the gravity comparison model, the position information of the multiple calculation points, and the second point distance, the second gravity anomaly integral quantity of the detection area is determined;

[0013] Based on the first gravity anomaly integral, the second gravity anomaly integral, and the width of the designated area, the crustal shortening of the detection area is determined.

[0014] In one possible implementation, determining the first gravity anomaly integral of the detection area based on the coordinate elevation and gravity value corresponding to the plurality of measuring points, and the first point distance, includes:

[0015] Based on the coordinate elevations corresponding to the multiple measuring points, Bouguer correction, terrain correction, and normal field correction are applied to the gravity values ​​corresponding to each measuring point to obtain the Bouguer gravity anomaly corresponding to each measuring point.

[0016] Based on the calculation of balanced gravity anomaly, filtering method or regional field removal method, and the Bouguer gravity anomaly corresponding to each measuring point, the residual gravity anomaly corresponding to each measuring point is determined, wherein the residual gravity anomaly corresponding to the measuring point located on the gravity background region is 0.

[0017] The first product of the residual gravity anomaly corresponding to each measuring point and the first distance from the first point is determined, and the first product corresponding to each measuring point is summed to obtain the first gravity anomaly integral of the detection area.

[0018] In one possible implementation, determining the second gravity anomaly integral of the detection area based on the gravity forward model of the horizontally infinitely extending prism model, the gravity comparison model, the position information of the multiple calculation points, and the second point distance includes:

[0019] Using the gravity forward modeling formula of a horizontally infinitely extending square column model and the position information of the multiple calculation points, the gravity comparison model is forward modeled at each calculation point to obtain the model gravity anomaly corresponding to each calculation point.

[0020] Determine the second product of the model gravity anomaly corresponding to each calculation point and the second point distance, and sum the second products corresponding to each calculation point to obtain the second gravity anomaly integral of the detection area.

[0021] In one possible implementation, determining the crustal shortening of the detection area based on the first gravity anomaly integral, the second gravity anomaly integral, and the width of the designated area includes:

[0022] Calculate the quotient of the first integral of gravity anomaly and the second integral of gravity anomaly;

[0023] Calculate the third product of the quotient and the width of the specified area, and determine the third product as the crustal shortening of the detection area.

[0024] Secondly, an apparatus for determining crustal shortening is provided, the apparatus comprising:

[0025] The first acquisition module is used to acquire the coordinate elevation and gravity value corresponding to multiple measuring points in the detection area. The detection area includes a compression structure area and a gravity background area. The multiple measuring points are linearly distributed at equal intervals and perpendicularly pass through the main body of the compression structure area. The distance between two adjacent measuring points is the first point distance.

[0026] The first determining module is used to determine the first gravity anomaly integral quantity of the detection area based on the coordinate elevation and gravity value corresponding to the plurality of measuring points and the first point distance;

[0027] The second acquisition module is used to acquire the location information of multiple calculation points in the detection area, the horizontal coordinates, depth, width and thickness of the specified area, the average density of the strata in the specified area and the average density of the strata in the structural background area, wherein the specified area is located within the extrusion structure area, the multiple calculation points are linearly distributed at equal intervals and perpendicularly pass through the main body of the specified area, at least three of the calculation points are located on the specified area, the distance between two adjacent calculation points is the second point distance, and the structural background area is an area located outside the extrusion structure area with a depth equal to the depth of the specified area and a thickness equal to the thickness of the specified area;

[0028] The calculation module is used to calculate the density difference between the average density of the strata in the specified area and the average density of the strata in the tectonic background area;

[0029] The second determining module is used to determine the horizontal coordinates, depth, width, thickness and extension direction of the specified area, as well as the density difference, as a gravity comparison model;

[0030] The third determining module is used to determine the second gravity anomaly integral quantity of the specified region based on the gravity forward model formula of the horizontally infinitely extending square column model, the gravity comparison model, the position information of the multiple calculation points, and the second point distance.

[0031] The fourth determining module is used to determine the crustal shortening of the detection area based on the first gravity anomaly integral, the second gravity anomaly integral, and the width of the designated area.

[0032] In one possible implementation, the first determining module is configured to:

[0033] Based on the coordinate elevations corresponding to the multiple measuring points, Bouguer correction, terrain correction, and normal field correction are applied to the gravity values ​​corresponding to each measuring point to obtain the Bouguer gravity anomaly corresponding to each measuring point.

[0034] Based on the calculation of balanced gravity anomaly, filtering method or regional field removal method, and the Bouguer gravity anomaly corresponding to each measuring point, the residual gravity anomaly corresponding to each measuring point is determined, wherein the residual gravity anomaly corresponding to the measuring point located on the gravity background region is 0.

[0035] The first product of the residual gravity anomaly corresponding to each measuring point and the first distance from the first point is determined, and the first product corresponding to each measuring point is summed to obtain the first gravity anomaly integral of the detection area.

[0036] In one possible implementation, the third determining module is configured to:

[0037] Using the gravity forward modeling formula of a horizontally infinitely extending square column model and the position information of the multiple calculation points, the gravity comparison model is forward modeled at each calculation point to obtain the model gravity anomaly corresponding to each calculation point.

[0038] Determine the second product of the model gravity anomaly corresponding to each calculation point and the second point distance, and sum the second products corresponding to each calculation point to obtain the second gravity anomaly integral of the detection area.

[0039] In one possible implementation, the fourth determining module is configured to:

[0040] Calculate the quotient of the first integral of gravity anomaly and the second integral of gravity anomaly;

[0041] Calculate the third product of the quotient and the width of the specified area, and determine the third product as the crustal shortening of the detection area.

[0042] Thirdly, a computer device is provided, the computer device including a processor and a memory, the memory storing at least one instruction, the instruction being loaded and executed by the processor to perform the operation performed by the method for determining crustal shortening.

[0043] Fourthly, a computer-readable storage medium is provided, wherein at least one instruction is stored in the storage medium, the instruction being loaded and executed by a processor to implement the operations performed by the method for determining crustal shortening.

[0044] Fifthly, a computer program product is provided, comprising computer program code, wherein when the computer program code is executed by a computer device, the computer device executes the method described in the first aspect and its possible implementations.

[0045] The beneficial effects of the technical solutions provided in this application are:

[0046] The scheme mentioned in this application embodiment can first determine the first gravity anomaly integral of the detection area based on gravity data of the detection area; then, based on geological data of the detection area, determine the gravity comparison model corresponding to the specified area, and further calculate the second gravity anomaly integral of the detection area; finally, based on the first gravity anomaly integral, the second gravity anomaly integral, and the width of the specified area, determine the crustal shortening corresponding to the detection area. In related technologies, for areas with complex stratigraphic structures, it is impossible to identify and track marker strata, making it impossible to obtain crustal shortening. However, using this method, based on the gravity anomaly of crustal compression structures, it is not necessary to identify or track marker strata in the compression structure area. It is only necessary to combine gravity data with the deformation of the compression structure to obtain crustal shortening simply, quickly, accurately, and reliably, thus providing reliable prior information for compression structure research and geophysical structural modeling. Attached Figure Description

[0047] To more clearly illustrate the technical solutions in the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0048] Figure 1 This is a schematic diagram of the stratigraphic structure before compression, provided in an embodiment of this application;

[0049] Figure 2 This is a schematic diagram of a stratigraphic structure after compression, provided in an embodiment of this application;

[0050] Figure 3 This is a schematic diagram of the structure of a terminal device provided in an embodiment of this application;

[0051] Figure 4 This is a flowchart illustrating a method for determining crustal shortening according to an embodiment of this application;

[0052] Figure 5 This is a schematic diagram of a gravity anomaly curve provided in an embodiment of this application;

[0053] Figure 6 This is a schematic diagram of a device for determining crustal shortening provided in an embodiment of this application;

[0054] Figure 7 This is a structural block diagram of a computer device provided in an embodiment of this application. Detailed Implementation

[0055] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0056] Most geological bodies in nature undergo deformation after formation due to stress, such as bending, folding, fracturing, and compression, causing them to move out of their original location. Sedimentary strata, when subjected to horizontal compressive stress, will fold, uplift, and fracture, resulting in a reduction in width; the amount of this reduction in stratum width is called crustal shortening. As an example, Figure 1 This is a schematic diagram of the stratigraphic structure before compression, provided in an embodiment of this application. Figure 2 This is a schematic diagram of a stratigraphic structure after compression, provided in an embodiment of this application.

[0057] First, several terms used in the embodiments of this application will be explained:

[0058] Compressional tectonic zone: A region in the Earth's crust where strata undergo significant deformation when subjected to compression.

[0059] Gravity background region: The region outside the extrusion structure where the gravity response of the extrusion structure is close to 0.

[0060] Structural orientation: The direction of the ridgeline of the structure formed after deformation in the compressional tectonic region.

[0061] The method for determining crustal shortening provided in this application embodiment can be executed by a terminal device. From a hardware perspective, such as... Figure 3 As shown, the terminal device may include a processor 310, a memory 320, a communication component 330, and a display component 340, etc.

[0062] The processor 310 can be a CPU (Central Processing Unit) or a SoC (System on Chip), etc. The processor 310 can be used to determine the first gravity anomaly integral in the detection area, to determine the density difference, to determine the gravity comparison model, to determine the second gravity anomaly integral in the detection area, to determine the crustal shortening in the detection area, and so on.

[0063] The memory 320 can be various volatile or non-volatile memories, such as SSD (Solid State Disk) or DRAM (Dynamic Random Access Memory). The memory 320 can be used to store pre-stored data, intermediate data, and result data in the process of determining crustal shortening. For example, it can store the coordinates, elevations, and gravity values ​​corresponding to multiple measuring points, the location information of multiple calculation points, the first and second point distances within the detection area, the average thickness of the strata in a specified area, and so on.

[0064] The communication component 330 can be a wired network connector, a WiFi (Wireless Fidelity) module, a Bluetooth module, a cellular network communication module, etc. The communication component 330 can be used to transmit data with other devices, such as gravimeters, theodolites, etc. The communication component 330 can be used to receive coordinate elevation and gravity values ​​corresponding to multiple measuring points, and can be used to receive location information of multiple calculation points, etc.

[0065] The display component 340 can be a display panel integrated with the terminal device, or it can be a display device separately from the terminal device but with a communication connection. When the display component is a display panel, it can be a TN (Twisted Nematic) panel, a VA (Vertical Alignment) panel, an IPS (In-Plane Switching) panel, etc. The display component 340 can be used to display the detection area, to display multiple measuring points and multiple calculation points in the detection area, to display the prismatic model corresponding to the gravity anomaly comparison data volume, etc.

[0066] This application provides embodiments such as Figure 4 The process flow for determining crustal shortening, as shown, includes the following steps.

[0067] S401. Obtain the coordinates, elevation, and gravity values ​​of multiple measuring points in the detection area.

[0068] After determining the area to be explored (including the extrusion structure area and the gravity background area), technicians can first determine the orientation of the main structure of the extrusion structure area. At the same time, a coordinate system corresponding to the exploration area can be established in the terminal equipment, and the coordinate information of each boundary point of the exploration area in the coordinate system can be determined. Any point in the exploration area can be used as the origin of the coordinate system.

[0069] Technicians can pre-define gravity survey lines for the detection area. These lines pass through the main body of the compression structure region, and both ends extend into the gravity background region. For example,... Figure 2 The gravity survey line O shown in the figure. Moreover, the aforementioned gravity survey line should be perpendicular to the structural trend of the main structure of the extrusion structure region.

[0070] Technicians can pre-set multiple measuring points along the gravity survey line, with the distance between any two adjacent measuring points (called the first point distance) being equal. Typically, the first point distance is less than or equal to 2 kilometers, such as 100 meters, 200 meters, 500 meters, 1 kilometer, 2 kilometers, etc. Technicians can use equipment such as gravimeters to measure the gravity value (in the field of geophysical exploration, this refers to gravitational acceleration) at each measuring point, and use equipment such as positioning instruments or theodolites to measure the corresponding coordinate elevation of each measuring point.

[0071] Finally, the technicians input the coordinates, elevation, and gravity values ​​corresponding to each measuring point through an input device (such as a mouse or keyboard) on the terminal device and send them to the processor 310. The processor 310 then obtains this data for subsequent processing. Optionally, the aforementioned gravimeter, positioning instrument, or theodolite can establish a communication connection with the terminal device. The processor 310 of the terminal device can directly obtain the coordinates, elevation, and gravity values ​​corresponding to each measuring point from the aforementioned instruments. The data transmission process will not be elaborated here.

[0072] S402. Based on the coordinate elevation and gravity values ​​corresponding to multiple measuring points, and the first point distance, determine the first gravity anomaly integral quantity of the detection area.

[0073] After the processor 310 of the terminal device obtains the coordinates, elevation, and gravity value corresponding to each measuring point, it can determine the first gravity anomaly integral in the detection area. The calculation process is as follows:

[0074] First, based on the coordinates and elevations of multiple measuring points, Bouguer correction, terrain correction, and normal field correction are applied to the gravity values ​​corresponding to each measuring point to obtain the Bouguer gravity anomaly corresponding to each measuring point.

[0075] Then, based on the Bouguer gravity anomaly corresponding to each measuring point, methods such as equilibrium gravity anomaly calculation, filtering, or region field removal (or a combination of methods) are used to remove the background gravity field, thereby determining the residual gravity anomaly for each measuring point. When calculating the residual gravity anomaly for each measuring point, it is necessary to ensure that the residual gravity anomaly corresponding to the measuring point in the gravity background region is 0. In other words, it is necessary to ensure that the residual gravity anomaly corresponding to the measuring point far from the specified region in the gravity background region approaches 0. The processor 310 can fit the residual gravity anomaly corresponding to each measuring point into a curve, as an example, such as... Figure 5 The curve Δg is shown in the figure.

[0076] Finally, the processor 310 calculates the first product of the remaining gravity anomaly corresponding to each measuring point and the first point distance in the detection area, and sums the first product corresponding to each measuring point to obtain the first gravity anomaly integral of the detection area.

[0077] S403. Obtain the location information of multiple calculation points in the detection area, the second point distance in the detection area, the horizontal coordinates, depth, width and thickness of the specified area, the average density of the strata in the specified area and the average density of the strata in the structural background area.

[0078] Technicians also need to select the specified area within the extrusion structure area, as an example, such as... Figure 2 The M region is shown in the figure, and data such as the horizontal coordinates, depth, width, thickness, and average density of the strata of the specified region are determined.

[0079] The specified area only needs to be located within the compressional structure region. For example, for this specified area, the horizontal coordinate can be the horizontal coordinate of the midpoint of the portion of the gravity survey line that passes through the compressional structure region; the depth (also known as the center burial depth) can be the average depth of the compressional structure region; the width can be 1 kilometer; the thickness can be the average thickness of the strata involved in the compressional structure region (i.e., the average thickness of the strata); and the average density of the strata can be the average density of the strata involved in the compressional structure region. The average depth is the average of the average depth of the upper surface of the strata involved in the background gravity field and the minimum depth of the upper surface of the strata involved in the compressional structure region.

[0080] The average thickness and average density of the strata can be obtained as follows: In the compressional structural region or its surrounding area, the thickness and density of the strata involved at each exploration point are determined through methods such as surface geological outcrops, well stratigraphic data, or geophysical data inversion. Then, the average thickness and average density of the strata involved at multiple exploration points are calculated, thus yielding the average thickness and average density of the strata involved in the compressional structural region. Optionally, technicians can obtain the average thickness and average density of the strata involved through analysis of pre-collected geological data; this is not limited here.

[0081] The width of the specified area can be set according to the actual calculation needs, such as 100 meters, 200 meters, 500 meters, 1 kilometer, 2 kilometers, etc. For ease of subsequent calculation, 1 kilometer is usually selected.

[0082] Technicians also need to define the distribution area of ​​calculation points along the gravity survey line and determine the distance between two adjacent calculation points (called the second point distance). Within this distribution area, the calculation points are uniformly and linearly distributed according to the second point distance, and at least three of the multiple calculation points are located within the designated area. The second point distance can be set according to the width of the designated area, such as 100 meters, 200 meters, 500 meters, 1 kilometer, etc., and is usually less than or equal to one-third of the width of the designated area.

[0083] After the technicians send the above-mentioned distribution area and the second point distance to the processor 310, the processor 310 can determine the number of calculation points and the location information of each calculation point. Furthermore, it can determine the positional relationship between each calculation point and the specified area to facilitate subsequent calculations.

[0084] After determining the designated area, the technicians also need to determine the corresponding construction background area. For example, the construction background area could be as follows: Figure 2 The region P shown is a tectonic background region located outside the compressional tectonic region, with a depth and thickness equal to the depth and thickness of the aforementioned specified region. Technicians can determine the average density of the strata in this tectonic background region using methods such as surface geological outcrops, well stratigraphic data, or geophysical data inversion; these methods will not be elaborated upon here.

[0085] Finally, the technicians input the relevant data mentioned in step S403 into the terminal device, where the processor 310 acquires the data and performs subsequent processing. The data transmission process will not be elaborated here.

[0086] S404. Calculate the density difference between the average density of the strata in the specified area and the average density of the strata in the tectonic background area.

[0087] S405. Determine the horizontal coordinates, depth, width, thickness, extension direction, and density difference of the specified area as a gravity comparison model.

[0088] The processor 310 can first calculate the density difference between the average density of the strata in the specified area and the average density of the strata in the tectonic background area.

[0089] The processor 310 can establish a gravity comparison model based on the horizontally infinitely extending prism model. The gravity comparison model can be considered as a data volume. This data volume is used to represent a specific horizontally infinitely extending prism established in the detection area based on the specified area. The data volume (i.e., the gravity comparison model) includes the horizontal coordinates, depth, width, thickness, extension direction, and density difference of the specified area. The extension direction is perpendicular to the gravity profile, and the gravity profile is a profile that passes perpendicularly through the specified area.

[0090] S406. Based on the gravity forward model of the horizontally infinitely extending square column model, the gravity comparison model, the position information corresponding to multiple calculation points, and the second point distance, determine the second gravity anomaly integral quantity of the detection area.

[0091] After determining the gravity anomaly comparison data volume, the processor 310 can determine the second gravity anomaly integral in the detection area. The calculation process is as follows:

[0092] First, the processor 310 can use the gravity forward modeling formula for a horizontally infinitely extending prism model and the position information of the multiple calculation points mentioned above to perform gravity forward modeling calculations at each calculation point of the gravity comparison model, obtaining the model gravity anomaly corresponding to each calculation point. The processor 310 can fit the model gravity anomaly corresponding to each calculation point into a curve, as shown in the example below. Figure 5 The curve Δg0 is shown in the figure.

[0093] As an example, the processor 310 can input each parameter and the coordinate information of each calculation point in the gravity comparison model into the gravity forward modeling formula of the horizontally infinitely extending prism model, thereby performing forward modeling calculations to obtain the model gravity anomaly corresponding to each calculation point.

[0094] Then, the processor 310 can calculate the second product of the model gravity anomaly and the second point distance corresponding to each calculation point, and sum the second products corresponding to each calculation point to obtain the second gravity anomaly integral of the above-mentioned detection area, that is, the sum of multiple second products is used as the second gravity anomaly integral of the detection area.

[0095] S407. Based on the first gravity anomaly integral, the second gravity anomaly integral, and the width of the specified area, determine the crustal shortening of the detection area.

[0096] After determining the first gravity anomaly integral quantity and the second gravity anomaly integral quantity of the specified area, the processor 310 can determine the crustal shortening amount corresponding to the detection area. The calculation process is as follows:

[0097] First, the processor 310 can calculate the quotient of the first gravity anomaly integral and the second gravity anomaly integral.

[0098] Then, calculate the third product of the above quotient and the width of the specified area, and determine the third product as the crustal shortening of the detection area.

[0099] As an example, the processor 310 can pre-establish a geological structure model of the exploration area based on one or more data obtained by technicians using methods such as gravity, magnetics, electrical resistivity, and seismic exploration, as well as the corresponding geological data. Then, the processor 310 can adjust the width parameters of the strata involved in the pre-established geological structure model according to the crustal shortening, thereby obtaining a more accurate geological structure model, and thus providing accurate reference for stratigraphic research and oil and gas exploration operations.

[0100] The scheme mentioned in this application embodiment can first determine the first gravity anomaly integral of the detection area based on gravity data of the detection area; then, based on geological data of the detection area, determine the gravity comparison model corresponding to the specified area, and further calculate the second gravity anomaly integral of the detection area; finally, based on the first gravity anomaly integral, the second gravity anomaly integral, and the width of the specified area, determine the crustal shortening corresponding to the detection area. In related technologies, for areas with complex stratigraphic structures, it is impossible to identify and track marker strata, making it impossible to obtain crustal shortening. However, using this method, based on the gravity anomaly of crustal compression structures, it is not necessary to identify or track marker strata in the compression structure area. It is only necessary to combine gravity data with the deformation of the compression structure to obtain crustal shortening simply, quickly, accurately, and reliably, thus providing reliable prior information for compression structure research and geophysical structural modeling.

[0101] This application provides a device for determining crustal shortening, which may be the processor 310 of the terminal device in the above embodiments, such as... Figure 6 As shown, the device includes: an acquisition module 610, a first determination module 620, a calculation module 630, a second determination module 640, a third determination module 650, and a fourth determination module 660.

[0102] The first acquisition module 610 is used to acquire the coordinate elevation and gravity value corresponding to multiple measuring points in the detection area. The detection area includes a compression structure area and a gravity background area. The multiple measuring points are linearly distributed at equal intervals and perpendicularly pass through the main body of the compression structure area. The distance between two adjacent measuring points is the first point distance.

[0103] The first determining module 620 is used to determine the first gravity anomaly integral of the detection area based on the coordinate elevation and gravity value corresponding to the plurality of measuring points and the first point distance;

[0104] The second acquisition module 630 is used to acquire the position information of multiple calculation points in the detection area, the horizontal coordinates, depth, width and thickness of the specified area, the average density of the strata in the specified area and the average density of the strata in the structural background area, wherein the specified area is located within the extrusion structure area, the multiple calculation points are linearly distributed at equal intervals and perpendicularly pass through the main body of the specified area, at least three of the calculation points are located on the specified area, the distance between two adjacent calculation points is the second point distance, and the structural background area is an area located outside the extrusion structure area with a depth equal to the depth of the specified area and a thickness equal to the thickness of the specified area;

[0105] The calculation module 640 is used to calculate the density difference between the average density of the strata in the specified area and the average density of the strata in the tectonic background area.

[0106] The second determining module 650 is used to determine the horizontal coordinates, depth, width, thickness and extension direction of the specified area, as well as the density difference, as a gravity comparison model;

[0107] The third determining module 660 is used to determine the second gravity anomaly integral quantity of the specified region based on the gravity forward model formula of the horizontally infinitely extending square column model, the gravity comparison model, the position information of the multiple calculation points, and the second point distance.

[0108] The fourth determining module 670 is used to determine the crustal shortening of the detection area based on the first gravity anomaly integral, the second gravity anomaly integral, and the width of the designated area.

[0109] In one possible implementation, the first determining module 620 is used for:

[0110] Based on the coordinate elevations corresponding to the multiple measuring points, Bouguer correction, terrain correction, and normal field correction are applied to the gravity values ​​corresponding to each measuring point to obtain the Bouguer gravity anomaly corresponding to each measuring point.

[0111] Based on the calculation of balanced gravity anomaly, filtering method or regional field removal method, and the Bouguer gravity anomaly corresponding to each measuring point, the residual gravity anomaly corresponding to each measuring point is determined, wherein the residual gravity anomaly corresponding to the measuring point located on the gravity background region is 0.

[0112] The first product of the residual gravity anomaly corresponding to each measuring point and the first distance from the first point is determined, and the first product corresponding to each measuring point is summed to obtain the first gravity anomaly integral of the detection area.

[0113] In one possible implementation, the third determining module 660 is used for:

[0114] Using the gravity forward modeling formula of a horizontally infinitely extending square column model and the position information of the multiple calculation points, the gravity comparison model is forward modeled at each calculation point to obtain the model gravity anomaly corresponding to each calculation point.

[0115] Determine the second product of the model gravity anomaly corresponding to each calculation point and the second point distance, and sum the second products corresponding to each calculation point to obtain the second gravity anomaly integral of the detection area.

[0116] In one possible implementation, the fourth determining module 670 is used for:

[0117] Calculate the quotient of the integral of the first gravity anomaly and the integral of the second gravity anomaly;

[0118] Calculate the third product of the quotient and the width of the specified area, and determine the third product as the crustal shortening of the detection area.

[0119] The scheme mentioned in this application embodiment can first determine the first gravity anomaly integral of the detection area based on gravity data of the detection area; then, based on geological data of the detection area, determine the gravity comparison model corresponding to the specified area, and further calculate the second gravity anomaly integral of the detection area; finally, based on the first gravity anomaly integral, the second gravity anomaly integral, and the width of the specified area, determine the crustal shortening corresponding to the detection area. In related technologies, for areas with complex stratigraphic structures, it is impossible to identify and track marker strata, making it impossible to obtain crustal shortening. However, using this method, based on the gravity anomaly of crustal compression structures, it is not necessary to identify or track marker strata in the compression structure area. It is only necessary to combine gravity data with the deformation of the compression structure to obtain crustal shortening simply, quickly, accurately, and reliably, thus providing reliable prior information for compression structure research and geophysical structural modeling.

[0120] It should be noted that the apparatus for determining crustal shortening provided in the above embodiments is only illustrated by the division of the functional modules described above. In practical applications, the functions described above can be assigned to different functional modules as needed, that is, the internal structure of the apparatus can be divided into different functional modules to complete all or part of the functions described above. Furthermore, the apparatus for determining crustal shortening provided in the above embodiments and the method embodiments for determining crustal shortening belong to the same concept; the specific implementation process is detailed in the method embodiments and will not be repeated here.

[0121] This application provides a computer device, which may be the terminal device described in the above embodiments. Figure 7This is a schematic diagram of the computer device 700. The computer device 700 can vary significantly due to differences in configuration or performance. It may include one or more CPUs 710 and one or more memories 720. The memories 720 store at least one instruction, which is loaded and executed by the processor 710 to implement the methods provided in the above-described embodiments. Of course, the computer device may also have wired or wireless network interfaces, a keyboard, and input / output interfaces for input and output. The computer device may also include other components for implementing device functions, which will not be elaborated upon here.

[0122] In an exemplary embodiment, a computer-readable storage medium is also provided, such as a memory including instructions executable by a processor in a terminal to perform the method for determining crustal shortening in the above embodiments. This computer-readable storage medium may be non-transitory. For example, the computer-readable storage medium may be a ROM (Read-Only Memory), RAM (Random Access Memory), CD-ROM, magnetic tape, floppy disk, and optical data storage device, etc.

[0123] It should be noted that the information (including but not limited to user equipment information, user personal information, etc.), data (including but not limited to data used for analysis, stored data, displayed data, etc.) and signals (including but not limited to signals transmitted between user terminals and other devices) involved in this application are all authorized by the user or fully authorized by all parties, and the collection, use and processing of related data must comply with the relevant laws, regulations and standards of the relevant countries and regions.

[0124] Those skilled in the art will understand that all or part of the steps of the above embodiments can be implemented by hardware or by a program instructing related hardware. The program can be stored in a computer-readable storage medium, such as a read-only memory, a disk, or an optical disk.

[0125] The above description is merely an optional embodiment of this application and is not intended to limit this application. Any modifications, equivalent substitutions, improvements, etc., made within the principles of this application should be included within the protection scope of this application.

Claims

1. A method for determining crustal shortening, characterized in that, The method includes: The coordinates, elevation, and gravity values ​​of multiple measuring points in the detection area are obtained. The detection area includes a compression structure area and a gravity background area. The multiple measuring points are linearly distributed at equal intervals and perpendicularly pass through the main body of the compression structure area. The distance between two adjacent measuring points is the first point distance. Based on the coordinate elevation and gravity value corresponding to the multiple measuring points, and the first point distance, the first gravity anomaly integral quantity of the detection area is determined; The location information of multiple calculation points in the detection area is obtained, along with the horizontal coordinates, depth, width, and thickness of the specified area, the average density of the strata in the specified area, and the average density of the strata in the tectonic background area. The specified area is located within the extrusion tectonic area. The multiple calculation points are linearly distributed at equal intervals and perpendicularly pass through the main body of the specified area. At least three of the calculation points are located on the specified area, and the distance between two adjacent calculation points is the second point distance. The tectonic background area is a region located outside the extrusion tectonic area with a depth equal to the depth of the specified area and a thickness equal to the thickness of the specified area. Calculate the density difference between the average density of the strata in the specified area and the average density of the strata in the tectonic background area; The horizontal coordinates, depth, width, thickness, and extension direction of the specified region, as well as the density difference, are determined as a gravity comparison model; Based on the gravity forward model formula of the horizontally infinitely extending square column model, the gravity comparison model, the position information of the multiple calculation points, and the second point distance, the second gravity anomaly integral quantity of the detection area is determined; Calculate the quotient of the first integral of gravity anomaly and the second integral of gravity anomaly; Calculate the third product of the quotient and the width of the specified area, and determine the third product as the crustal shortening of the detection area.

2. The method according to claim 1, characterized in that, The step of determining the first gravity anomaly integral of the detection area based on the coordinate elevation and gravity value corresponding to the multiple measuring points and the first point distance includes: Based on the coordinate elevations corresponding to the multiple measuring points, Bouguer correction, terrain correction, and normal field correction are applied to the gravity values ​​corresponding to each measuring point to obtain the Bouguer gravity anomaly corresponding to each measuring point. Based on the calculation of balanced gravity anomaly, filtering method or regional field removal method, and the Bouguer gravity anomaly corresponding to each measuring point, the residual gravity anomaly corresponding to each measuring point is determined, wherein the residual gravity anomaly corresponding to the measuring point located on the gravity background region is 0. The first product of the residual gravity anomaly corresponding to each measuring point and the first distance from the first point is determined, and the first product corresponding to each measuring point is summed to obtain the first gravity anomaly integral of the detection area.

3. The method according to claim 1, characterized in that, The determination of the second gravity anomaly integral quantity in the detection area based on the gravity forward model of the horizontally infinitely extending prism model, the gravity comparison model, the position information of the multiple calculation points, and the second point distance includes: Using the gravity forward modeling formula of a horizontally infinitely extending square column model and the position information of the multiple calculation points, the gravity comparison model is forward modeled at each calculation point to obtain the model gravity anomaly corresponding to each calculation point. Determine the second product of the model gravity anomaly corresponding to each calculation point and the second point distance, and sum the second products corresponding to each calculation point to obtain the second gravity anomaly integral of the detection area.

4. A device for determining crustal shortening, characterized in that, The device includes: The first acquisition module is used to acquire the coordinate elevation and gravity value corresponding to multiple measuring points in the detection area. The detection area includes a compression structure area and a gravity background area. The multiple measuring points are linearly distributed at equal intervals and perpendicularly pass through the main body of the compression structure area. The distance between two adjacent measuring points is the first point distance. The first determining module is used to determine the first gravity anomaly integral quantity of the detection area based on the coordinate elevation and gravity value corresponding to the plurality of measuring points and the first point distance; The second acquisition module is used to acquire the location information of multiple calculation points in the detection area, the horizontal coordinates, depth, width and thickness of the specified area, the average density of the strata in the specified area and the average density of the strata in the structural background area, wherein the specified area is located within the extrusion structure area, the multiple calculation points are linearly distributed at equal intervals and perpendicularly pass through the main body of the specified area, at least three of the calculation points are located on the specified area, the distance between two adjacent calculation points is the second point distance, and the structural background area is an area located outside the extrusion structure area with a depth equal to the depth of the specified area and a thickness equal to the thickness of the specified area; The calculation module is used to calculate the density difference between the average density of the strata in the specified area and the average density of the strata in the tectonic background area; The second determining module is used to determine the horizontal coordinates, depth, width, thickness and extension direction of the specified area, as well as the density difference, as a gravity comparison model; The third determining module is used to determine the second gravity anomaly integral quantity of the specified region based on the gravity forward model formula of the horizontally infinitely extending square column model, the gravity comparison model, the position information of the multiple calculation points, and the second point distance. The fourth determining module is used to calculate the quotient of the first gravity anomaly integral and the second gravity anomaly integral; calculate the third product of the quotient and the width of the specified area, and determine the third product as the crustal shortening of the detection area.

5. The apparatus according to claim 4, characterized in that, The first determining module is used for: Based on the coordinate elevations corresponding to the multiple measuring points, Bouguer correction, terrain correction, and normal field correction are applied to the gravity values ​​corresponding to each measuring point to obtain the Bouguer gravity anomaly corresponding to each measuring point. Based on the calculation of balanced gravity anomaly, filtering method or regional field removal method, and the Bouguer gravity anomaly corresponding to each measuring point, the residual gravity anomaly corresponding to each measuring point is determined, wherein the residual gravity anomaly corresponding to the measuring point located on the gravity background region is 0. The first product of the residual gravity anomaly corresponding to each measuring point and the first distance from the first point is determined, and the first product corresponding to each measuring point is summed to obtain the first gravity anomaly integral of the detection area.

6. The apparatus according to claim 4, characterized in that, The third determining module is used for: Using the gravity forward modeling formula of a horizontally infinitely extending square column model and the position information of the multiple calculation points, the gravity comparison model is forward modeled at each calculation point to obtain the model gravity anomaly corresponding to each calculation point. Determine the second product of the model gravity anomaly corresponding to each calculation point and the second point distance, and sum the second products corresponding to each calculation point to obtain the second gravity anomaly integral of the detection area.

7. A computer device, characterized in that, The computer device includes a processor and a memory, the memory storing at least one instruction, which is loaded and executed by the processor to perform the operations performed by the method for determining crustal shortening as described in any one of claims 1 to 3.

8. A computer-readable storage medium, characterized in that, The storage medium stores at least one instruction, which is loaded and executed by a processor to perform the operations of the method for determining crustal shortening as described in any one of claims 1 to 3.

Images