CN118295031A - Coastal zone stratigraphic structure detection method, device, electronic equipment and storage medium - Google Patents

Abstract

The present application discloses a method, device, electronic device and storage medium for detecting stratigraphic structure in a coastal zone, and relates to the field of geological exploration technology. The method comprises: obtaining detection data collected by multiple seismometers, and obtaining position data of each seismometer; wherein each seismometer is sequentially and continuously arranged on the land side and the ocean side of the target coastal zone according to a preset spacing; performing cross-correlation calculation on the detection data collected by two adjacent seismometers to obtain the Green's function corresponding to the two adjacent seismometers; performing surface wave dispersion curve imaging according to each Green's function and each position data to obtain the corresponding surface wave dispersion curve energy map; and obtaining the stratigraphic structure profile of the target coastal zone by inverting the energy map of each surface wave dispersion curve. The present application can obtain the shear wave velocity structure information of the shallow strata of the target coastal zone by arranging multiple seismometers on both the land side and the ocean side of the target coastal zone, thereby realizing seamless sea-land joint measurement.

Description

Coastal zone stratigraphic structure detection method, device, electronic equipment and storage medium

Technical Field

The present application relates to the field of geological exploration technology, and in particular to a method, device, electronic equipment and storage medium for detecting stratigraphic structures in coastal zones.

Background technique

For a long time, obtaining the geological structure and stratigraphic information of the shallow strata on both sides of the coast usually requires two different types of geophysical exploration methods, land and sea. The existing shallow strata survey technology in the coastal zone has the following problems: 1) The existing land/sea seismic exploration technology cannot directly achieve seamless land-sea joint survey; 2) The two existing land/sea seismic exploration technologies have different abilities to reveal the stratigraphic structure. Due to the differences in benchmarks, accuracy and resolution, the data obtained by the two methods cannot be effectively spliced.

Summary of the invention

The main purpose of the embodiments of the present application is to propose a method, device, electronic equipment and storage medium for detecting the stratigraphic structure of the coastal zone, so as to carry out seamless sea-land joint measurement of the coastal zone.

To achieve the above-mentioned purpose, one aspect of an embodiment of the present application provides a method for detecting stratigraphic structure in a coastal zone, the method comprising:

Acquire detection data collected by a plurality of seismographs, and acquire position data of each of the seismographs; wherein each of the seismographs is sequentially and continuously deployed on the land side and the sea side of the target coastal zone at a preset interval;

Performing cross-correlation calculation on the detection data collected by two adjacent seismographs to obtain Green's functions corresponding to the two adjacent seismographs;

Perform surface wave dispersion curve imaging according to each of the Green's functions and each of the position data to obtain a corresponding surface wave dispersion curve energy diagram;

The stratigraphic structure profile of the target coastal zone is obtained by inverting the energy diagrams of the surface wave dispersion curves.

In some embodiments, before performing cross-correlation calculation on the detection data acquired by two adjacent seismographs to obtain Green's functions corresponding to the two adjacent seismographs, the method further includes:

Preprocessing the detection data to obtain preprocessed data;

The cross-correlation calculation is performed on the detection data collected by two adjacent seismographs to obtain Green's functions corresponding to the two adjacent seismographs, including:

A cross-correlation calculation is performed on the pre-processed data corresponding to two adjacent seismographs to obtain the Green's functions corresponding to the two adjacent seismographs.

In some embodiments, preprocessing the detection data to obtain preprocessed data includes:

The detection data is subjected to at least one of instrument response removal, mean removal, trend removal, filtering, noise signal segmentation, time domain normalization, and spectrum whitening to obtain the preprocessed data.

In some embodiments, the acquiring of detection data acquired by a plurality of seismographs includes:

Acquire waveform data in seed format collected by a plurality of seismographs;

The seed format waveform data is converted into sac format waveform data as the detection data.

In some embodiments, performing cross-correlation calculation on the detection data collected by two adjacent seismographs to obtain Green's functions corresponding to the two adjacent seismographs includes:

Cross-correlation calculation is performed on the continuous waveform data in the detection data collected by two adjacent seismographs to obtain the Green's functions corresponding to the two adjacent seismographs.

In some embodiments, performing surface wave dispersion curve imaging according to each of the Green's functions and each of the position data to obtain a corresponding surface wave dispersion curve energy diagram includes:

