Deep sea rare earth shallow stratum profile imaging method, device and equipment and storage medium

Abstract

The invention discloses a deep sea rare earth shallow stratum profile imaging method, device and equipment and a storage medium, and relates to the technical field of data processing, and the method comprises the steps: carrying out the preprocessing of the initial seismic data of a deep sea rare earth shallow stratum profile; performing correction and noise removal on the preprocessed initial seismic data as refined seismic data; and determining a thickness distribution image of the deep sea rare earth shallow stratum profile according to the refined seismic data. According to the method, the initial seismic data is preprocessed, corrected and denoised to obtain the refined seismic data, and then the refined seismic data is utilized to determine the distribution image of the deep sea rare earth shallow stratum profile, so that the exploration accuracy and exploration efficiency of the deep sea rare earth shallow stratum profile are greatly improved.

Description

Technical Field

[0001] This application relates to the field of data processing technologies, and in particular to a deep-sea rare earth shallow stratum profile imaging method, device, equipment, and storage medium. Background Art

[0002] As a new type of marine mineral resource rich in medium- and heavy-rare earth elements, the resource volume of deep-sea rare earth far exceeds the land rare earth reserves and has great potential. At present, the investigation of deep-sea rare earth resources mainly adopts a method combining sediment core sampling and shallow stratum detection. According to relevant research, the deepest high-grade rare earth enrichment layer in the seabed sediment can reach hundreds of meters below the surface. However, the currently commonly used box corer, grab sampler, and gravity piston coring technologies cannot meet the requirements for deep-sea sediment sampling depths exceeding more than ten meters. The obtained cores mostly do not penetrate the bottom boundary of the enrichment layer, and deep-sea drilling requires a professional drilling ship, with high costs and long time consumption. Summary of the Invention

[0003] The main purpose of the embodiments of this application is to propose a deep-sea rare earth shallow stratum profile imaging method, device, equipment, and storage medium to accurately and efficiently determine the spatial distribution of the deep-sea rare earth shallow stratum profile.

[0004] To achieve the above object, on the one hand, the embodiments of this application propose a deep-sea rare earth shallow stratum profile imaging method, and the method includes the following steps:

[0005] Preprocess the initial seismic data of the deep-sea rare earth shallow stratum profile;

[0006] Correct and remove noise from the preprocessed initial seismic data to obtain refined seismic data;

[0007] Determine the thickness distribution image of the deep-sea rare earth shallow stratum profile according to the refined seismic data.

[0008] In some embodiments, the preprocessing of the initial seismic data of the deep-sea rare earth shallow stratum profile includes the following steps:

[0009] Determine the sound velocity of the initial seismic data, the excitation frequency, and the number of segmented record files of each segmented survey line of the deep-sea rare earth shallow stratum profile;

[0010] Reset the serial number of each segmented survey line according to the acquisition record time of each segmented record file;

[0011] Unify the lengths of each segmented record file of each segmented survey line;

[0012] Merge each segmented record file of each segmented survey line according to the serial number of each segmented survey line to obtain the preprocessed initial seismic data.

[0013] In some embodiments, the step of correcting the preprocessed initial seismic data includes the following steps:

[0014] Perform delay correction on the preprocessed initial seismic data so that the corresponding seabed information in the preprocessed initial seismic data is continuous, and the starting value of the file number of the segmented record file in the preprocessed initial seismic data corresponds one-to-one with the reset serial numbers of each segmented survey line;

[0015] Perform amplitude correction on the preprocessed initial seismic data.

[0016] In some embodiments, the step of removing noise from the preprocessed initial seismic data includes the following steps:

[0017] Remove data exceeding a preset amplitude threshold from the seismic traces of the preprocessed initial seismic data;

[0018] Remove data outside the target wave frequency band from the preprocessed initial seismic data.

[0019] In some embodiments, determining the thickness distribution image of the deep-sea rare earth shallow stratum profile according to the refined seismic data includes the following steps:

[0020] Extract amplitude attributes from the refined seismic data;

[0021] Determine the top layer and bottom boundary of the deep-sea rare earth shallow stratum profile according to the amplitude attributes;

[0022] Determine the thickness distribution image of the deep-sea rare earth shallow stratum profile according to the top layer and bottom boundary.

[0023] In some embodiments, determining the thickness distribution image of the deep-sea rare earth shallow stratum profile according to the top layer and bottom boundary includes the following steps:

[0024] Perform time-depth conversion according to the time curves of the top layer and bottom boundary, the deep-sea drilling plan data of the deep-sea rare earth shallow stratum profile, and the empirical velocity around the deep-sea rare earth shallow stratum profile, so as to obtain the vertical thickness distribution image of the deep-sea rare earth shallow stratum profile.

[0025] In some embodiments, the method further includes the following steps:

[0026] Convert the geographical angular coordinates in the refined seismic data into longitude and latitude coordinates as the geographical coordinates of the deep-sea rare earth shallow stratum profile.

[0027] To achieve the above object, on the other hand, an imaging device for the deep-sea rare earth shallow stratum profile is proposed in an embodiment of the present application. The device includes:

