Deep-sea rare earth sub-bottom profile imaging method, device, equipment and storage medium

Abstract

The application discloses a deep-sea rare earth subsoil profile imaging method, device and equipment and a storage medium, relates to the technical field of data processing, and comprises the following steps: preprocessing initial seismic data of a deep-sea rare earth subsoil profile; correcting and removing noise from the initial seismic data after preprocessing to obtain refined seismic data; and determining a thickness distribution image of the deep-sea rare earth subsoil profile according to the refined seismic data. The application preprocesses, corrects and denoises the initial seismic data to obtain refined seismic data, and then determines the distribution image of the deep-sea rare earth subsoil profile by using the refined seismic data, so that the exploration accuracy and efficiency of the deep-sea rare earth subsoil profile are greatly improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of data processing, and in particular to a deep-sea rare earth sub-bottom profile imaging method, device, equipment and storage medium. BACKGROUND

[0002] Deep-sea rare earth is a new type of marine mineral resource rich in heavy rare earth elements, and its resource reserves far exceed those on land, with great potential. At present, the method of combining sediment core sampling and sub-bottom detection is mainly used for deep-sea rare earth resource investigation. According to relevant research, the high-grade rare earth enrichment layer in seabed sediments can be as deep as several hundred meters below the surface. However, the commonly used box, grab and gravity piston column sampling techniques cannot meet the requirements of deep-sea sediment sampling over 10 meters in depth, and the obtained cores do not penetrate the bottom boundary of the enrichment layer. Deep-sea drilling requires a professional drilling ship, which is costly and time-consuming. SUMMARY

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

[0004] To achieve the above-mentioned purpose, one aspect of the embodiments of the present application proposes a deep-sea rare earth sub-bottom profile imaging method, which comprises the following steps:

[0005] Pretreating the initial seismic data of the deep-sea rare earth sub-bottom profile;

[0006] Correcting and removing noise from the pretreated initial seismic data as refined seismic data;

[0007] Determining the thickness distribution image of the deep-sea rare earth sub-bottom profile according to the refined seismic data.

[0008] In some embodiments, the pretreatment of the initial seismic data of the deep-sea rare earth sub-bottom profile comprises the following steps:

[0009] Determining the sound speed, excitation frequency of the initial seismic data and the number of segment record files of each segment line of the deep-sea rare earth sub-bottom profile;

[0010] Resetting the serial number of each segment line according to the acquisition record time of each segment record file;

[0011] Uniforming the length of each segment record file of each segment line;

[0012] Merging each segment record file of each segment line according to the serial number of each segment line to obtain the pretreated initial seismic data.

[0013] In some embodiments, the step of correcting the pre-processed initial seismic data comprises the following steps:

[0014] The pre-processed initial seismic data is subjected to a time delay correction, so that the corresponding seabed information in the pre-processed initial seismic data is continuous, and the file number start value of the segmented record file in the pre-processed initial seismic data corresponds to the reset serial number of each segmented survey line.

[0015] The pre-processed initial seismic data is subjected to an amplitude correction.

[0016] In some embodiments, the step of removing noise from the pre-processed initial seismic data comprises the following steps:

[0017] Data exceeding a preset amplitude threshold is removed from the seismic trace of the pre-processed initial seismic data.

[0018] Data outside the target wave frequency band in the pre-processed initial seismic data is removed.

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

[0020] An amplitude attribute is extracted from the refined seismic data.

[0021] The top layer and bottom boundary of the deep-sea rare earth shallow stratigraphic profile are determined according to the amplitude attribute.

[0022] The thickness distribution image of the deep-sea rare earth shallow stratigraphic profile is determined according to the top layer and bottom boundary.

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

[0024] Time-depth conversion is performed according to the time curve of the top layer and bottom boundary, the deep-sea drilling plan data of the deep-sea rare earth shallow stratigraphic profile, and the empirical velocity around the deep-sea rare earth shallow stratigraphic profile, thereby obtaining the vertical thickness distribution image of the deep-sea rare earth shallow stratigraphic profile.

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

[0026] The geographic angular coordinates in the refined seismic data are converted into latitude and longitude coordinates as the geographic coordinates of the deep-sea rare earth shallow stratigraphic profile.

[0027] To achieve the above object, another aspect of the embodiments of the present application provides a deep-sea rare earth shallow stratigraphic profile imaging device, which comprises:

[0028] a data preprocessing unit configured to preprocess initial seismic data of a deep-sea rare earth shallow stratigraphic profile;

[0029] a data refinement unit configured to correct and remove noise from the preprocessed initial seismic data as refined seismic data;

