A method and device for quantitatively identifying marine magnetic anomaly stripes

Abstract

The present invention discloses a method and device for quantitatively identifying marine magnetic anomaly stripes. The method includes: interpolating each positive-polarity anomaly and negative-polarity anomaly of the actually observed marine magnetic anomaly stripes and the forward-modeled marine magnetic anomaly stripes into data points with equal numbers; dividing the data points of each positive-polarity anomaly and negative-polarity anomaly into blocks with the same number and calculating the area of each block; sliding the forward-modeled marine magnetic anomaly stripes along the actually observed marine magnetic anomaly stripes, calculating the cosine similarity between the forward-modeled magnetic anomaly stripes and the observed magnetic anomaly stripes according to the area of each block, and the position corresponding to the maximum value of the cosine similarity is the position of the identified magnetic anomaly stripes. The advantages of the present invention are: realizing the quantitative and automatic identification of marine magnetic anomaly stripes, solving the problem that the existing marine magnetic anomaly stripes are only identified by the human eye, which is subjective and cannot provide a quantitative judgment basis for the identification results.

Description

Technical Field

[0001] The present invention relates to the field of marine magnetics, and more particularly to a method and device for quantitatively identifying marine magnetic anomaly stripes. Background Art

[0002] Marine magnetic anomaly stripes are formed by the thermoremanent magnetization recording the geomagnetic field intensity at that time when the magma near the mid-ocean ridge cools down to the Curie point of the internal ferromagnetic substances during the sea-floor spreading process. They carry information such as the age of the oceanic crust and the change of the paleogeomagnetic field intensity.

[0003] Marine magnetic anomaly stripes are often used to determine the spreading age of the ocean basin. The method is to compare the forward-simulated marine magnetic anomaly stripes with the actually observed marine magnetic anomaly stripes by the human eye, and observe the similarity between the two. If the similarity between the actually observed marine magnetic anomaly stripes and the forward-simulated marine magnetic anomaly stripes is high, it is considered that the ages of the actually observed marine magnetic anomaly stripes and the forward-simulated marine magnetic anomaly stripes are the same. This method of relying on the human eye to observe and compare to identify magnetic anomaly stripes has been used until now.

[0004] However, the actually observed marine magnetic anomaly stripes are often affected by seamounts, transform faults, late magmatic activities, etc., and their shapes become irregular. By observing and comparing with the human eye, a quantitative evaluation of the similarity between the actually observed magnetic anomaly stripes and the forward-simulated magnetic anomaly stripes cannot be given, and subjectivity is inevitable in the identification process. Summary of the Invention

[0005] The technical problem to be solved by the present invention is that the existing technology for identifying marine magnetic anomaly stripes mainly uses the human eye to observe and compare, and a quantitative evaluation of the similarity between the actually observed magnetic anomaly stripes and the forward-simulated magnetic anomaly stripes cannot be given, and subjectivity is inevitable in the identification process.

[0006] The present invention solves the above technical problem by the following technical means: A method for quantitatively identifying marine magnetic anomaly stripes, the method comprising:

[0007] Step 1: Interpolate each positive-polarity anomaly and negative-polarity anomaly of the actually observed marine magnetic anomaly stripes and the forward-simulated marine magnetic anomaly stripes into data points with equal numbers.

[0008] Step 2: Divide each positive-polarity anomaly and negative-polarity anomaly into the same number of blocks along the data points and calculate the area of each block.

[0009] Step 3: Slide the forward-simulated marine magnetic anomaly stripes along the actually observed marine magnetic anomaly stripes, calculate the cosine similarity between the forward-simulated magnetic anomaly stripes and the observed magnetic anomaly stripes according to the area of each block, and the position corresponding to the maximum value of the cosine similarity is the position of the identified magnetic anomaly stripes.

[0010] The present invention calculates the cosine similarity between the forward magnetic anomaly strip and the observed magnetic anomaly strip, and the position corresponding to the maximum value of the cosine similarity is the identified position of the magnetic anomaly strip, realizing the quantitative automatic identification of marine magnetic anomaly strips, and solving the problem that the existing marine magnetic anomaly strips are only identified by the human eye, which is subjective and cannot provide a quantitative judgment basis for the identification results.

[0011] Further, the first step includes:

[0012] Both the actually observed marine magnetic anomaly strip and the forward marine magnetic anomaly strip include positive polarity anomalies and negative polarity anomalies. Each positive polarity anomaly and negative polarity anomaly of the actually observed marine magnetic anomaly strip and the forward marine magnetic anomaly strip are interpolated into the same number of data points by cubic spline interpolation, and the number of data points after interpolation of positive polarity anomalies and negative polarity anomalies in different years is equal.