[0028] A data preprocessing unit for preprocessing the initial seismic data of the shallow stratum profile of deep-sea rare earths;

[0029] A data refinement unit for correcting and removing noise from the preprocessed initial seismic data to obtain refined seismic data;

[0030] A thickness distribution determination unit for determining the thickness distribution image of the shallow stratum profile of deep-sea rare earths based on the refined seismic data.

[0031] To achieve the above object, on the other hand, an embodiment of the present application provides an electronic device, which includes a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the above-mentioned imaging method for the shallow stratum profile of deep-sea rare earths is implemented.

[0032] To achieve the above object, on the other hand, an embodiment of the present application provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the above-mentioned imaging method for the shallow stratum profile of deep-sea rare earths is implemented.

[0033] The embodiments of the present application at least include the following beneficial effects:

[0034] The present application can preprocess the initial seismic data of the shallow stratum profile of deep-sea rare earths; correct and remove noise from the preprocessed initial seismic data to obtain refined seismic data; determine the thickness distribution image of the shallow stratum profile of deep-sea rare earths based on the refined seismic data. By preprocessing, correcting, and denoising the initial seismic data, refined seismic data is obtained, and then the distribution image of the shallow stratum profile of deep-sea rare earths is determined using the refined seismic data, greatly improving the exploration accuracy and efficiency of the shallow stratum profile of deep-sea rare earths. Description of the Drawings

[0035] To more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings required for the description of the embodiments. Obviously, the following drawings are only some embodiments of the present application. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0036] Figure 1 It is a schematic flowchart of the imaging method for the shallow stratum profile of deep-sea rare earths provided by the embodiments of the present application;

[0037] Figure 2 It is an example flowchart of the imaging method for the shallow stratum profile of deep-sea rare earths provided by the embodiments of the present application;

[0038] Figures 3(a), 3(b), and 3(c) are the original shallow stratum profile segmented survey line data diagrams provided by the embodiments of the present application;

[0039] Figure 4 is the original shallow stratum profile diagram after data merging provided by the embodiments of the present application;

[0040] Figure 5 is the shallow stratum profile diagram after fine processing and delay correction provided by the embodiments of the present application;

[0041] Figures 6(a) and 6(b) are the shallow stratum profile diagrams before and after amplitude correction provided by the embodiments of the present application;

[0042] Figure 7 is the shallow stratum profile diagram after burst noise suppression provided by the embodiments of the present application;

[0043] Figure 8 is the shallow stratum profile diagram after band-pass filtering provided by the embodiments of the present application;

[0044] Figures 9(a) and 9(b) are the shallow stratum profile diagrams before and after fine processing provided by the embodiments of the present application;

[0045] Figure 10 is an example diagram of the instantaneous amplitude attribute of the shallow stratum profile data provided by the embodiments of the present application;

[0046] Figure 11 is a schematic diagram of the top and bottom interfaces of the deep-sea rare earth-rich sedimentary layer in the shallow stratum profile data provided by the embodiments of the present application;

[0047] Figure 12 is the thickness distribution diagram of the deep-sea rare earth-rich sedimentary layer in a single shallow profile survey line provided by the embodiments of the present application;

[0048] Figure 13 is the thickness plane distribution diagram of the deep-sea rare earth-rich sedimentary layer in multiple shallow profile survey lines provided by the embodiments of the present application;

[0049] Figure 14 is the structural schematic diagram of the deep-sea rare earth shallow stratum profile imaging device provided by the embodiments of the present application;

[0050] Figure 15 is the hardware structural schematic diagram of an electronic device provided by the embodiments of the present application. Detailed implementation manners

[0051] In order to make the objectives, technical solutions, and advantages of the present application more clearly understood, the present application will be further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application. When the following description refers to the accompanying drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the embodiments of the present application. They are only examples of devices and methods consistent with some aspects of the embodiments of the present application as detailed in the appended claims.

[0052] It can be understood that the terms "first", "second", etc. used in the present application may be used herein to describe various concepts, but unless otherwise specified, these concepts are not limited by these terms. These terms are only used to distinguish one concept from another. For example, without departing from the scope of the embodiments of the present application, the first information may also be referred to as the second information. Similarly, the second information may also be referred to as the first information. Depending on the context, the words "if", "when" as used herein may be interpreted as "when...", "while...", or "in response to determining".

[0053] The terms "at least one", "multiple", "each", "any one", etc. used in the present application, at least one includes one, two, or more than two, multiple includes two or more than two, each refers to each one of the corresponding multiple, and any one refers to any one of the multiple.

[0054] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled in the technical field to which the present application belongs. The terms used herein are only for the purpose of describing the embodiments of the present application and are not intended to limit the present application.

[0055] Before elaborating on the embodiments of the present application in detail, some related technologies involved in the embodiments of the present application will be described first, as follows:

[0056] In view of the problems of inaccuracy, high cost, and long time consumption in the exploration of deep-sea rare earth shallow strata profiles in the prior art, it is urgent to develop a fast and feasible deep-sea rare earth exploration plan. At present, the shallow strata exploration technology plan (shallow strata profile measurement) has been widely used in the field of marine mid-shallow strata research and is mainly used to study the sediment distribution characteristics of deep-sea basins and seamounts in the field of ocean mineral resources investigation. It is a conventional technology for detecting shallow strata of submarine sediments. However, at present, the use of shallow strata exploration technology in the field of rare earth resources investigation is still in the exploratory stage, and there is still no standardized and effective processing plan for improving the imaging accuracy of deep-sea rare earth shallow strata profile data.

[0057] Terrestrial rare earth ores are mainly light rare earths, and medium and heavy rare earths are still extremely scarce key metal resources. In recent years, geologists have discovered sediments rich in medium and heavy rare earths in the deep sea. Their metallogenic characteristics are significantly different from those of terrestrial ion-adsorbed rare earth ores. However, there has always been controversy about the occurrence state and extraordinary enrichment mechanism of deep-sea rare earths. In addition, there is a lack of rapid exploration and evaluation technologies for deep-sea rare earths. At present, the investigation of deep-sea rare earth resources mainly adopts a method combining sediment core sampling and shallow stratum detection. According to relevant technical research, the deepest high-grade rare earth enrichment layer in submarine sediments can reach more than a hundred meters below the surface. At present, the commonly used box corer, grab sampler, and gravity piston coring technologies cannot meet the sampling requirements of deep-sea rare earth-rich sediments deeper than ten meters. The obtained cores basically do not penetrate the bottom boundary of the enrichment. Deep-sea drilling requires a professional drilling ship, with high costs and low exploration efficiency. Therefore, it is urgent to develop feasible shallow stratum exploration technologies for deep-sea rare earths. At present, there is no systematic acquisition - processing - interpretation process for using shallow stratum exploration technologies in rare earth resource investigations in China, and the application effect is poor, still in the exploratory stage.

[0058] In view of the problem of how to efficiently analyze the top and bottom interfaces of deep-sea rare earth enrichment horizons, the inventors of the present application proposed a processing solution to improve the imaging accuracy of deep-sea rare earth shallow stratum profile data. The solution of the present application will be carried out from two aspects: fine processing of deep-sea rare earth shallow stratum profiles and attribute analysis. The fine processing will mainly focus on denoising and improving the signal-to-noise ratio to improve the fidelity; the attribute analysis is intended to deeply mine the attribute information, characteristic information, etc. closely related to deep-sea rare earth resources contained in the deep-sea acoustic information on the basis of the fine processing of shallow stratum profile data, and realize the picking of the top and bottom interfaces of deep-sea rare earth enrichment horizons and the spatial distribution information.

[0059] The embodiments of the present application provide a method, device, equipment, and storage medium for imaging deep-sea rare earth shallow stratum profiles. The method includes: preprocessing the initial seismic data of deep-sea rare earth shallow stratum profiles; correcting and removing noise from the preprocessed initial seismic data to obtain refined seismic data; and determining the thickness distribution image of deep-sea rare earth shallow stratum profiles according to the refined seismic data. By preprocessing, correcting, and denoising the initial seismic data, and then using the refined seismic data to determine the distribution image of deep-sea rare earth shallow stratum profiles, the present application greatly improves the exploration accuracy and exploration efficiency of deep-sea rare earth shallow stratum profiles.

[0060] The embodiments of the present application provide a deep-sea rare earth shallow subsurface profile imaging method, device, equipment and storage medium, which relate to the technical field of data processing. The deep-sea rare earth shallow subsurface profile imaging method, device, equipment and storage medium provided by the embodiments of the present application can be applied to a terminal, can also be applied to a server, or can be software running on a terminal or a server. In some embodiments, the terminal can be a smart phone, a tablet computer, a laptop computer, a desktop computer, a smart speaker, a smart watch, a vehicle-mounted terminal, etc., but is not limited thereto; the server side can be configured as an independent physical server, can also be configured as a server cluster or a distributed system composed of multiple physical servers, or can be configured as a cloud server providing basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN, and big data and artificial intelligence platforms. The server can also be a node server in a blockchain network; the software can be an application implementing the deep-sea rare earth shallow subsurface profile imaging method, etc., but is not limited to the above forms.

[0061] The present application can be used in many general or special computer system environments or configurations. For example: personal computers, server computers, handheld or portable devices, tablet devices, multi-processor systems, microprocessor-based systems, set-top boxes, programmable consumer electronic devices, network PCs, minicomputers, mainframe computers, distributed computing environments including any of the above systems or devices, and so on. The present application can be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types. The present application can also be practiced in a distributed computing environment where tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules can be located in local and remote computer storage media including storage devices.

[0062] Referring to Figure 1 , the embodiments of the present application provide a deep-sea rare earth shallow subsurface profile imaging method, which may include but is not limited to S100 to S120, as follows:

[0063] S100: Preprocess the initial seismic data of the deep-sea rare earth shallow subsurface profile.

[0064] Further, S100 may include the following steps S101 to S104:

[0065] S101: Determine the sound velocity of the initial seismic data, the excitation frequency, and the number of segmented record files of each segmented survey line of the deep-sea rare earth shallow subsurface profile;

[0066] S102: Reset the sequence number of each segmented survey line according to the acquisition record time of each segmented record file;