[0030] a thickness distribution determination unit configured to determine a thickness distribution image of the deep-sea rare earth shallow stratigraphic profile according to the refined seismic data.

[0031] To achieve the above object, another aspect of the embodiment of the present application provides an electronic device, which comprises a memory and a processor, the memory stores a computer program, and the processor implements the deep-sea rare earth shallow stratigraphic profile imaging method when executing the computer program.

[0032] To achieve the above object, another aspect of the embodiment of the present application provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the deep-sea rare earth shallow stratigraphic profile imaging method.

[0033] The embodiment of the present application at least has the following beneficial effects:

[0034] The present application can preprocess initial seismic data of a deep-sea rare earth shallow stratigraphic profile, correct and remove noise from the preprocessed initial seismic data as refined seismic data, and determine a thickness distribution image of the deep-sea rare earth shallow stratigraphic profile according to the refined seismic data. The present application preprocesses, corrects and removes noise from the initial seismic data to obtain refined seismic data, and then determines a distribution image of the deep-sea rare earth shallow stratigraphic profile by using the refined seismic data, which greatly improves the exploration accuracy and efficiency of the deep-sea rare earth shallow stratigraphic profile. BRIEF DESCRIPTION OF DRAWINGS

[0035] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings needed to be used in the embodiment description. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can also be obtained by those skilled in the art without any creative effort based on these drawings.

[0036] Figure 1 A flowchart of the deep-sea rare earth shallow stratigraphic profile imaging method provided by the embodiment of the present application is shown in the figure.

[0037] Figure 2 An example flowchart of the deep-sea rare earth shallow stratigraphic profile imaging method provided by the embodiment of the present application is shown in the figure.

[0038] Fig. 3(a), Fig. 3(b), Fig. 3(c) are original shallow stratigraphic profile segment trace data diagrams provided by the embodiment of the present application;

[0039] Figure 4 Fig. 4 is a diagram of data after merging provided by the embodiment of the present application;

[0040] Figure 5 Fig. 5 is a diagram of shallow stratigraphic profile after fine processing and delay correction provided by the embodiment of the present application;

[0041] Fig. 6(a), Fig. 6(b) are diagrams of shallow stratigraphic profile before and after amplitude correction provided by the embodiment of the present application;

[0042] Figure 7 Fig. 7 is a diagram of shallow stratigraphic profile after burst noise suppression provided by the embodiment of the present application;

[0043] Figure 8 Fig. 8 is a diagram of shallow stratigraphic profile after band-pass filtering provided by the embodiment of the present application;

[0044] Fig. 9(a), Fig. 9(b) are diagrams of shallow stratigraphic profile before and after fine processing provided by the embodiment of the present application;

[0045] Figure 10 Fig. 10 is an example diagram of instantaneous amplitude attribute of shallow stratigraphic profile data provided by the embodiment of the present application;

[0046] Figure 11 Fig. 11 is a diagram of top and bottom interfaces of deep sea rare earth rich sediment layer in shallow stratigraphic profile data provided by the embodiment of the present application;

[0047] Figure 12 Fig. 12 is a diagram of thickness distribution of deep sea rare earth rich sediment layer in single shallow profile trace provided by the embodiment of the present application;

[0048] Figure 13 Fig. 13 is a diagram of thickness planar distribution of deep sea rare earth rich sediment layer in multiple shallow profile traces provided by the embodiment of the present application;

[0049] Figure 14 Fig. 14 is a structural schematic diagram of deep sea rare earth shallow stratigraphic profile imaging device provided by the embodiment of the present application;

[0050] Figure 15 Fig. 15 is a hardware structural schematic diagram of an electronic device provided by the embodiment of the present application. DETAILED DESCRIPTION

[0051] In order to make the purposes, technical solutions and advantages of the present application clearer, the present application will be further described in detail below with reference to the drawings and embodiments. It should be understood that the specific embodiments described herein are only intended to explain the present application and not to limit the present application. When the following description refers to the drawings, the same numbers in different drawings represent the same or similar elements unless otherwise indicated. The implementations described in the following exemplary embodiments do not represent all implementations consistent with embodiments of the present application. They are merely examples of apparatuses 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", and the like used in the present application can be used herein to describe various concepts, but unless specifically stated, 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 can also be referred to as the second information, and similarly, the second information can also be referred to as the first information. Depending on the context, the word "if" as used herein can be interpreted as "when" or "when" or "in response to determining".

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

[0054] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as understood by those skilled in the art 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 the embodiments of the present application are described in detail, first, some related technologies involved in the embodiments of the present application are described as follows:

[0056] In view of the problems of inaccuracy, high cost and long time-consuming in the exploration of deep-sea rare earth sub-bottom profile in the prior art, it is urgent to develop a fast and feasible deep-sea rare earth exploration scheme. The current sub-bottom exploration technology scheme (sub-bottom profile measurement) has been widely applied in the field of marine sub-bottom research. In the field of ocean mineral resource investigation, it is mainly used for studying the sediment distribution characteristics of deep-sea basins and seamounts, and is a conventional technology for detecting sub-bottom sediments on the seabed. However, the use of sub-bottom exploration technology in the field of rare earth resource investigation is still in the exploratory stage, and there is still no effective processing scheme for improving the imaging accuracy of deep-sea rare earth sub-bottom profile data.

[0057] Land-based rare earth ore is mainly light rare earth, and medium and heavy rare earth is still a key metal resource in short supply. In recent years, geologists have discovered sediments rich in medium and heavy rare earth in the deep sea, and the ore-forming characteristics are obviously different from land ion adsorption type rare earth ore, but there has been controversy about the occurrence state and super enrichment mechanism of deep sea rare earth, and there is also a lack of rapid exploration and evaluation technology for deep sea rare earth. At present, the method of combining sediment core sampling with shallow stratigraphic exploration is mainly used for deep sea rare earth resource investigation. According to the relevant technical research, the high-grade rare earth enrichment layer in the seabed sediment can be as deep as hundreds of meters below the surface. However, the commonly used box, grab and gravity piston column sampling techniques cannot meet the requirements of sampling deep sea rare earth rich sediment more than ten meters deep. The obtained cores basically do not penetrate the bottom boundary of the enrichment layer. Deep sea drilling requires a professional drilling ship, which is high in cost and low in exploration efficiency. Therefore, it is urgent to develop a feasible deep sea rare earth shallow stratigraphic exploration technology. At present, there is no systematic collection-processing-interpretation process for rare earth resource investigation using shallow stratigraphic exploration technology in China, and the application effect is poor, and it is still in the exploration stage.

[0058] The inventors of the present application propose a processing scheme for improving the imaging accuracy of deep sea rare earth shallow stratigraphic profile data in view of how to efficiently analyze the top and bottom interfaces of the deep sea rare earth enrichment layer. The scheme of the present application will be developed from two aspects of fine processing and attribute analysis of deep sea rare earth shallow stratigraphic profile. The fine processing will mainly improve the fidelity by denoising and improving the signal-to-noise ratio. The attribute analysis will further dig the attribute information and characteristic information closely related to deep sea rare earth resources from the deep sea acoustic information on the basis of fine processing of shallow stratigraphic profile data, so as to realize the picking and spatial distribution information of the top and bottom interfaces of the deep sea rare earth enrichment layer.

[0059] The embodiments of the present application provide a deep sea rare earth shallow stratigraphic profile imaging method, device, equipment and storage medium. The method comprises: preprocessing the initial seismic data of the deep sea rare earth shallow stratigraphic profile; correcting and removing noise from the preprocessed initial seismic data as fine seismic data; and determining the thickness distribution image of the deep sea rare earth shallow stratigraphic profile according to the fine seismic data. The initial seismic data is preprocessed, corrected and denoised to obtain fine seismic data, and then the distribution image of the deep sea rare earth shallow stratigraphic profile is determined by using the fine seismic data, which greatly improves the exploration accuracy and efficiency of the deep sea rare earth shallow stratigraphic profile.

[0060] The embodiment of the present application provides a deep-sea rare earth sub-bottom profile imaging method, device, equipment and storage medium, relates to the technical field of data processing. The deep-sea rare earth sub-bottom profile imaging method, device, equipment and storage medium provided by the embodiment of the present application can be applied to a terminal, can be applied to a server, and can also be software running in the terminal or the server. In some embodiments, the terminal can be a smart phone, a tablet computer, a notebook computer, a desktop computer, a smart speaker, a smart watch and a vehicle-mounted terminal, but is not limited thereto; the server end can be configured as a stand-alone physical server, can be configured as a server cluster or a distributed system formed by multiple physical servers, can be configured as a cloud server providing cloud service, cloud database, cloud computing, cloud function, cloud storage, network service, cloud communication, middleware service, domain name service, security service, CDN and basic cloud computing services such as big data and artificial intelligence platform, and the server can also be a node server in a blockchain network; and the software can be an application for implementing the deep-sea rare earth sub-bottom profile imaging method, 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, multiprocessor systems, microprocessor-based systems, set-top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, distributed computing environments including any of the above systems or devices, etc. 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, in which tasks are performed by remote processing devices connected by a communication network. In a distributed computing environment, program modules can be located in local and remote computer storage media, including storage devices.