extracting surface wave information from each of the Green's functions;

The surface wave dispersion curve is imaged according to each of the surface wave information and each of the position data by using the frequency-Bessel transform method to obtain an energy diagram corresponding to the surface wave dispersion curve.

In some embodiments, the step of obtaining a stratigraphic structure profile of the target coastal zone according to the inversion of each of the surface wave dispersion curve energy diagrams includes:

A two-dimensional shear wave velocity structure diagram of the target coastal zone is obtained as the stratigraphic structure profile diagram based on the inversion of each of the surface wave dispersion curve energy diagrams.

To achieve the above-mentioned purpose, another aspect of the embodiment of the present application provides a coastal zone stratigraphic structure detection device, the device comprising:

A data acquisition unit, used to acquire detection data collected by a plurality of seismographs, and to acquire position data of each of the seismographs; wherein each of the seismographs is sequentially and continuously deployed on the land side and the sea side of the target coastal zone at a preset interval;

A cross-correlation calculation unit, used for performing cross-correlation calculation on the detection data collected by two adjacent seismographs to obtain Green's functions corresponding to the two adjacent seismographs;

A surface wave imaging unit, used for imaging a surface wave dispersion curve according to each of the Green's functions and each of the position data, to obtain an energy diagram of the corresponding surface wave dispersion curve;

The structural inversion unit is used to obtain the stratigraphic structural profile of the target coastal zone by inverting the energy diagrams of the surface wave dispersion curves.

To achieve the above objective, another aspect of an embodiment of the present application provides an electronic device, the electronic device comprising a memory and a processor, the memory storing a computer program, and the processor implementing the above method when executing the computer program.

To achieve the above objective, another aspect of an embodiment of the present application provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and the computer program implements the above method when executed by a processor.

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

The present application obtains detection data collected by multiple seismometers and obtains the position data of each seismometer; wherein each seismometer is sequentially and continuously deployed on the land side and the ocean side of the target coastal zone at a preset spacing; performs cross-correlation calculation on the detection data collected by two adjacent seismometers to obtain the Green's function corresponding to the two adjacent seismometers; performs surface wave dispersion curve imaging according to each Green's function and each position data to obtain the corresponding surface wave dispersion curve energy map; and obtains the stratigraphic structure profile of the target coastal zone by inverting the energy map of each surface wave dispersion curve. The present application can obtain the shear wave velocity structure information of the shallow strata of the target coastal zone by deploying multiple seismometers on the land side and the ocean side of the target coastal zone, and uses the surface wave information in the detection data collected by the seismometers to perform surface wave dispersion curve imaging, and then obtains the stratigraphic structure profile of the target coastal zone by inverting the energy map of the surface wave dispersion curve, thereby realizing seamless land and sea joint measurement.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG1 is a schematic diagram of the structure of a seismograph provided in an embodiment of the present application;

FIG2 is a schematic diagram of a seismograph deployment provided in an embodiment of the present application;

FIG3 is a schematic diagram of another deployment of a seismograph provided in an embodiment of the present application;

FIG4 is a schematic diagram of a flow chart of a coastal zone stratigraphic structure detection method provided in an embodiment of the present application;

FIG5 is an example diagram of a seismograph deployment scenario provided in an embodiment of the present application;

FIG6 is a schematic diagram of data collected by a seismograph provided in an embodiment of the present application;

FIG7 is an energy diagram of a surface wave dispersion curve provided in an embodiment of the present application;

FIG8 is a cross-sectional view of a two-dimensional shallow stratum shear wave velocity structure in a coastal zone provided by an embodiment of the present application;

FIG9 is an example flow chart of a coastal zone stratigraphic structure detection method provided by an embodiment of the present application;

FIG10 is a schematic diagram of the structure of a coastal zone stratigraphic structure detection device provided in an embodiment of the present application;

FIG. 11 is a schematic diagram of the hardware structure of an electronic device provided in an embodiment of the present application.

Detailed ways

In order to make the purpose, technical solutions and advantages of the present application clearer, the present application is further described in detail below in conjunction with the accompanying drawings and examples. It should be understood that the specific embodiments described herein are only used to explain the present application and are not intended to limit the present application. When the following description refers to the 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 attached claims.