[0067] S103: Unify the lengths of the segmented record files of each segmented survey line;

[0068] S104: Merge the segmented record files of each segmented survey line according to the sequence number of each segmented survey line to obtain the preprocessed initial seismic data.

[0069] S110: Correct and remove noise from the preprocessed initial seismic data to obtain refined seismic data.

[0070] As an optional implementation manner, the step of correcting the preprocessed initial seismic data in S110 includes the following steps S111 - S112:

[0071] S111: Perform delay correction on the preprocessed initial seismic data so that the corresponding seabed information in the preprocessed initial seismic data is continuous, and the starting value of the file number of the segmented record file in the preprocessed initial seismic data corresponds one-to-one with the reset sequence number of each segmented survey line;

[0072] S112: Perform amplitude correction on the preprocessed initial seismic data.

[0073] As an optional implementation manner, the step of removing noise from the preprocessed initial seismic data in S110 includes the following steps S113 - S114:

[0074] S113: Remove the data exceeding the preset amplitude threshold from the seismic traces of the preprocessed initial seismic data;

[0075] S114: Remove the data outside the target wave frequency band from the preprocessed initial seismic data.

[0076] S120: Determine the thickness distribution image of the deep-sea rare earth shallow stratum profile according to the refined seismic data.

[0077] Further, S120 may include the following steps S121 - S123:

[0078] S121: Extract the amplitude attribute from the refined seismic data;

[0079] S122: Determine the top layer and bottom boundary of the deep-sea rare earth shallow stratum profile according to the amplitude attribute;

[0080] S123: Determine the thickness distribution image of the deep-sea rare earth shallow stratum profile according to the top layer and bottom boundary.

[0081] More specifically, S123 may include the following steps:

[0082] Perform time-depth conversion based on the time curves of the top layer and the bottom boundary, the deep-sea drilling plan data of the deep-sea rare earth shallow stratum profile, and the empirical velocity around the deep-sea rare earth shallow stratum profile, so as to obtain the vertical thickness distribution image of the deep-sea rare earth shallow stratum profile.

[0083] As a further implementation manner, the embodiments of the present application may further include the following steps:

[0084] S130: Convert the geographical angular coordinates in the refined seismic data into longitude and latitude coordinates as the geographical coordinates of the deep-sea rare earth shallow stratum profile.

[0085] Next, specific application examples will be combined to introduce and explain the solutions of the embodiments of the present application in detail.

[0086] The solution of this embodiment will be carried out from two aspects: fine processing of the deep-sea rare earth shallow stratum profile and multi-attribute extraction and inversion. The fine processing will mainly focus on denoising and improving the signal-to-noise ratio to improve the fidelity. The attribute analysis is planned to carry out amplitude attribute extraction based on the fine processing of the shallow stratum profile data, so as to obtain the top and bottom boundary layer data of the high-quality rare earth-rich sedimentary layer.

[0087] 1. Taking the two-dimensional survey line of the deep-sea rare earth shallow stratum profile as an example, a processing solution for improving the imaging accuracy of the deep-sea rare earth shallow stratum profile data is proposed, making full use of the two major technologies of fine processing of the shallow stratum profile and multi-attribute extraction and inversion: effectively combining the two and applying them to the imaging of the deep-sea rare earth shallow stratum profile data, so as to better improve the imaging accuracy of the deep-sea rare earth shallow stratum profile data and explore the top and bottom interfaces of the deep-sea rare earth enrichment layer.

[0088] 2. The fine processing technology of the deep-sea rare earth shallow stratum profile, including resetting the serial numbers of segmented survey lines, loading segmented survey line data, resetting the recording lengths of segmented survey lines, merging and quality control of segmented survey line data, amplitude correction, burst noise elimination, frequency-wavenumber domain filtering, coordinate conversion, delay correction, and SEGY format output.

[0089] 3. The multi-attribute extraction and inversion technology. On the shallow stratum profile data after fine processing, deeply excavate the attribute information, characteristic information, etc. closely related to mineral resources contained in the deep-sea acoustic information to realize the picking of the top and bottom interfaces of the deep-sea rare earth enrichment layer and the spatial distribution information. Then, perform time-thickness conversion according to the time curves of the top and bottom interfaces of the obtained deep-sea rare earth enrichment layer in combination with the ocean drilling data and empirical velocity of the study area, and finally determine the vertical thickness distribution of the deep-sea rare earth-rich sedimentary layer.

[0090] The purpose of this embodiment is to propose a processing solution for improving the imaging accuracy of deep-sea rare earth shallow stratum profile data, aiming to deeply explore the attribute information, characteristic information, etc. closely related to mineral resources contained in deep-sea acoustic information, realize the picking up of the top and bottom interfaces of the deep-sea rare earth enrichment horizons and their spatial distribution information, so as to better and efficiently analyze the problems of the top and bottom interfaces of the deep-sea rare earth enrichment horizons, and track the vertical thickness and spatial distribution of the deep-sea rare earth-rich sedimentary layer.