[0062] Reference Figure 1 The embodiment of the present application provides a deep-sea rare earth sub-bottom profile imaging method, which can include but is not limited to S100 to S120, and specifically as follows:

[0063] S100: preprocessing initial seismic data of a deep-sea rare earth sub-bottom profile.

[0064] Further, S100 can include the following steps S101-S104:

[0065] S101: determining the sound velocity, excitation frequency of the initial seismic data and the number of segmented record files of each segmented measuring line of the deep-sea rare earth sub-bottom profile;

[0066] S102: resetting the serial number of each segmental survey line according to the acquisition record time of each segmental record file;

[0067] S103: unifying the length of each segmental record file of each segmental survey line;

[0068] S104: merging each segmental record file of each segmental survey line according to the serial number of each segmental survey line to obtain the pretreated initial seismic data.

[0069] S110: correcting and removing noise from the pretreated initial seismic data as refined seismic data.

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

[0071] S111: performing delay correction on the pretreated initial seismic data so that the corresponding seabed information in the pretreated initial seismic data is continuous, and the file number starting value of the segmental record file in the pretreated initial seismic data corresponds to the reset serial number of each segmental survey line one by one;

[0072] S112: performing amplitude correction on the pretreated initial seismic data.

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

[0074] S113: removing data exceeding a preset amplitude threshold from the seismic trace of the pretreated initial seismic data;

[0075] S114: removing data outside the target wave frequency band in the pretreated initial seismic data.

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

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

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

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

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

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

[0082] According to the time curve of the top layer and the bottom boundary, the deep-sea drilling plan data of the deep-sea rare earth shallow stratigraphic profile, and the empirical velocity around the deep-sea rare earth shallow stratigraphic profile, time-depth conversion is performed to obtain a vertical thickness distribution image of the deep-sea rare earth shallow stratigraphic profile.

[0083] As a further implementation, the embodiment of the application can further include the following steps:

[0084] S130: converting the geographic angle coordinates in the refined seismic data into latitude and longitude coordinates as geographic coordinates of the deep-sea rare earth shallow stratigraphic profile.

[0085] Next, the scheme of the embodiment of the application will be described and explained in detail in combination with specific application examples.

[0086] The embodiment scheme will be developed from two aspects of deep-sea rare earth shallow stratigraphic profile refinement processing and multi-attribute extraction and inversion. The refinement processing will mainly improve the fidelity by denoising and signal-to-noise ratio improvement; and the attribute analysis will be carried out on the basis of the refined shallow stratigraphic profile data to extract amplitude attributes and then obtain high-quality top and bottom boundary data of the rare earth-rich sedimentary layer.

[0087] 1. Taking a two-dimensional line of deep-sea rare earth shallow stratigraphic profile data as an example, a processing scheme for improving the imaging accuracy of deep-sea rare earth shallow stratigraphic profile data is proposed. The two technologies of shallow stratigraphic profile refinement processing and multi-attribute extraction and inversion are effectively combined and applied to the imaging of deep-sea rare earth shallow stratigraphic profile data, so as to better improve the imaging accuracy of deep-sea rare earth shallow stratigraphic profile data and prove the top and bottom boundaries of deep-sea rare earth-rich layers.

[0088] 2. Deep-sea rare earth shallow stratigraphic profile refinement processing technology. The refinement processing includes segmented line number resetting, segmented line data loading, segmented line record length resetting, segmented line data merging and quality control, amplitude correction, burst noise elimination, frequency-wave number domain filtering, coordinate conversion, delay correction, and SEGY format output.

[0089] 3. Multi-attribute extraction and inversion technology. On the basis of the refined shallow stratigraphic profile data, the attribute information and characteristic information closely related to mineral resources contained in deep-sea acoustic information are deeply mined to realize the picking and spatial distribution information of the top and bottom boundaries of deep-sea rare earth-rich layers. Then, according to the time curve of the obtained top and bottom boundaries of deep-sea rare earth-rich layers, combined with the ocean drilling data and empirical velocity of the study area, time-thickness conversion is performed to finally determine the vertical thickness distribution of deep-sea rare earth-rich sedimentary layers.

[0090] The purpose of the embodiment is to propose a processing scheme for improving the imaging accuracy of deep-sea rare earth sub-bottom profile data, the purpose being to deeply mine the attribute information and characteristic information related to mineral resources contained in deep-sea acoustic information, to realize the picking and spatial distribution information of the top and bottom interfaces of the deep-sea rare earth enrichment horizon, and to better and efficiently analyze the problem of the top and bottom interfaces of the deep-sea rare earth enrichment horizon, and to track the vertical thickness and spatial distribution of the deep-sea rare earth enrichment layer.