It is understood that the terms "first", "second", etc. used in this application can be used to describe various concepts in this article, but unless otherwise specified, these concepts are not limited by these terms. These terms are only used to distinguish one concept from another concept. For example, without departing from the scope of the embodiment of the present application, the first information may also be referred to as the second information, and similarly, the second information may also be referred to as the first information. Depending on the context, the words "if" and "if" as used herein can be interpreted as "at the time of" or "when" or "in response to determination".

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

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

Before describing the embodiments of the present application in detail, some related technologies involved in the embodiments of the present application are first described as follows:

The coastal zone refers to the zone where the sea and land interact, that is, the intertidal zone affected by the daily rise and fall of the tides and the sea-land transition zone of a certain range of land and shallow sea on both sides, which has a very special geological and tectonic background. The coastal zone is also one of the areas with the most active, concentrated and intensive development and utilization of human economic activities. Its environmental evolution is directly related to human living space, quality of life and sustainable development of society. In recent years, it has become a hot research area. my country's eastern coast is close to the sea, with many islands and a total coastline length of 32,000 kilometers. Systematically and comprehensively revealing the spatial distribution and environmental change characteristics of underground space resources in the coastal zone plays a very important role in coastal planning, development and environmental protection. With the improvement of my country's economic development level, the demand for coastal zone survey technology methods in terms of bedrock surface burial depth, underground velocity structure detection, and shallow layer detection accuracy and breadth has gradually increased. The existing shallow layer survey in the coastal zone needs to use two geophysical survey methods, land and ocean, respectively, and cannot be carried out simultaneously. Generally speaking, the land survey adopts the refraction/reflection seismic exploration method, and the ocean survey adopts the single-channel seismic or shallow layer profile exploration method.

Seismic exploration refers to the use of seismological principles to obtain information about underground structures and properties based on the propagation characteristics of seismic waves in underground media. In the process of seismic exploration, seismic wavelets are first excited by the source, and the time and amplitude of the seismic waves arriving at different locations are recorded by seismic instruments. Then, these data are processed and interpreted using the theories and methods of seismology. By analyzing the path of seismic waves propagating underground, parameters such as the thickness, velocity, and density of the strata can be inverted to understand the underground structure and rock characteristics. The seismic refraction/reflection method on land and the single-channel seismic method on the ocean are based on the above theory. The seismic refraction/reflection method obtains underground structures by extracting the refraction and reflection travel time information of seismic waves in underground media. Single-channel seismic is currently recognized as one of the commonly used methods for comprehensive marine geophysical surveys in the world. Its acquisition system mainly consists of four parts: recording system, receiving cable, excitation source, and navigation and positioning system. Conventional single-channel seismic acquisition adopts a one-transmit-one-receive mode. The mainstream source types include air guns and sparks. The air gun source is mainly GI guns, and the excitation frequency is relatively low (effective frequency band range 40-300Hz), so it has good formation penetration ability. The spark source has a higher excitation frequency (effective frequency band range 300-5000Hz), high resolution and significantly lower formation penetration ability than the air gun source. Shallow formation profiles mainly use the principle of echo sounding and the propagation and reflection characteristics of sound waves in seawater and seabed sediments to continuously detect the seabed sediment structure, thereby obtaining a more intuitive seabed shallow formation velocity structure profile.

However, the existing technology has the following shortcomings: For a long time, obtaining the geological structure and stratigraphic information of the shallow strata on both sides of the coast usually requires two different types of geophysical exploration methods, land and sea. Due to the limitations of exploration technology development and the limitations of land and sea exploration methods, as well as the irreversible impact of different technical means on the ecological environment of the coastal zone, the shallow strata detection technology in the coastal zone and some coastal waters has always been a technical problem in the fields of engineering geology and environmental engineering at home and abroad, and is also a research hotspot of marine exploration technology at home and abroad. Existing seismic exploration technologies, such as single-channel seismic and shallow stratum profiles, have different capabilities in revealing the stratigraphic structure of land and sea. At the same time, the transmission frequency of existing marine seismic exploration technology is too high and is not suitable for seismic imaging of the seabed with hard bottom. High-density electrical methods are limited by environmental and human interference. The ability of drilling to obtain shallow information is very limited. In areas with shallow water depths, ships cannot pass and cannot carry out marine survey operations. There is a blind spot for the last 1 kilometer of investigation in the coastal zone. In terms of data processing, there are also large differences in accuracy and resolution between land and marine seismic exploration data, and effective splicing cannot be achieved. Therefore, there is an urgent need for new non-invasive geophysical exploration technologies and methods that are non-destructive, anti-interference, non-disturbing, and easy to construct, and are applicable to both land and sea, so as to obtain continuous high-resolution underground geological structure distribution characteristic information in the coastal zone.