[0091] Referring to Figure 2 , the fine processing technology of the deep-sea rare earth shallow stratum profile in this embodiment, and the fine processing includes (1) resetting the serial numbers of segmented survey lines, (2) loading the data of segmented survey lines, (3) resetting the record lengths of segmented survey lines, (4) merging and quality control of the data of segmented survey lines, (5) amplitude correction, (6) eliminating burst noise, (7) frequency-wavenumber domain filtering, (8) coordinate transformation, and (9) SEGY format output.

[0092] The establishment of the multi-attribute extraction and inversion technology includes: (1) the layer boundaries of multi-attribute extraction; (2) multi-attribute inversion; and (3) the vertical thickness of the deep-sea rare earth-rich sedimentary layer.

[0093] Specifically, a processing solution for improving the imaging accuracy of deep-sea rare earth shallow stratum profile data is as follows:

[0094] Step 1: Taking the deep-sea rare earth shallow stratum profile survey line as an example, conduct a tabular report collation on the original shallow stratum profile data collected in the field at sea, and determine the sound velocity, excitation frequency, and the number of segmented record files for each segmented survey line. The sound velocity is 1531 m / s, the excitation frequency is 6 kHz, and the number of segmented files is 3.

[0095] Exemplarily, FIGS. 3(a), 3(b), and 3(c) are example diagrams of the segmented survey lines of an actual deep-sea rare earth shallow stratum profile collected in the field. This segmented survey line recorded three segmented file numbers during field collection, namely FIGS. 3(a), 3(b), and 3(c). In FIGS. 3(a), 3(b), and 3(c), the horizontal direction of the profile is the trace number, the vertical direction is the time (milliseconds), the effective recording range of the data is 0 - 500 ms, and there is a recording delay in the data.

[0096] Step 2: Reset the serial numbers of the segmented survey lines. Since the serial numbers of the data record files of each segmented survey line start from 1, it is easy to get confused during the data merging process. Reset the serial numbers of each segmented survey line in sequence according to the acquisition record time of the record file names.

[0097] Step 3: Loading of segmented survey line data and resetting of record lengths. When loading segmented survey line data, ensure the accuracy of relevant position information at the trace headers (X and Y coordinates, recorded delay trace headers, shot numbers, file numbers). During field acquisition, due to recorded delays, the record lengths of data from different segments of the same survey line may vary. To ensure normal subsequent data merging, unify the record lengths of the loaded segmented survey line data.

[0098] Step 4: Merging and quality control of segmented survey line data. After Step 3, merge the segmented survey line data and simultaneously perform quality control on the merged data. This can be used for quality control of delay correction. After delay correction, the seabed of the survey line data should be continuous, and the starting value of the file number after merged recording should correspond one by one to the sequence number after resetting of the segmented survey line.

[0099] Exemplarily, Figure 4 is the original shallow stratigraphic profile after data merging obtained from a segmented survey line through Steps 2, 3, and 4. Figure 4 In the middle profile, Source is the trace number, Pick2 is the recorded delay trace header, the vertical axis is time, and the effective recording range of the data is 0 - 300 ms (milliseconds).

[0100] Exemplarily, Figure 5 is a shallow stratigraphic profile after delay correction. In this profile, Source is the range of trace numbers, Pick2 is the recorded delay trace header, the vertical axis is time, and the data recording time range after delay correction is 6700 - 7200 ms (milliseconds).

[0101] Step 5: Amplitude correction. Perform different types of abnormal amplitude corrections on the data.

[0102] Exemplarily, Figures 6(a) and 6(b) are the shallow stratigraphic profiles before and after amplitude correction respectively. Perform different types of abnormal amplitude corrections on the merged shallow stratigraphic profile data in Step 4. It can be seen from the comparison between Figures 6(a) and 6(b) that after amplitude correction, some types of abnormal amplitudes are suppressed, but there are still burst noises and other noises.

[0103] Step 6: Elimination of burst noises. Used to remove high - amplitude burst noises from seismic traces.

[0104] Exemplarily, Figure 7 is the shallow stratigraphic profile after burst noise suppression. After Step 5, there are still individual high - amplitude burst noises. After burst noise suppression, the high - amplitude burst noises in the seismic traces are effectively suppressed. Compared with Figure 6(b), Figure 7 the overall signal - to - noise ratio of the profile is improved.

[0105] Step 7: Frequency-wavenumber domain filtering. Frequency filtering has an obvious effect on eliminating interference outside the effective wave frequency band. After filtering, the background random noise is effectively suppressed.

[0106] Exemplarily, Figure 8 is the shallow strata profile after band-pass filtering. Frequency filtering has an obvious effect on eliminating interference outside the effective wave frequency band. After Step 6, the background random noise is effectively suppressed, and the overall signal-to-noise ratio of the profile is improved.

[0107] Step 8: Coordinate transformation. The geographic coordinates are recorded in the field data, and the unit is angular unit (ms). Theoretically, the interpretation system (Geoframe) can directly load the geographic coordinates in degrees when loading seismic data. Therefore, when loading seismic data, only the conversion of ms to degrees is required, and then the seismic data and coordinate information can be loaded simultaneously.