[0013] Further, the second step includes:

[0014] Each positive polarity anomaly and negative polarity anomaly are divided into the same number of blocks, and the area of each block is obtained by integrating the internal area of the block

[0015]

[0016] where f(x) is the magnetic anomaly amplitude, k is the number of subdivided regions inside the block, x0 is the starting point coordinate of each block, x k is the ending point coordinate of each block, and x j is the interpolation subdivision point coordinate inside the block.

[0017] Further, the third step includes:

[0018] The forward marine magnetic anomaly strip is slid along the actually observed marine magnetic anomaly strip, and the sliding step size each time is the number n of blocks after the subdivision of the positive polarity anomaly or negative polarity anomaly. After each sliding, the positive polarity anomaly of the forward magnetic anomaly strip exactly corresponds to the positive polarity anomaly or negative polarity anomaly of the actually observed magnetic anomaly strip.

[0019] Further, the third step also includes:

[0020] By the formula calculate the cosine similarity between each positive polarity anomaly and negative polarity anomaly between the forward magnetic anomaly strip and the actually observed magnetic anomaly strip;

[0021] where P s ={S p1 , S p2 ,..., S pn} is the set of areas of each block after the positive or negative polarity anomalies of the marine magnetic anomaly strip in the forward modeling are dissected, Q s ={S q1 , S q2 ,..., S qn} is the set of areas of each block after the positive or negative polarity anomalies of the actually observed marine magnetic anomaly strip are dissected. are the average values of P s and Q s respectively, n is the number of blocks divided for each positive or negative polarity anomaly, and the position of the maximum similarity is the position of the identified magnetic anomaly strip.

[0022] The present invention also provides a device for quantitatively identifying marine magnetic anomaly strips, and the device includes:

[0023] A data point interpolation module, configured to interpolate each positive and negative polarity anomaly of the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip into an equal number of data points;

[0024] A block division module, configured to divide the data points of each positive and negative polarity anomaly into the same number of blocks and calculate the area of each block;

[0025] A magnetic anomaly strip position identification module, configured to slide the forward modeled marine magnetic anomaly strip along the actually observed marine magnetic anomaly strip, calculate the cosine similarity between the forward modeled magnetic anomaly strip and the observed magnetic anomaly strip according to the area of each block, and the position corresponding to the maximum cosine similarity is the position of the identified magnetic anomaly strip.

[0026] Further, the data point interpolation module is further configured to:

[0027] Both the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip include positive and negative polarity anomalies. Each positive and negative polarity anomaly of the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip is interpolated into the same number of data points by cubic spline interpolation, and the number of data points after interpolation of positive and negative polarity anomalies in different years is equal.

[0028] Further, the block division module is further configured to:

[0029] Divide each positive and negative polarity anomaly into the same number of blocks, and integrate the internal area of the block to obtain the area of each block

[0030]

[0031] where f(x) is the magnetic anomaly amplitude, k is the number of dissected regions inside the block, x0 is the starting point coordinate of each block, x kThe coordinate of the termination point of each block, x j The coordinates of the interpolation and dissection points inside the block.

[0032] Furthermore, the magnetic anomaly strip position recognition module is also used for:

[0033] Sliding the forward modeled marine magnetic anomaly strip along the actually observed marine magnetic anomaly strip, with the step size for each slide being the number of blocks n after dissection of the positive or negative polarity anomaly. After each slide, the positive polarity anomaly of the forward modeled magnetic anomaly strip exactly corresponds to the positive or negative polarity anomaly of the actually observed magnetic anomaly strip.

[0034] Furthermore, the magnetic anomaly strip position recognition module is also used for:

[0035] By the formula Calculate the cosine similarity between each positive and negative polarity anomaly between the forward modeled magnetic anomaly strip and the actually observed magnetic anomaly strip;

[0036] Where, P s ={S p1 , S p2 ,..., S pn} is the area set of each block after dissection of the positive or negative polarity anomaly of the forward modeled marine magnetic anomaly strip, Q s ={S q1 , S q2 ,..., S qn} is the area set of each block after dissection of the positive or negative polarity anomaly of the actually observed marine magnetic anomaly strip, are the average values of P s and Q s respectively, n is the number of blocks divided for each positive or negative polarity anomaly, and the position of the maximum similarity value is the recognized position of the magnetic anomaly strip.