In summary, through a comprehensive analysis of the characteristics of the marine and land environments, detection technology, data processing, and output methods, the existing coastal shallow strata survey technology has the following problems: 1) The existing land/marine seismic exploration technology cannot directly achieve seamless land-sea joint surveys; 2) The existing marine seismic exploration technology has a too high transmission frequency and is not suitable for seismic imaging of the seabed with a hard bottom; 3) The two existing land/marine seismic exploration technologies have different abilities to reveal the stratigraphic structure, and due to differences in benchmarks, accuracy, and resolution, effective data splicing cannot be achieved.

Therefore, the embodiment of the present application provides a method, device, electronic device and storage medium for detecting the stratigraphic structure of a coastal zone. The scheme obtains the detection data collected by multiple seismometers and obtains the position data of each seismometer; wherein each seismometer is continuously arranged on the land side and the sea side of the target coastal zone according to a preset spacing; performs cross-correlation calculation on the detection data collected by two adjacent seismometers to obtain the Green's function corresponding to the two adjacent seismometers; performs surface wave dispersion curve imaging according to each Green's function and each position data to obtain the corresponding surface wave dispersion curve energy map; obtains the stratigraphic structure profile of the target coastal zone by inverting the energy map of each surface wave dispersion curve. The present application can obtain the shear wave velocity structure information of the shallow strata of the target coastal zone by deploying multiple seismometers on the land side and the sea side of the target coastal zone, using the surface wave information in the detection data collected by the seismometers to perform surface wave dispersion curve imaging, and then inverting the stratigraphic structure profile of the target coastal zone according to the energy map of the surface wave dispersion curve, thereby realizing seamless sea-land joint measurement.

The embodiment of the present application provides a coastal zone stratigraphic structure detection method, which relates to the field of geological exploration technology. The coastal zone stratigraphic structure detection method provided in the embodiment of the present application can be applied to a terminal, can also be applied to a server, and can also be software running in 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, and a car terminal, etc., but is not limited to this; the server side can be configured as an independent physical server, or it can be configured as a server cluster or distributed system composed of multiple physical servers, and can also be configured as a cloud server that provides 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 that implements the coastal zone stratigraphic structure detection method, etc., but is not limited to the above forms.

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 distributed computing environments, in which 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.

Before describing the embodiments of the present application in detail, the seismograph and the process of deploying the seismograph are described first.

The seismograph of this embodiment may be any type of seismograph capable of collecting seismic wave data. For example, this embodiment may adopt an optional seismograph as shown in FIG1 , which is a broadband node seismograph, including a land node seismograph and an ocean node seismograph. The upper part of FIG1 shows a land node seismograph, and the lower part of FIG1 shows an ocean node seismograph. The land node seismograph may include:

1——4G antenna interface;

2——Left: power and communication interface, right: external GNSS interface;

3——Acquisition storage unit & timing control unit;

4——Sensor signal conditioning unit;

5 – Protective housing;

6——Rechargeable large-capacity lithium battery pack;

7——seismometer core;

8——Coupling foot.

Marine node seismometers may include:

9——Power supply, GNSS interface or external hydrophone interface;

10 – hydrophone;

11——wire hole;

12——Circuit main board (collection storage unit & timing control unit & sensor signal conditioning unit);

13——Stainless steel pressure cabin;

14——Rechargeable large-capacity lithium battery pack;

15——seismometer core;

16——Coupling foot.