[0108] Step 9: Extraction of attributes and characteristic information. After Steps 2 to 8, the finely processed shallow strata profile is obtained. The instantaneous amplitude attribute is extracted from the finely processed shallow strata profile, and the top and bottom interfaces of the deep-sea rare earth enrichment horizon (i.e., the transparent layer) are picked on the amplitude attribute profile, so as to obtain the top and bottom interfaces of the deep-sea rare earth-rich sedimentary layer.

[0109] Exemplarily, Figures 9(a) and 9(b) are the shallow strata profiles before and after fine processing. After Steps 1 to 7, various noises in the profile are effectively suppressed, and the overall signal-to-noise ratio of the shallow strata profile data is greatly improved.

[0110] Exemplarily, Figure 10 is the instantaneous amplitude attribute of the shallow strata profile data. The instantaneous amplitude attribute reflects the magnitude of the reflection coefficient of the underground reflection interface. Generally, a high instantaneous amplitude is related to a strong reflection interface underground, which reflects a strong change in the lithology of adjacent underground rock layers or the presence of a gas-bearing layer. The gradual lateral change of the instantaneous amplitude usually reflects the change in the thickness of the rock layer or the lateral change in lithofacies. From Figure 10 the instantaneous amplitude attribute, it can be seen that the thickness of the deep-sea rare earth-rich muddy layer (transparent layer) from the western part to the eastern part of the survey line generally shows a trend of thick-thin-thick. The western part of the survey line is a homogeneous thin muddy layer, and the main feature below the seabed is transparent reflection. Near and at the seamount, it shows non-penetrating and strong reflectivity characteristics, corresponding to the outcrop of the bedrock without soft sediment cover.

[0111] Exemplarily, Figure 11 is a schematic diagram of the top and bottom interfaces of the deep-sea rare earth-rich sedimentary layer in the shallow strata profile data.

[0112] Using the characteristics of shallow stratigraphic profiles to distinguish pelagic sediment types is of great significance in the investigation of deep-sea rare earth-rich sediments. The results of recent investigations and studies have shown that rare earth elements are mainly hosted in pelagic clay and zeolite clay sediment layers in the Pacific Ocean, and the REY content therein is generally higher than 400x10 -6 . Generally, calcareous sediments deposited earlier are below the opaque chert layer. In the shallow stratigraphic profile diagram, several meters below the seabed is a continuously distributed acoustically transparent layer, corresponding to a pelagic clay layer about 2 m thick at the top of the sediment core; below it is a dark acoustically opaque layer, corresponding to a chert mixed layer with a high acoustic impedance. Accordingly, the transparent layer can be regarded as the target hosting horizon of deep-sea rare earth sediments; the siliceous and calcareous layered layers are horizons that are not conducive to the enrichment of rare earth elements; the opaque layer can be regarded as the lower boundary of the rare earth-rich sediment layer.

[0113] Based on the fact that the transparent layer is regarded as the target hosting horizon of deep-sea rare earth sediments and the opaque layer is regarded as the bottom boundary of the rare earth-rich sediment layer, combined with the instantaneous amplitude attribute of the shallow stratigraphic profile data, the top and bottom interfaces of the shallow stratigraphic profile data are determined by picking.

[0114] Step Ten: For the vertical thickness of the deep-sea rare earth-rich sediment layer, based on the top and bottom time curves of the top and bottom boundaries of the extracted deep-sea rare earth-rich sediment layer, combined with DSDP data and the empirical velocity around the study area, time-depth conversion is carried out, and then the vertical thickness distribution map of the deep-sea rare earth-rich sediments is obtained.

[0115] Exemplarily, Figure 12 is the thickness distribution map of the deep-sea rare earth-rich sediment layer in a single shallow profile line.

[0116] The layer values of the top and bottom interfaces of the deep-sea rare earth-rich sediment layer are respectively output as text. By performing subtraction on these two values of the top and bottom interfaces of the deep-sea rare earth-rich sediment layer, the time difference between the top and bottom of the deep-sea rare earth-rich sediment layer can be obtained. Combined with DSDP data and the empirical velocity around the study area, based on the seismic velocity of 1.74 km / s for calculation, after simple operations, the deposition thickness of the deep-sea rare earth-rich sediment layer can be calculated as Figure 10 shown. From Figure 10 it can be seen that the thickness range of the deep-sea rare earth-rich sediment layer is roughly 2 m - 20 m, with a minimum of 2 m and a maximum of 20 m.

[0117] Exemplarily, Figure 13 is the thickness plane distribution map of the deep-sea rare earth-rich sediment layer for multiple shallow profile lines.

[0118] Repeat the above steps and apply the above method to other shallow profiling lines in the study area. In the entire study area, considering the transparent layer as the target occurrence horizon of deep-sea rare-earth-rich sediments and the opaque layer as the bottom boundary of the rare-earth-rich sediment layer, the thickness range and planar distribution characteristics of the rare-earth-rich sediment layer in the study area can be identified. It can be seen that the thickness range of the rare-earth-rich sediment layer in the study area is approximately 6m - 36m, with a minimum of 6m and a maximum of 36m. Generally, it shows an east-west zoning feature, being thinner on both the east and west sides, thicker in the west than in the east, and thickest in the middle, with its thickness mainly ranging between 12m - 36m and the maximum thickness reaching 36m. Among them, the thickness of the relatively thick belt in the middle mainly ranges between 30m - 36m, which will be the key area for future rare-earth resource exploration.