[0037] The advantages of the present invention are as follows: The present invention calculates the cosine similarity between the forward modeled magnetic anomaly strip and the observed magnetic anomaly strip, and the position corresponding to the maximum cosine similarity value is the recognized position of the magnetic anomaly strip, realizing the quantitative and automatic recognition of marine magnetic anomaly strips, and solving the problem that the existing marine magnetic anomaly strips are only recognized by the human eye, which is subjective and cannot provide a quantitative judgment basis for the recognition results. Description of the Drawings

[0038] Figure 1 It is a schematic diagram of the dissection of the actually observed magnetic anomaly strip and the forward modeled magnetic anomaly strip in a method for quantitatively recognizing marine magnetic anomaly strips disclosed in Embodiment 1 of the present invention;

[0039] Figure 2Schematic diagram of the forward magnetic anomaly stripe sliding recognition process in a method for quantitatively identifying marine magnetic anomaly stripes disclosed in Embodiment 1 of the present invention;

[0040] Figure 3 Schematic diagram of the forward marine magnetic anomaly stripes in the southwestern Pacific Ocean in a method for quantitatively identifying marine magnetic anomaly stripes disclosed in Embodiment 1 of the present invention;

[0041] Figure 4 Schematic diagram of the actually observed marine magnetic anomaly stripes in the southwestern Pacific Ocean in a method for quantitatively identifying marine magnetic anomaly stripes disclosed in Embodiment 1 of the present invention;

[0042] Figure 5 Recognition result of C32 when identifying the polarity of the actually observed magnetic anomaly stripes in the southwestern Pacific Ocean in a method for quantitatively identifying marine magnetic anomaly stripes disclosed in Embodiment 1 of the present invention;

[0043] Figure 6 Recognition result of C31 - 30 when identifying the polarity of the actually observed magnetic anomaly stripes in the southwestern Pacific Ocean in a method for quantitatively identifying marine magnetic anomaly stripes disclosed in Embodiment 1 of the present invention;

[0044] Figure 7 Recognition result of C29 - 27 when identifying the polarity of the actually observed magnetic anomaly stripes in the southwestern Pacific Ocean in a method for quantitatively identifying marine magnetic anomaly stripes disclosed in Embodiment 1 of the present invention. Detailed implementation manner

[0045] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0046] Embodiment 1

[0047] A method for quantitatively identifying marine magnetic anomaly stripes, the method includes:

[0048] S1: Interpolate each positive - polarity anomaly and negative - polarity anomaly of the actually observed marine magnetic anomaly stripes and the forward marine magnetic anomaly stripes into data points with equal numbers; the specific process is as follows:

[0049] As Figure 1The figure shows a schematic diagram of the dissection of the actually observed magnetic anomaly strip and the forward modeled magnetic anomaly strip in an embodiment of the present invention. Both the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip include positive polarity anomalies and negative polarity anomalies. Each positive polarity anomaly and negative polarity anomaly of the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip are interpolated into the same number of data points using cubic splines, and the number of data points after interpolation of the positive polarity anomalies and negative polarity anomalies in different years is equal.

[0050] S2: Divide the data points of each positive polarity anomaly and negative polarity anomaly into the same number of blocks and calculate the area of each block; the specific process is as follows:

[0051] Divide each positive polarity anomaly and negative polarity anomaly into the same number of blocks, and integrate the internal area of the blocks to obtain the area of each block

[0052]

[0053] where f(x) is the magnetic anomaly amplitude, k is the number of dissection regions inside the block, x0 is the starting point coordinate of each block, x k is the ending point coordinate of each block, and x j is the interpolation dissection point coordinate inside the block.

[0054] S3: Slide the forward modeled marine magnetic anomaly strip along the actually observed marine magnetic anomaly strip, and calculate the cosine similarity between the forward modeled magnetic anomaly strip and the observed magnetic anomaly strip according to the area of each block. The position corresponding to the maximum value of the cosine similarity is the identified position of the magnetic anomaly strip; the specific process is as follows:

[0055] As Figure 2 shown, it is a schematic diagram of the sliding identification process of the forward modeled magnetic anomaly strip in an embodiment of the present invention. Slide the forward modeled marine magnetic anomaly strip along the actually observed marine magnetic anomaly strip, and the step size of each slide is the number n of blocks after the dissection of the positive polarity anomaly or negative polarity anomaly. After each slide, the positive polarity anomaly of the forward modeled magnetic anomaly strip exactly corresponds to the positive polarity anomaly or negative polarity anomaly of the actually observed magnetic anomaly strip.