The two seismograph models in Figure 1 are TVG-63 land node seismograph (working frequency band: 10s-200Hz) and TDO-74N ocean node seismograph (10s-100Hz, can be connected to an external hydrophone). The land/sea node seismograph uses satellite timing, and the absolute time service accuracy is 10μs. The sampling rate of the seismograph is optional at 2000Hz, 1000Hz, 500Hz, 250Hz, 200Hz, 100Hz, and 50Hz. The specific sampling rate can be evaluated according to the formula f=v/4h, where f represents the main frequency, v represents the average shear wave velocity of the formation, and h represents the detection depth. For example, the underground shear wave velocity is estimated to be 300m/s. If it is necessary to detect a target with an underground burial depth of about 5m, according to the above formula, the main frequency is about 15Hz. The sampling rate of the seismograph is generally at least twice the main frequency, so the sampling rate of the seismograph should be set to at least 30Hz. The spacing of seismographs should be comparable to the formation thickness (or related to the depth of the research target). The spacing determines the ultimate resolution of the final inversion shear wave model, and the arrangement length determines the total detectable formation thickness.

Then, the land/sea node seismometers are placed one by one in different environments (hillsides, beaches and seabeds) to collect micro-seismic data (see Figure 2 for a schematic diagram, where the yellow seismometers are land node seismometers and the gray seismometers are ocean node seismometers). When placing land node seismometers, for better coupling, the seismometers can be buried in pits (or compacted to the ground without burying them). The deployment of ocean node seismometers requires the use of frogmen. Compared with the existing direct deployment method, frogmen can more accurately place node seismometers at the designed locations, providing guarantees for fine imaging of later data processing (see Figure 3 for a schematic diagram).

4 , an embodiment of the present application provides a method for detecting stratigraphic structures in a coastal zone, which may include but is not limited to steps S400 to S430, as follows:

S400: Acquire detection data collected by a plurality of seismographs, and acquire position data of each of the seismographs; wherein each of the seismographs is sequentially and continuously deployed on the land side and the sea side of the target coastal zone at a preset interval.

Specifically, in this embodiment, a certain number of land/sea node seismometers can be deployed in different environments such as hillsides, beaches, and seabeds through a specific deployment method, and each seismometer can be deployed at equal intervals. For example, a schematic diagram of the deployment position of each seismometer can be referred to Figure 5. When deploying the marine node seismometer, first connect them in series with ropes, and then use frogmen to carry them to the seabed for sequential deployment, and keep the seismometer posture horizontal during the deployment process.

After each seismograph is deployed, the present embodiment can obtain the micromotion data recorded by each seismograph as detection data, and obtain the positioning data collected by the positioning device. The present application embodiment is described by taking the seismographs deployed as shown in FIG5 as an example, and the collected micromotion data (i.e., detection data) can refer to FIG6, which is the micromotion data collected by the vertical component of the land/sea node seismograph for 24 hours, where the horizontal axis time is UTC time.

Optionally, the positioning device can use a Trimble R9s GPS receiver, which uses RTK technology and has an accuracy of up to centimeters, providing accurate positioning information for subsequent data processing and fine imaging. Different observation systems are set up according to different target coastal zones (including the setting of seismometer deployment spacing, sampling rate and other information), and the appropriate deployment time is selected.

As an optional implementation, S400 may include: acquiring waveform data in a seed format acquired by a plurality of seismographs;

The seed format waveform data is converted into sac format waveform data as the detection data.

In order to improve data accuracy and improve data processing efficiency, the embodiment of the present application may further include: preprocessing the detection data to obtain preprocessed data.

Furthermore, the steps of data preprocessing can be more specific as follows:

The detection data is subjected to at least one of instrument response removal, mean removal, trend removal, filtering, noise signal segmentation, time domain normalization, and spectrum whitening to obtain the preprocessed data.

S410: performing cross-correlation calculation on the detection data collected by two adjacent seismographs to obtain Green's functions corresponding to the two adjacent seismographs.

Further, S410 may include:

Cross-correlation calculation is performed on the continuous waveform data in the detection data collected by two adjacent seismographs to obtain the Green's functions corresponding to the two adjacent seismographs.

S420: Perform surface wave dispersion curve imaging according to each of the Green's functions and each of the position data to obtain a corresponding surface wave dispersion curve energy diagram.

Further, S420 may include:

extracting surface wave information from each of the Green's functions;

The surface wave dispersion curve is imaged according to each of the surface wave information and each of the position data by using the frequency-Bessel transform method to obtain an energy diagram corresponding to the surface wave dispersion curve.