[0091] With reference to Figure 2 The deep-sea rare earth sub-bottom profile fine processing technology of the embodiment includes (1) segmented line number resetting, (2) segmented line data loading, (3) segmented line record length resetting, (4) segmented line data merging and quality control, (5) amplitude correction, (6) burst noise elimination, (7) frequency-wavenumber domain filtering, (8) coordinate conversion, and (9) SEGY format output.

[0092] The multi-attribute extraction and inversion technology includes (1) multi-attribute extraction layer interface, (2) multi-attribute inversion, and (3) vertical thickness of the deep-sea rare earth enrichment layer.

[0093] Specifically, a processing scheme for improving the imaging accuracy of deep-sea rare earth sub-bottom profile data has the following specific implementation steps:

[0094] Step one: Taking a deep-sea rare earth sub-bottom profile line as an example, the original sub-bottom profile data collected in the field is arranged, the sound velocity, the excitation frequency, and the number of segmented record files of each segmented line are determined. 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 a segmented line of a deep-sea rare earth sub-bottom profile actually collected in the field, which records three segmented files in the field collection, which are FIGS. 3(a), 3(b), and 3(c). In FIGS. 3(a), 3(b), and 3(c), the profile is transversely numbered, and longitudinally time (milliseconds) is recorded, the effective data recording range is 0-500 ms, and there is a recording delay in the data.

[0096] Step two: The number of segmented lines is reset. Since the number of data record files of each segmented line starts from 1, it is easy to be confused in the data merging process. According to the recording time of the record file name, the number of each segmented line is reset in turn.

[0097] Step three: segmented line data loading and record length resetting, when loading segmented line data, ensure that the relevant position information of the trace head is accurate (X, Y coordinates, record delay trace head, shot number, file number); during field acquisition, due to record delay, the record lengths of different segments of the same line are different, in order to ensure normal subsequent data merging, the record lengths of the loaded segmented line data are unified.

[0098] Step four: segmented line data merging and quality control, after step three, the segmented line data is merged, and the merged data is quality controlled, which can be delay corrected, after delay correction, the seabed of the line data should be continuous, and the starting value of the file number after merging the record should correspond to the sequence number after resetting the segmented line.

[0099] Exemplarily, Figure 4 is a data-merged original shallow seismic profile obtained after a segmented line passes through steps two, three, and four. Figure 4 The profile has a profile Source as a trace number, a Pick2 as a record delay trace head, a vertical direction as time, and a data effective record range of 0-300 ms (milliseconds).

[0100] Exemplarily, Figure 5 is a delay-corrected shallow seismic profile. The profile has a profile Source as a trace number range, a Pick2 as a record delay trace head, a vertical direction as time, and a delay-corrected data record time range of 6700-7200 ms (milliseconds).

[0101] Step five: amplitude correction, different types of abnormal amplitude correction are performed on the data.

[0102] Exemplarily, FIG. 6(a) and FIG. 6(b) are shallow seismic profiles before and after amplitude correction, respectively, different types of abnormal amplitude correction are performed on the merged shallow seismic profile data in step four. As can be seen from the comparison between FIG. 6(a) and FIG. 6(b), after amplitude correction, part of the types of abnormal amplitudes are suppressed, but there are still burst noises and other noises.

[0103] Step six: burst noise elimination, used to remove high-amplitude burst noises from seismic traces.

[0104] Exemplarily, Figure 7 is a shallow seismic profile after burst noise suppression, after step five, 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 FIG. 6(b), Figure 7 the overall signal-to-noise ratio of the profile is improved.

[0105] Step seven: frequency wave number domain filtering, frequency filtering has obvious effect on eliminating interference outside the effective wave frequency band, after filtering, the background random noise is effectively suppressed.

[0106] Exemplarily, Figure 8 The shallow stratum profile after band-pass filtering, frequency filtering has obvious effect on eliminating interference outside the effective wave frequency band, after step six, the background random noise is effectively suppressed, and the overall signal-to-noise ratio of the profile is improved.

[0107] Step eight: coordinate conversion, the geographic coordinates are recorded in the field data, and the unit is angle unit (ms). In theory, the interpretation system (Geoframe) can directly load the geographic coordinates in degrees when loading the seismic data, so when loading the seismic data, only the ms needs to be converted to degrees, that is, the seismic data and coordinate information can be loaded at the same time.

[0108] Step nine: attribute and feature information extraction, after steps two to eight, the fine-processed shallow stratum profile is obtained, the instantaneous amplitude attribute of the fine-processed shallow stratum profile is extracted, and the top and bottom interfaces of the deep-sea rare earth enrichment layer (i.e. transparent layer) are picked up on the amplitude attribute profile, so as to obtain the top and bottom interfaces of the deep-sea rare earth enrichment layer.