[0056] Calculate the cosine similarity between each positive polarity anomaly and negative polarity anomaly between the forward modeled magnetic anomaly strip and the actually observed magnetic anomaly strip through the formula ;

[0057] where P s ={S p1 , S p2 ,..., S pn} is the set of areas of each block after the dissection of the positive polarity anomaly or negative polarity anomaly of the forward modeled marine magnetic anomaly strip, and Q s ={S q1 , Sq2 ,...,S qn} is the area set of each block after the positive or negative polarity anomalies of the actually observed marine magnetic anomaly strips are dissected. are the average values of P s and Q s respectively, n is the number of blocks divided for each positive or negative polarity anomaly, and the similarity between the forward marine magnetic anomaly strip and the actually observed marine magnetic anomaly strip is the average value of the similarities of each positive and negative polarity anomaly. The position of the maximum similarity value is the position of the identified magnetic anomaly strip.

[0058] As Figure 3 shown, it is the forward marine magnetic anomaly strip in the southwestern Pacific Ocean in the embodiment of the present invention.

[0059] The forward marine magnetic anomaly strip can be calculated by establishing a magnetization intensity model according to the geomagnetic polarity time scale and using the magnetic anomalies of the polyhedron. The polarity time of the forward marine magnetic anomaly strip in the southwestern Pacific Ocean is C32 - C27. The forward marine magnetic anomaly strip is divided into three polarity time segments, namely polarity time C32, polarity time C31 - 30, and polarity time C29 - 27. The dashed line is the boundary of different polarity time segments.

[0060] See Figure 4 the actually observed marine magnetic anomaly strip in the southwestern Pacific Ocean in the embodiment of the present invention. The actually observed marine magnetic anomaly strip in the southwestern Pacific Ocean is polarity time C32 - C27. The dashed line is the boundary of different polarity time segments.

[0061] See Figure 5 the identification result of the polarity time C32 of the actually observed magnetic anomaly strip in the southwestern Pacific Ocean in the embodiment of the present invention. When the forward marine magnetic anomaly strip with polarity time C32 slides past the actually observed marine magnetic anomaly strip, the similarities between the forward magnetic anomaly strip and the actually observed magnetic anomaly strip at different polarity times are calculated. When the forward magnetic anomaly strip with polarity time C32 passes above the actually observed magnetic anomaly strip with polarity time C32, the similarity between the forward magnetic anomaly strip and the actually observed magnetic anomaly strip is the highest, and the similarity value is 0.8, correctly identifying the polarity time C32.

[0062] See Figure 6In the embodiment of the present invention, the recognition result of C31-30 when observing the polarity of the actual magnetic anomaly strip in the southwestern Pacific Ocean. When the forward magnetic anomaly strip C31-30 slides past the actually observed ocean magnetic anomaly strip during polarity determination, the similarity between the forward magnetic anomaly strip and the actually observed magnetic anomaly strip at different polarities is calculated. When the forward magnetic anomaly strip C31-30 passes above the actually observed magnetic anomaly strip C31-30 during polarity determination, the similarity between the forward magnetic anomaly strip and the actually observed magnetic anomaly strip is the highest, with a similarity value of 0.57, correctly identifying the polarity of C31-30.

[0063] See Figure 7 In the embodiment of the present invention, the recognition result of C29-27 when observing the polarity of the actual magnetic anomaly strip in the southwestern Pacific Ocean. When the forward magnetic anomaly strip C29-27 slides past the actually observed ocean magnetic anomaly strip during polarity determination, the similarity between the forward magnetic anomaly strip and the actually observed magnetic anomaly strip at different polarities is calculated. When the forward magnetic anomaly strip C29-27 passes above the actually observed magnetic anomaly strip C29-27 during polarity determination, the similarity between the forward magnetic anomaly strip and the actually observed magnetic anomaly strip is the highest, with a similarity value of 0.86, correctly identifying the polarity of C29-27.

[0064] Through the above technical solutions, the present invention provides a method for quantitatively and automatically identifying ocean magnetic anomaly strips. Each positive and negative polarity anomaly of the actually observed ocean magnetic anomaly strip and the forward ocean magnetic anomaly strip is divided into an equal number of data points by cubic spline interpolation. Each positive and negative polarity anomaly is divided into the same number of blocks, and the area of each block is obtained by integrating the internal area of the block. The forward ocean magnetic anomaly strip is slid along the actually observed ocean magnetic anomaly strip, and the cosine similarity between the actually observed ocean magnetic anomaly strip and the forward ocean magnetic anomaly strip is calculated. According to the magnitude of the cosine similarity, the similarity value between the forward magnetic anomaly strip and the actual magnetic anomaly strip is given. The position of the maximum similarity value is the position of the identified magnetic anomaly strip. This solves the problem that it is difficult to avoid subjectivity when using the human eye to observe and compare for the identification of ocean magnetic anomaly strips and there is no quantitative similarity judgment basis.