Exemplarily, Figure 7 shows the energy diagrams of surface wave dispersion curves of land/sea node seismometers in different environments, which are, from top to bottom, the energy diagrams of surface wave dispersion curves of hillsides, beaches, and seabeds. The horizontal axis is the surface wave frequency, and the vertical axis is the phase velocity spectrum. The black lines in the energy diagrams of each surface wave dispersion curve in Figure 7 are manually picked fundamental-order dispersion curves.

S430: Obtaining a stratigraphic structure profile of the target coastal zone by inverting energy diagrams of each of the surface wave dispersion curves.

Further, S430 may include:

A two-dimensional shear wave velocity structure diagram of the target coastal zone is obtained as the stratigraphic structure profile diagram based on the inversion of each of the surface wave dispersion curve energy diagrams.

Exemplarily, FIG8 shows a two-dimensional shallow stratum shear wave velocity structure profile of a coastal zone, wherein the horizontal axis is distance and the vertical axis is depth. FIG8 is obtained by inverting the energy diagram of the surface wave dispersion curve shown in FIG7 .

The embodiments of the present application can solve the following technical problems:

1. Creatively propose the use of an observation system that combines sea/land node seismometers to achieve seamless sea-land joint measurement and solve the technical problem of the last 1 km survey blind spot in the coastal zone.

2. Expandability: Apply surface wave exploration technology to the sea, use the surface wave information in the micro-motion data recorded by the seismograph, extract the dispersion information of the surface wave for inversion, and solve the technical problem that land and ocean exploration data cannot be spliced together.

The embodiments of the present application have the following beneficial effects:

Based on the land/sea node seismometer, the surface wave exploration technology is used to achieve seamless connection of land and sea joint measurement, solving the technical problem of the last 1 km detection blind spot in the coastal zone. The embodiment of the present application can save manpower, material resources and financial resources for the detection of shallow strata structures in the coastal zone, solve the technical problem that ships cannot pass through shallow water areas and data collection is difficult, and has important guiding significance in prospecting, underground structure detection (basement, fracture), noise detection and positioning, disaster warning, etc., and has important practical significance for improving the operational investigation capabilities in shallow water areas.

Next, the solution of the embodiment of the present application will be introduced and explained in conjunction with specific application examples:

9 , this embodiment provides an example flow chart of a method for detecting stratigraphic structures in a coastal zone.

Specifically, this embodiment may include two processes: data collection and data processing.

Data collection: Still taking the deployment position shown in Figure 5 as an example, in a certain period of time, researchers deployed 22 stations in the coastal zone of an island, each station deployed a seismograph, and each seismograph was used as a station. The selection of the station has three different geological backgrounds, among which 5 stations (B01-B05) were deployed on the hillside, 7 stations (B06-B12) were deployed on the beach, and 10 stations (T01-T05, S01-S05) were deployed on the seabed. The data acquisition time of this embodiment is 1 day, the station spacing is 5m or 20m (laying different spacings can use a limited number of seafloor node seismographs to obtain shallow stratum structure detection imaging results of different water depths in a data acquisition process), the total length of the two-dimensional survey line (B01-S05) is about 142m, and the sampling rate of each seismograph is set to 50Hz. The positioning device of this embodiment is a Trimble R9s GPS receiver, which adopts RTK technology, and the positioning accuracy can reach centimeter level. The specific location and elevation of the station are shown in Table 1.

Table 1 Location data of each seismograph

Data processing: The original data collected by the land/sea node seismometers are in the international seed format file, and the micromotion data (i.e., detection data) collected by each seismometer is a data file generated every hour. First, the seed file format is converted into a visual sac file format, and then the micromotion data is preprocessed, including removing instrument response, removing mean, removing trend, filtering, noise signal segmentation, time domain normalization, and spectrum whitening; secondly, the continuous waveform records of the two stations in the micromotion data are cross-correlated to obtain the Green's function of the two station pairs; finally, the frequency-Bessel transform method (F-J method) is used to image the energy map of the surface wave dispersion curve, and the shear wave velocity structure map of the shallow layer below the survey line composed of each seismometer is obtained by inversion.

Compared with the existing coastal exploration methods, this embodiment has the advantages of being non-destructive, anti-interference, non-disturbance to residents, green economy, and low construction difficulty. It does not require the use of professional survey ships for operations, solves the technical problems of shallow water areas where ships cannot pass and data collection is difficult, and realizes seamless sea and land joint measurement.