[0119] Referring to Figure 14 , the embodiment of the present application also provides a deep-sea rare-earth shallow-stratum profile imaging device, which can implement the above-mentioned deep-sea rare-earth shallow-stratum profile imaging method. The device includes:

[0120] A data preprocessing unit for preprocessing the initial seismic data of the deep-sea rare-earth shallow-stratum profile;

[0121] A data refinement unit for correcting and removing noise from the preprocessed initial seismic data to obtain refined seismic data;

[0122] A thickness distribution determination unit for determining the thickness distribution image of the deep-sea rare-earth shallow-stratum profile based on the refined seismic data.

[0123] It can be understood that the content in the above method embodiments is applicable to the device embodiments of the present application. The functions specifically implemented by the device embodiments of the present application are the same as those in the above method embodiments, and the beneficial effects achieved are also the same as those in the above method embodiments.

[0124] The embodiment of the present application also provides an electronic device. The electronic device includes a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, it implements the above-mentioned deep-sea rare-earth shallow-stratum profile imaging method. The electronic device can be any intelligent terminal including a tablet computer, in-vehicle computer, etc.

[0125] It can be understood that the content in the above method embodiments is applicable to the device embodiments of the present application. The functions specifically implemented by the device embodiments of the present application are the same as those in the above method embodiments, and the beneficial effects achieved are also the same as those in the above method embodiments.

[0126] Please refer to Figure 15 , Figure 15 which shows the hardware structure of an electronic device in another embodiment. The electronic device includes:

[0127] The processor 1501 can be implemented in the form of a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, etc., and is used to execute relevant programs to implement the technical solutions provided in the embodiments of the present application;

[0128] The memory 1502 can be implemented in the form of a read-only memory (ROM), a static storage device, a dynamic storage device, or a random access memory (RAM), etc. The memory 1502 can store an operating system and other application programs. When implementing the technical solutions provided in the embodiments of this specification through software or firmware, the relevant program codes are stored in the memory 1502 and are called by the processor 1501 to execute the deep-sea rare earth shallow subsurface profiling method of the embodiments of the present application;

[0129] The input / output interface 1503 is used to implement information input and output;

[0130] The communication interface 1504 is used to implement communication interaction between this device and other devices, and can achieve communication through a wired method (such as USB, network cable, etc.) or through a wireless method (such as mobile network, WIFI, Bluetooth, etc.);

[0131] The bus 1505 transmits information between the various components of the device (such as the processor 1501, the memory 1502, the input / output interface 1503, and the communication interface 1504);

[0132] Among them, the processor 1501, the memory 1502, the input / output interface 1503, and the communication interface 1504 achieve communication connections with each other inside the device through the bus 1505.

[0133] The embodiments of the present application also provide a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the above-mentioned deep-sea rare earth shallow subsurface profiling method is implemented.

[0134] It can be understood that the content in the above method embodiments is applicable to the embodiments of this storage medium. The functions specifically implemented by the embodiments of this storage medium are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those of the above method embodiments.

[0135] The memory, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs and non-transitory computer-executable programs. In addition, the memory may include high-speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, a flash memory device, or other non-transitory solid-state storage devices. In some embodiments, the memory may optionally include memories remotely disposed relative to the processor, and these remote memories can be connected to the processor through a network. Examples of the above networks include but are not limited to the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0136] The embodiments described in the embodiments of the present application are for more clearly illustrating the technical solutions of the embodiments of the present application, and do not constitute a limitation on the technical solutions provided by the embodiments of the present application. Those skilled in the art will know that with the evolution of technology and the emergence of new application scenarios, the technical solutions provided by the embodiments of the present application are equally applicable to similar technical problems.

[0137] Those skilled in the art can understand that the technical solutions shown in the figures do not constitute a limitation on the embodiments of the present application, and may include more or fewer steps than those shown in the figures, or combine certain steps, or different steps.

[0138] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated, that is, they may be located in one place, or may be distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0139] Those of ordinary skill in the art can understand that all or some of the steps in the methods disclosed above, and the functional modules / units in the systems and devices, can be implemented as software, firmware, hardware, and appropriate combinations thereof.

[0140] The terms "first", "second", "third", "fourth", etc. (if any) in the specification of the present application and the above-mentioned drawings are used to distinguish similar objects, and do not have to be used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged under appropriate circumstances, so that the embodiments of the present application described here can be implemented in an order other than those illustrated or described here. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or units does not have to be limited to those steps or units clearly listed, but may include other steps or units not clearly listed or inherent to these processes, methods, products, or devices.

[0141] It should be understood that in this application, "at least one (item)" means one or more, and "a plurality" means two or more. "And / or" is used to describe the association relationship of associated objects, indicating that there can be three relationships. For example, "A and / or B" can mean: only A exists, only B exists, and both A and B exist at the same time. Among them, A and B can be singular or plural. The character " / " generally represents an "or" relationship between the associated objects before and after. "At least one (item) of the following" or its similar expressions refer to any combination of these items, including any combination of single items (items) or plural items (items). For example, at least one (item) of a, b, or c can mean: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", where a, b, and c can be single or multiple.