[0065] Embodiment 2

[0066] Based on Embodiment 1, Embodiment 2 of the present invention further provides a device for quantitatively identifying ocean magnetic anomaly strips, and the device includes:

[0067] A data point interpolation module, configured to interpolate each positive and negative polarity anomaly of the actually observed ocean magnetic anomaly strip and the forward ocean magnetic anomaly strip into an equal number of data points;

[0068] The block division module is used to divide the data points of each positive-polarity anomaly and negative-polarity anomaly into the same number of blocks and calculate the area of each block;

[0069] The magnetic anomaly stripe position recognition module is used to slide the forward-modeled marine magnetic anomaly stripe along the actually observed marine magnetic anomaly stripe, calculate the cosine similarity between the forward-modeled magnetic anomaly stripe and the observed magnetic anomaly stripe according to the area of each block, and the position corresponding to the maximum cosine similarity is the recognized magnetic anomaly stripe position.

[0070] Specifically, the data point interpolation module is further used for:

[0071] Both the actually observed marine magnetic anomaly stripe and the forward-modeled marine magnetic anomaly stripe include positive-polarity anomalies and negative-polarity anomalies. Each positive-polarity anomaly and negative-polarity anomaly of the actually observed marine magnetic anomaly stripe and the forward-modeled marine magnetic anomaly stripe are interpolated into the same number of data points by cubic spline interpolation, and the number of data points after interpolation of positive-polarity anomalies and negative-polarity anomalies in different years is equal.

[0072] Specifically, the block division module is further used for:

[0073] Divide each positive-polarity anomaly and negative-polarity anomaly into the same number of blocks, and integrate the internal area of the block to obtain the area of each block

[0074]

[0075] Among them, f(x) is the magnetic anomaly amplitude, k is the number of internal subdivision regions of the block, x0 is the starting point coordinate of each block, and x k is the ending point coordinate of each block, and x j is the interpolation subdivision point coordinate inside the block.

[0076] Specifically, the magnetic anomaly stripe position recognition module is further used for:

[0077] Slide the forward-modeled marine magnetic anomaly stripe along the actually observed marine magnetic anomaly stripe. The step size of each slide is the number of blocks n after subdivision of the positive-polarity anomaly or negative-polarity anomaly. After each slide, the positive-polarity anomaly of the forward-modeled magnetic anomaly stripe exactly corresponds to the positive-polarity anomaly or negative-polarity anomaly of the actually observed magnetic anomaly stripe.

[0078] Specifically, the magnetic anomaly stripe position recognition module is further used for:

[0079] Through the formula Calculate the cosine similarity between each positive-polarity anomaly and negative-polarity anomaly between the forward-modeled magnetic anomaly stripe and the actually observed magnetic anomaly stripe;

[0080] Among them, P s ={Sp1 , S p2 ,..., S pn}, which is the area set of each block after the positive or negative polarity anomalies of the marine magnetic anomaly strips in the forward modeling are dissected. Q s = {S q1 , S q2 ,..., S qn}, which is the area set of each block after the positive or negative polarity anomalies of the actually observed marine magnetic anomaly strips are dissected. are respectively the averages of P s and Q s . n is the number of blocks divided for each positive or negative polarity anomaly. The position of the maximum similarity is the position of the identified magnetic anomaly strip.

[0081] The above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements on some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for quantitatively identifying marine magnetic anomaly stripes, characterized in that, The method includes: Step 1: Interpolate each positive and negative polarity anomaly of the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip into an equal number of data points. Step 2: Divide each positive and negative polarity anomaly along the data points into the same number of blocks and calculate the area of each block. Step 3: Slide the forward modeled marine magnetic anomaly strip along the actually observed marine magnetic anomaly strip. The sliding step size each time is the number n of blocks after the dissection of the positive or negative polarity anomaly. After each slide, the positive polarity anomaly of the forward modeled magnetic anomaly strip exactly corresponds to the positive or negative polarity anomaly of the actually observed magnetic anomaly strip. Calculate the cosine similarity between the forward modeled magnetic anomaly strip and the observed magnetic anomaly strip according to the area of each block. The position corresponding to the maximum value of the cosine similarity is the identified position of the magnetic anomaly strip.