10 , the embodiment of the present application further provides a coastal zone stratigraphic structure detection device, which can implement the above-mentioned coastal zone stratigraphic structure detection method, and the device includes:

A data acquisition unit, used to acquire detection data collected by a plurality of seismographs, and to acquire position data of each of the seismographs; wherein each of the seismographs is sequentially and continuously deployed on the land side and the sea side of the target coastal zone at a preset interval;

A cross-correlation calculation unit, used for performing cross-correlation calculation on the detection data collected by two adjacent seismographs to obtain Green's functions corresponding to the two adjacent seismographs;

A surface wave imaging unit, used for performing surface wave dispersion curve imaging according to each of the Green's functions and each of the position data, to obtain a corresponding surface wave dispersion curve energy diagram;

The structural inversion unit is used to obtain the stratigraphic structural profile of the target coastal zone by inverting the energy diagrams of the surface wave dispersion curves.

It can be understood that the contents of the above method embodiments are all applicable to the present device embodiments, the functions specifically implemented by the present device embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

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 the processor implements the above coastal zone stratigraphic structure detection method when executing the computer program. The electronic device can be any smart terminal including a tablet computer, a car computer, etc.

It can be understood that the contents of the above method embodiments are all applicable to the present device embodiments, the functions specifically implemented by the present device embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

Please refer to FIG. 11 , which schematically shows the hardware structure of an electronic device according to another embodiment. The electronic device includes:

The processor 1101 may be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (Application Specific Integrated Circuit, ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of the present application;

The memory 1102 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). The memory 1102 can store an operating system and other application programs. When the technical solution provided in the embodiment of this specification is implemented by software or firmware, the relevant program code is stored in the memory 1102, and the processor 1101 calls and executes the coastal zone stratigraphic structure detection method of the embodiment of this application;

Input/output interface 1103, used to implement information input and output;

The communication interface 1104 is used to realize the communication interaction between the device and other devices. The communication can be realized through a wired manner (such as USB, network cable, etc.) or a wireless manner (such as mobile network, WIFI, Bluetooth, etc.);

A bus 1105 that transmits information between various components of the device (e.g., the processor 1101, the memory 1102, the input/output interface 1103, and the communication interface 1104);

The processor 1101 , the memory 1102 , the input/output interface 1103 and the communication interface 1104 are connected to each other in communication within the device via the bus 1105 .

An embodiment of the present application also provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the above-mentioned coastal zone stratigraphic structure detection method is implemented.

It can be understood that the contents of the above method embodiments are all applicable to the present storage medium embodiments, the functions specifically implemented by the present storage medium embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

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

The embodiments described in the embodiments of the present application are intended 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 in the embodiments of the present application. Those skilled in the art will appreciate that with the evolution of technology and the emergence of new application scenarios, the technical solutions provided in the embodiments of the present application are also applicable to similar technical problems.

Those skilled in the art will appreciate 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 shown in the figures, or a combination of certain steps, or different steps.

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

Those skilled in the art will appreciate that all or some of the steps in the methods disclosed above, and the functional modules/units in the systems and devices may be implemented as software, firmware, hardware, or a suitable combination thereof.

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 are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

It should be understood that in the present application, "at least one (item)" means one or more, and "plurality" means two or more. "And/or" is used to describe the association relationship of associated objects, indicating that three relationships may exist. For example, "A and/or B" can mean: only A exists, only B exists, and A and B exist at the same time, where A and B can be singular or plural. The character "/" generally indicates that the objects associated before and after are in an "or" relationship. "At least one of the following" or similar expressions refers to any combination of these items, including any combination of single or plural items. For example, at least one 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, c can be single or multiple.

In the several embodiments provided in the present 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 schematic. For example, the division of the above units is only a logical function division. There may be other division methods in actual implementation, such as 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 mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

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

In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

If 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 this understanding, the technical solution of the present application, or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including multiple instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of various embodiments of the present application. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (Read-Only Memory, referred to as ROM), random access memory (Random Access Memory, referred to as RAM), disk or optical disk and other media that can store programs.

The preferred embodiments of the present application are described above with reference to the accompanying drawings, but the scope of the rights of the present application is not limited thereto. Any modification, equivalent substitution and improvement made by a person skilled in the art without departing from the scope and essence of the present application should be within the scope of the rights of the present application.

Claims (10)