[0109] Exemplarily, FIG. 9(a) and FIG. 9(b) are shallow stratum profile graphs before and after fine processing, after steps one to seven, various noises in the profile are effectively suppressed, and the overall signal-to-noise ratio of the shallow stratum profile data is greatly improved.

[0110] Exemplarily, Figure 10 The instantaneous amplitude attribute of the shallow stratum profile data, the instantaneous amplitude attribute reflects the size of the reflection coefficient of the underground reflection interface. Generally, high instantaneous amplitude is related to strong reflection interface underground, which reflects the strong change of the lithology of adjacent rock layers underground, or the existence of gas-bearing layers. The gradual change of the instantaneous amplitude in the horizontal direction usually reflects the change of the thickness of the rock layer, or the horizontal change of the lithofacies. From Figure 10 From the instantaneous amplitude attribute of the profile, it can be seen that the thickness of the deep-sea rare earth argillaceous layer (transparent layer) in the west to the east of the survey line as a whole presents a thick-thin-thick trend, the west of the survey line is a homogeneous thin mud layer, and the transparent reflection is the main feature below the sea bottom, near and at the seamount, it shows the characteristics of non-penetration and strong reflection, which corresponds to the bedrock outcrop without soft sediment coverage.

[0111] Exemplarily, Figure 11 The top and bottom interfaces of the deep-sea rare earth enrichment layer in the shallow stratum profile data are shown in the schematic diagram.

[0112] It is of great significance to distinguish the types of pelagic sediments by using the characteristics of shallow seismic profiles. The results of recent investigations show that REY is mainly hosted in the Pacific Ocean pelagic clay and zeolite clay sediment layers, and the REY content is generally higher than 400x10 -6 Under the opaque chert layer, there are generally early deposited calcareous sediments. In the shallow seismic profile, the acoustic transparent layer is continuously distributed in the first several meters below the sea floor, corresponding to the pelagic clay layer with a thickness of about 2m at the top of the sediment core; and the dark acoustic opaque layer is below the transparent layer, corresponding to the high acoustic impedance chert mixed layer. Accordingly, the transparent layer can be regarded as the target hosting layer of deep-sea rare earth sediment; the siliceous and drill layers are the layers that are not conducive to the enrichment of rare earth elements; and the opaque layer can be regarded as the lower boundary of the rare earth-rich sediment layer.

[0113] According to the transparent layer regarded as the target hosting layer of deep-sea rare earth sediment and the opaque layer regarded as the bottom boundary of the rare earth-rich sediment layer, the top and bottom boundaries of the shallow seismic profile data are determined by picking up the instantaneous amplitude attribute of the shallow seismic profile data.

[0114] Step ten: The vertical thickness of the deep-sea rare earth-rich sediment layer is obtained by time-depth conversion based on the top and bottom time curves of the deep-sea rare earth-rich sediment layer, combined with the DSDP data and the empirical velocity of the surrounding area of the study area, and then the vertical thickness distribution map of the deep-sea rare earth-rich sediment 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 survey line.

[0116] The top and bottom boundary values of the deep-sea rare earth-rich sediment layer are output as text, and the time difference between the top and bottom boundaries of the deep-sea rare earth-rich sediment layer is obtained by subtracting the two values. Based on the seismic velocity calculation of 1.74km / s, the deposition thickness of the deep-sea rare earth-rich sediment layer can be calculated by simple operation, as shown in Figure 10 It can be seen from Figure 10 that the thickness of the deep-sea rare earth-rich sediment layer is generally in the range of 2m-20m, with a thickness of 2m in the thinnest place and a thickness of 20m in the thickest place.

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

[0118] The above steps are repeated, and the above method is applied to other shallow profile survey lines in the study area. In the entire study area, according to the transparent layer being regarded as a target occurrence layer of deep-sea rare earth-rich sediments and the opaque layer being regarded 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 are identified. It can be seen that the thickness range of the rare earth-rich sediment layer in the study area is approximately 6m-36m, the thinnest is 6m, the thickest is 36m, and the overall performance shows east-west zonation characteristics, with the west being relatively thicker than the east, and the middle being thicker. The thickness of the middle part is mainly between 12m and 36m, and the maximum thickness can reach 36m. The zonal region with greater thickness in the middle part, mainly between 30m and 36m, will be the key area for future rare earth resource exploration.

[0119] Reference Figure 14 The embodiment of the present application also provides a deep-sea rare earth shallow stratigraphic profile imaging device, which can implement the deep-sea rare earth shallow stratigraphic profile imaging method described above. The device comprises:

[0120] A data preprocessing unit is configured to preprocess initial seismic data of a deep-sea rare earth shallow stratigraphic profile.