[0142] In several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the above division of units is only a logical function division. In actual implementation, there can be other division methods. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed coupling or direct coupling or communication connection between each other can be through some interfaces. The indirect coupling or communication connection of devices or units can be in electrical, mechanical or other forms.

[0143] The units described above as separate components may or may not be physically separated. The components displayed as units may or may not be physical units, that is, they can be located in one place, or they can be distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0144] In addition, each functional unit in various embodiments of this application can be integrated in a processing unit, or each unit can exist physically alone, or two or more units can be integrated in one unit. The above integrated units can be implemented in the form of hardware or in the form of software functional units.

[0145] When the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on such understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or all or part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes multiple instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods of various embodiments of this application. The aforementioned storage medium includes: various media that can store programs such as USB flash drives, mobile hard disks, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical discs.

[0146] The preferred embodiments of the embodiments of this application have been described above with reference to the accompanying drawings, which does not limit the scope of rights of the embodiments of this application. Any modifications, equivalent replacements, and improvements made by those skilled in the art without departing from the scope and essence of the embodiments of this application shall be within the scope of rights of the embodiments of this application.

Claims

1. Deep-sea rare earth shallow stratum profile imaging method, characterized in that, The method includes the following steps: Preprocess the initial seismic data of the shallow stratum profile of deep-sea rare earth. Correct and remove noise from the preprocessed initial seismic data to obtain refined seismic data. Determine the thickness distribution image of the shallow stratum profile of deep-sea rare earth according to the refined seismic data.

2. The deep-sea rare earth shallow subsurface profiling imaging method according to claim 1, wherein The preprocessing of the initial seismic data of the shallow stratum profile of deep-sea rare earth includes the following steps: Determine the sound velocity of the initial seismic data, the excitation frequency, and the number of segmented record files of each segmented survey line of the shallow stratum profile of deep-sea rare earth. Reset the serial number of each segmented survey line according to the acquisition record time of each segmented record file. Unify the lengths of each segmented record file of each segmented survey line. Merge each segmented record file of each segmented survey line according to the serial number of each segmented survey line to obtain the preprocessed initial seismic data.

3. The deep-sea rare earth shallow subsurface profiling imaging method according to claim 1, wherein The steps for correcting the preprocessed initial seismic data include the following steps: Perform delay correction on the preprocessed initial seismic data so that the corresponding seabed information in the preprocessed initial seismic data is continuous, and the starting value of the file number of the segmented record file in the preprocessed initial seismic data corresponds one by one to the reset serial number of each segmented survey line. Perform amplitude correction on the preprocessed initial seismic data.

4. The deep-sea rare earth shallow subsurface profiling imaging method according to claim 1, wherein The steps for removing noise from the preprocessed initial seismic data include the following steps: Remove the data exceeding the preset amplitude threshold from the seismic traces of the preprocessed initial seismic data. Remove the data outside the target wave frequency band in the preprocessed initial seismic data.

5. The deep-sea rare earth shallow subsurface profile imaging method according to claim 1, characterized in that The determination of the thickness distribution image of the shallow stratum profile of deep-sea rare earth according to the refined seismic data includes the following steps: Extract the amplitude attribute from the refined seismic data. Determine the top layer and bottom boundary of the shallow stratum profile of deep-sea rare earth according to the amplitude attribute. Determine the thickness distribution image of the shallow stratum profile of deep-sea rare earth according to the top layer and bottom boundary.

6. The deep-sea rare earth shallow subsurface profile imaging method according to claim 5, characterized in that, The determination of the thickness distribution image of the shallow stratum profile of deep-sea rare earth according to the top layer and bottom boundary includes the following steps: Perform time-depth conversion according to the time curves of the top layer and bottom boundary, the deep-sea drilling plan data of the shallow stratum profile of deep-sea rare earth, and the empirical velocity around the shallow stratum profile of deep-sea rare earth, so as to obtain the vertical thickness distribution image of the shallow stratum profile of deep-sea rare earth.

7. The deep-sea rare earth shallow subsurface profile imaging method according to any one of claims 1 to 6, characterized in that, The method further includes the following steps: Convert the geographical angular coordinates in the refined seismic data into longitude and latitude coordinates as the geographical coordinates of the shallow stratum profile of deep-sea rare earth.

8. Deep-sea rare earth shallow stratum profile imaging device, characterized in that, The device includes: A data preprocessing unit for preprocessing the initial seismic data of the shallow stratum profile of deep-sea rare earth. A data refinement unit for correcting and removing noise from the preprocessed initial seismic data to obtain refined seismic data. A thickness distribution determination unit for determining the thickness distribution image of the shallow stratum profile of deep-sea rare earth according to the refined seismic data.

9. An electronic device, characterized in that, The electronic device includes a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, it implements the deep-sea rare earth shallow subsurface profiling method according to any one of claims 1 to 7.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the deep-sea rare earth shallow subsurface profiling method according to any one of claims 1 to 7.

Images