1. A method for detecting stratigraphic structure in a coastal zone, characterized in that the method comprises: Acquire detection data collected by a plurality of seismographs, and acquire position data of each of the seismographs; wherein each of the seismographs is sequentially and continuously deployed on the land side and the sea side of the target coastal zone at a preset interval; Performing cross-correlation calculation on the detection data collected by two adjacent seismographs to obtain Green's functions corresponding to the two adjacent seismographs; Perform surface wave dispersion curve imaging according to each of the Green's functions and each of the position data to obtain a corresponding surface wave dispersion curve energy diagram; The stratigraphic structure profile of the target coastal zone is obtained by inverting the energy diagrams of the surface wave dispersion curves. 2. A method for detecting stratigraphic structure in a coastal zone according to claim 1, characterized in that before performing cross-correlation calculation on the detection data collected by two adjacent seismographs to obtain the Green's functions corresponding to the two adjacent seismographs, the method further comprises: Preprocessing the detection data to obtain preprocessed data; The cross-correlation calculation is performed on the detection data collected by two adjacent seismographs to obtain Green's functions corresponding to the two adjacent seismographs, including: A cross-correlation calculation is performed on the pre-processed data corresponding to two adjacent seismographs to obtain the Green's functions corresponding to the two adjacent seismographs. 3. A coastal zone stratigraphic structure detection method according to claim 2, characterized in that the preprocessing of the detection data to obtain preprocessed data comprises: The detection data is subjected to at least one of instrument response removal, mean removal, trend removal, filtering, noise signal segmentation, time domain normalization, and spectrum whitening to obtain the preprocessed data. 4. A coastal zone stratigraphic structure detection method according to claim 1, characterized in that the acquisition of detection data collected by multiple seismic instruments comprises: Acquire waveform data in seed format collected by a plurality of seismographs; The seed format waveform data is converted into sac format waveform data as the detection data. 5. A method for detecting stratigraphic structure in a coastal zone according to claim 1, characterized in that the cross-correlation calculation is performed on the detection data collected by two adjacent seismographs to obtain the Green's function corresponding to the two adjacent seismographs, comprising: Cross-correlation calculation is performed on the continuous waveform data in the detection data collected by two adjacent seismographs to obtain the Green's functions corresponding to the two adjacent seismographs. 6. A method for detecting stratigraphic structure in a coastal zone according to claim 1, characterized in that the imaging of surface wave dispersion curves is performed according to each of the Green's functions and each of the position data to obtain a corresponding surface wave dispersion curve energy diagram, comprising: extracting surface wave information from each of the Green's functions; The surface wave dispersion curve is imaged according to each of the surface wave information and each of the position data by using the frequency-Bessel transform method to obtain an energy diagram corresponding to the surface wave dispersion curve. 7. A method for detecting stratigraphic structure in a coastal zone according to any one of claims 1 to 6, characterized in that the stratigraphic structure profile of the target coastal zone is obtained by inverting the energy diagrams of the surface wave dispersion curves, comprising: A two-dimensional shear wave velocity structure diagram of the target coastal zone is obtained as the stratigraphic structure profile diagram based on the inversion of each of the surface wave dispersion curve energy diagrams. 8. A device for detecting stratigraphic structure in a coastal zone, characterized in that the device comprises: A data acquisition unit, used to acquire detection data collected by a plurality of seismographs, and to acquire position data of each of the seismographs; wherein each of the seismographs is sequentially and continuously deployed on the land side and the sea side of the target coastal zone at a preset interval; A cross-correlation calculation unit, used for performing cross-correlation calculation on the detection data collected by two adjacent seismographs to obtain Green's functions corresponding to the two adjacent seismographs; A surface wave imaging unit, used for performing surface wave dispersion curve imaging according to each of the Green's functions and each of the position data, to obtain a corresponding surface wave dispersion curve energy diagram; The structural inversion unit is used to obtain the stratigraphic structural profile of the target coastal zone by inverting the energy diagrams of the surface wave dispersion curves. 9. An electronic device, characterized in that the electronic device comprises a memory and a processor, the memory stores a computer program, and the processor implements the method according to any one of claims 1 to 7 when executing the computer program. 10. A computer-readable storage medium storing a computer program, wherein the computer program implements the method according to any one of claims 1 to 7 when executed by a processor.