[0121] A data refinement unit is configured to correct and remove noise from the preprocessed initial seismic data as refined seismic data.

[0122] A thickness distribution determination unit is configured to determine a thickness distribution image of the deep-sea rare earth shallow stratigraphic profile according to the refined seismic data.

[0123] It can be understood that the contents in the above method embodiments are all applicable to the present device embodiment. The present device embodiment specifically implements the same functions as the above method embodiments, and achieves the same beneficial effects as the above method embodiments.

[0124] The embodiment of the present application also provides an electronic device, which comprises a memory and a processor. The memory stores a computer program, and the processor implements the deep-sea rare earth shallow stratigraphic profile imaging method described above when executing the computer program. The electronic device can be any intelligent terminal, such as a tablet computer or a vehicle-mounted computer.

[0125] It can be understood that the contents in the above method embodiments are all applicable to the present device embodiment. The present device embodiment specifically implements the same functions as the above method embodiments, and achieves the same beneficial effects as the above method embodiments.

[0126] Please refer to Figure 15 , Figure 15 The hardware structure of the electronic device of another embodiment is illustrated, which comprises:

[0127] The processor 1501 can be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an ASIC (Application Specific Integrated Circuit), or one or more integrated circuits, and is configured to execute related programs to implement the technical solutions provided by the embodiments of the present application.

[0128] The memory 1502 can be implemented by a ROM (Read Only Memory), a static storage device, a dynamic storage device, or a RAM (Random Access Memory), and the like. The memory 1502 can store an operating system and other application programs. When the technical solutions provided by the embodiments of the present application are implemented by software or firmware, the related program codes are stored in the memory 1502 and are called and executed by the processor 1501 to implement the deep-sea rare earth subsoil profile imaging method of the embodiments of the present application.

[0129] The input / output interface 1503 is configured to implement information input and output.

[0130] The communication interface 1504 is configured to implement the communication interaction between the device and other devices. The communication can be realized by a wired manner (for example, a USB, a network cable, or the like) or a wireless manner (for example, a mobile network, WIFI, Bluetooth, or the like).

[0131] The bus 1505 is configured to transmit information between various components (for example, the processor 1501, the memory 1502, the input / output interface 1503, and the communication interface 1504) of the device.

[0132] The processor 1501, the memory 1502, the input / output interface 1503, and the communication interface 1504 are connected to each other by the bus 1505 to realize the communication connection between the devices.

[0133] The embodiments of the present application further provide a computer readable storage medium, which stores a computer program. The computer program is executed by a processor to implement the deep-sea rare earth subsoil profile imaging method.

[0134] It can be understood that the contents of the above method embodiments are applicable to the present storage medium embodiments. The present storage medium embodiments specifically implement the functions of the above method embodiments, and achieve the same beneficial effects as the above method embodiments.

[0135] 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 can include a high-speed random access memory and can also include a non-transitory memory, such as at least one magnetic disk storage device, a flash memory device, or other non-transitory solid-state memory device. In some embodiments, the memory can optionally include a memory that is remotely disposed relative to the processor, and these remote memories can be connected to the processor through a network. Examples of the above network 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 used to more clearly illustrate 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 can know that, with the evolution of technology and the appearance of new application scenarios, the technical solutions provided by the embodiments of the present application are also 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 can include more or fewer steps than shown in the figures, or combine certain steps, or different steps.

[0138] The device embodiments described above are only schematic, and units described as separate components can or can not be physically separate, i.e., can be located in one place, or can be distributed on multiple network units. Part or all of the modules can be selected according to actual needs to achieve the purpose of the embodiments of the present application.

[0139] Those skilled in the art can understand that all or some of the steps in the above disclosed method, the functional modules / units in the system and the device can be implemented as software, firmware, hardware and their appropriate combinations.

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

[0141] It should be understood that, in the application, "at least one" refers to one or more, and "multiple" refers to two or more. "And / or" is used to describe the association relationship of the associated objects, which means that there can be three relationships, for example, "A and / or B" can represent three cases of only A, only B, and A and B existing at the same time, wherein A and B can be singular or plural. The character " / " generally represents an "or" relationship between the associated objects before and after it. "At least one of the following" or similar expressions means any combination of these items, including any combination of single or multiple items. For example, at least one of a, b or c can represent a, b, c, "a and b", "a and c", "b and c", or "a and b and c", wherein a, b, and c can be single or multiple.

[0142] In several embodiments provided in the 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 only illustrative, for example, the division of the above units is only a logical function division, and actual implementation can have another division manner, for example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the displayed or discussed units can be indirect coupling or communication connection through some interfaces, devices or units, which can be electrical, mechanical or other forms.