2. The quantitative identification method of marine magnetic anomaly stripes according to claim 1, characterized in that The said Step 1 includes: Both the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip include positive and negative polarity anomalies. Interpolate each positive and negative polarity anomaly of the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip into the same number of data points using cubic spline interpolation. The number of data points after interpolation of positive and negative polarity anomalies in different years is equal.

3. The quantitative identification method of marine magnetic anomaly stripes according to claim 1, characterized in that The said Step 2 includes: Divide each positive and negative polarity anomaly into the same number of blocks, and integrate the internal area of the block to obtain the area of each block. Among them, f(x) is the magnetic anomaly amplitude, k is the number of internal dissection regions of the block, x0 is the starting point coordinate of each block, and x k is the ending point coordinate of each block, and x j is the interpolation dissection point coordinate inside the block.

4. A method for quantitatively identifying marine magnetic anomaly stripes according to claim 1, characterized in that, The said Step 3 further includes: Through the formula Calculate the cosine similarity between each positive and negative magnetic anomaly between the forward magnetic anomaly strip and the actual observed magnetic anomaly strip; Among them, P s ={S p1 , S p2 ,..., S pn} is the area set of each block after the positive or negative polarity anomaly dissection of the forward marine magnetic anomaly strip. Q s ={S q1 , S q2 ,..., S qn} is the area set of each block after the positive or negative polarity anomaly dissection of the actually observed marine magnetic anomaly strip. are the average values of P s and Q s respectively. n is the number of blocks divided by each positive or negative polarity anomaly, and the position of the maximum similarity value is the position of the identified magnetic anomaly strip.

5. An apparatus for quantitatively identifying marine magnetic anomaly stripes, characterized in that, The device includes: A data point interpolation module, which is used to interpolate each positive and negative polarity anomaly of the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip into an equal number of data points. A block division module, which is used to divide each positive and negative polarity anomaly along the data points into the same number of blocks and calculate the area of each block. A magnetic anomaly strip position identification module, which is used to slide the forward modeled marine magnetic anomaly strip along the actually observed marine magnetic anomaly strip. The sliding step size each time is the number n of blocks after the dissection of the positive or negative polarity anomaly. After each slide, the positive polarity anomaly of the forward modeled magnetic anomaly strip exactly corresponds to the positive or negative polarity anomaly of the actually observed magnetic anomaly strip. Calculate the cosine similarity between the forward modeled magnetic anomaly strip and the observed magnetic anomaly strip according to the area of each block. The position corresponding to the maximum value of the cosine similarity is the identified position of the magnetic anomaly strip.

6. The quantitative identification device for marine magnetic anomaly stripes according to claim 5, characterized in that The said data point interpolation module is also used for: Both the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip include positive and negative polarity anomalies. Interpolate each positive and negative polarity anomaly of the actually observed marine magnetic anomaly strip and the forward modeled marine magnetic anomaly strip into the same number of data points using cubic spline interpolation. The number of data points after interpolation of positive and negative polarity anomalies in different years is equal.

7. The quantitative identification device for marine magnetic anomaly stripes according to claim 5, characterized in that, The said block division module is also used for: Divide each positive and negative polarity anomaly into the same number of blocks, and integrate the internal area of the block to obtain the area of each block. Among them, f(x) is the magnetic anomaly amplitude, k is the number of internal dissection regions of the block, x0 is the starting point coordinate of each block, and x k is the ending point coordinate of each block, and x j is the interpolation dissection point coordinate inside the block.

8. A device for quantitatively identifying marine magnetic anomaly stripes according to claim 5, characterized in that, The said magnetic anomaly strip position identification module is also used for: Through the formula Calculate the cosine similarity between each positive and negative magnetic anomaly between the forward magnetic anomaly strip and the actual observed magnetic anomaly strip; Among them, P s = {S p1 , S p2 ,..., S pn} is the set of areas of each block after the positive or negative polarity anomalies of the marine magnetic anomaly strips in the forward modeling are dissected. Q s = {S q1 , S q2 ,..., S qn} is the set of areas of each block after the positive or negative polarity anomalies of the actually observed marine magnetic anomaly strips are dissected. p and q are the average values of P s and Q s respectively. n is the number of blocks divided for each positive or negative polarity anomaly. The position of the maximum similarity value is the position of the identified magnetic anomaly strip.

Images