[0143] The units described above as separate components can or can not be physically separated, and the components shown as units can or can not be physical units, that is, they can be located in one place, or can be distributed on multiple network units. According to actual needs, part or all of the units can be selected to achieve the purpose of the embodiment scheme.

[0144] In addition, the functional units in each embodiment of the application can be integrated in one processing unit, or each unit can be physically present separately, or two or more units can be integrated in one unit. The integrated unit can be realized in the form of hardware or in the form of a software functional unit.

[0145] The integrated unit, if implemented in the form of a software function unit and sold or used as an independent product, can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application, essentially or in other words, the part that contributes to the prior art or the whole or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes multiple instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods of the various embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, and various program storage media.

[0146] The preferred embodiments of the embodiments of the present application are described above with reference to the accompanying drawings, and are not limited to the scope of the embodiments of the present 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 the present application shall be within the scope of the embodiments of the present application.

Claims

1. A method for deep-sea rare earth shallow stratigraphic profile imaging, characterized in that, The method comprises the following steps: Pretreatment of initial seismic data of deep-sea rare earth subsoil profile; Correction and noise removal of the pretreated initial seismic data as refined seismic data; Determination of thickness distribution image of the deep-sea rare earth subsoil profile according to the refined seismic data; The pretreatment of initial seismic data of deep-sea rare earth subsoil profile comprises the following steps: Determination of sound velocity, excitation frequency of the initial seismic data and number of segment record files of each segment line of the deep-sea rare earth subsoil profile; Resetting the serial number of each segment line according to the acquisition record time of each segment record file; Uniforming the length of each segment record file of each segment line; Merging each segment record file of each segment line according to the serial number of each segment line to obtain the pretreated initial seismic data; The determination of thickness distribution image of the deep-sea rare earth subsoil profile according to the refined seismic data comprises the following steps: Extracting amplitude attribute from the refined seismic data; Determination of top layer and bottom boundary of the deep-sea rare earth subsoil profile according to the amplitude attribute; wherein the top layer is a transparent layer; Determination of thickness distribution image of the deep-sea rare earth subsoil profile according to the top layer and bottom boundary; The determination of thickness distribution image of the deep-sea rare earth subsoil profile according to the top layer and bottom boundary comprises the following steps: Time-depth conversion according to the time curve of the top layer and bottom boundary, deep-sea drilling plan data of the deep-sea rare earth subsoil profile and empirical velocity around the deep-sea rare earth subsoil profile, and then obtaining the vertical thickness distribution image of the deep-sea rare earth subsoil profile.

2. The deep ocean rare earth subbottom profile imaging method of claim 1, wherein, The step of correction of the pretreated initial seismic data comprises the following steps: Delay correction of the pretreated initial seismic data, so that the corresponding seabed information in the pretreated initial seismic data is continuous, and the file number starting value of the segment record file in the pretreated initial seismic data corresponds to the reset serial number of each segment line one by one; Amplitude correction of the pretreated initial seismic data.

3. The deep ocean rare earth subbottom profile imaging method of claim 1, wherein, The step of noise removal of the pretreated initial seismic data comprises the following steps: Removing data exceeding the preset amplitude threshold from the seismic trace of the pretreated initial seismic data; Removing data outside the target wave frequency band in the pretreated initial seismic data.

4. The deep ocean rare earth subbottom profiling method according to any one of claims 1 to 3, characterized in that, The method further comprises the following steps: Converting the geographic angle coordinates in the refined seismic data into latitude and longitude coordinates as the geographic coordinates of the deep-sea rare earth subsoil profile.

5. A deep-sea rare earth sub-bottom profiler apparatus, characterized by, The device is applied to the deep-sea rare earth subsoil profile imaging method as claimed in claim 1, and the device comprises: A data pretreatment unit for pretreating initial seismic data of a deep-sea rare earth subsoil profile; A data refinement unit for correcting and removing noise of the pretreated initial seismic data as refined seismic data; A thickness distribution determination unit for determining a thickness distribution image of the deep-sea rare earth subsoil profile according to the refined seismic data.

6. An electronic device, comprising: The electronic device comprises a memory and a processor, the memory stores a computer program, and the processor implements the deep-sea rare earth sub-bottom profile imaging method according to any one of claims 1 to 4 when executing the computer program.

7. A computer-readable storage medium storing a computer program, wherein the computer program comprises the following steps of: receiving a request for a resource from a client; determining whether the client is authorized to access the resource; and if the client is authorized to access the resource, providing the resource to the client. The computer program is executed by the processor to implement the deep-sea rare earth sub-bottom profile imaging method according to any one of claims 1 to 4.

Images