A joint processing method for land-sea interface zone comprehensive geophysical data

Abstract

The present application belongs to the technical field of data processing of land-sea interface zone, and particularly relates to a joint processing method of comprehensive geophysical data of land-sea interface zone, which is used for processing and interpreting data of geophysical prospecting methods at the land-sea interface zone, finding out geological structure features and distribution of adverse geological bodies in the study area, including data processing and comprehensive interpretation of different methods, wherein the method comprises A: conical field source transient electromagnetic method, B: geological radar data processing method, C: high-density resistivity method, and D: micro-motion data processing method, and the specific steps of the conical field source transient electromagnetic method are A1: data arrangement and preprocessing are performed on original measurement data, and the arranged data are converted into.TEM format, so that the data become editable data traces. The present application has higher sensitivity to faults and fracture zones, has strong adaptability to urban environment, and to some extent, overcomes the problem of multi-solution, and has high efficiency, low cost, and higher accuracy and reliability.

Description

Technical Field

[0001] This invention belongs to the field of land-sea boundary zone data processing technology, specifically relating to a method for joint processing of comprehensive geophysical data of the land-sea boundary zone. Background Technology

[0002] The land-sea boundary zone refers to the area where land and sea meet, including coastal plains, estuaries, coastlines, and coastal mountains. This area has complex geological structures, topography, hydrogeochemistry, and other characteristics, and is an important research object in resource development, environmental protection, urban construction, and other fields. In order to detect and evaluate the underground structure, resource distribution, and environmental conditions of the land-sea boundary zone, geophysical methods are required for exploration.

[0003] Geophysical methods utilize the interaction between geophysical fields (such as electromagnetic fields, seismic waves, and gravitational fields) and subsurface media. By measuring the distribution and changes of geophysical fields, the physical parameters and spatial distribution of the subsurface media can be deduced. Geophysical methods have advantages such as large detection depth, wide range, low cost, and non-destructive nature, making them one of the important tools in geological exploration.

[0004] However, geophysical exploration in the land-sea boundary zone faces many difficulties and challenges, mainly in the following aspects:

[0005] The exploration environment in coastal cities is complex, with various anthropogenic and natural sources of interference, such as power lines, pipelines, metal objects, groundwater, high-velocity layers on the surface, and urban vibrations, which affect the measurement and interpretation of geophysical fields.

[0006] The subsurface media in the land-sea boundary zone exhibit strong heterogeneity, nonlinearity, and instability, such as variations in soil and rock layers, the development of fault structures, and the influence of saltwater intrusion, leading to the complexity and uncertainty of the geophysical field.

[0007] The geophysical field at the land-sea interface is affected by factors such as the conductivity of seawater, changes in the coastline, and the effects of tides, resulting in variations and anomalies in the geophysical field.

[0008] Geophysical exploration in the land-sea boundary zone involves the coupling and transformation of various geophysical fields, such as electromagnetic fields and gravitational fields, seismic waves and electromagnetic waves, and requires consideration of various physical effects and interactions.

[0009] Therefore, in order to effectively detect and evaluate the underground structure, resource distribution, and environmental conditions of the land-sea boundary zone, it is necessary to adopt a reasonable combination of geophysical methods and data processing and interpretation methods to reduce environmental interference, improve coupling conditions, accurately display anomalies, and improve detection accuracy and reliability. Summary of the Invention

[0010] The purpose of this invention is to provide a joint processing method for comprehensive geophysical data of the land-sea boundary zone, which has high sensitivity to faults and fracture zones, strong adaptability to urban environments, and overcomes the problem of multiple solutions to a certain extent. It is efficient, low-cost, and has higher accuracy and reliability.

[0011] The specific technical solution adopted by this invention is as follows:

[0012] A method for joint processing of comprehensive geophysical data in the land-sea boundary zone, wherein geophysical data processing and interpretation are carried out at the land-sea boundary zone to ascertain the geological structural characteristics and distribution of adverse geological bodies in the study area, including data processing and comprehensive interpretation of different methods;

[0013] Data processing and comprehensive interpretation using various methods

[0014] The data processing for each method includes the following steps:

[0015] A: Transient electromagnetic method using a conical field source;

[0016] B: Ground-penetrating radar data processing method;

[0017] C: High-density resistivity method;

[0018] D: Micro-motion data processing method.

[0019] The specific steps of the transient electromagnetic method using the conical field source are as follows:

[0020] A1: Organize and preprocess the raw measurement data, and convert the organized data into .TEM format to make the data editable.

[0021] A2: Remove or correct interference from the data channels, and test the pre-processed data using Fraser filtering, Kalman filtering, and Kalrous filtering to determine the effective filtering method and select the corresponding filtering result.

[0022] A3: Based on the full delay theory, according to the instrument's recorded turn-off time, mutual inductance influence correction is performed within the full delay range, and effective data from the mid-to-late stage data channels is extracted.

[0023] A4: Calculate the mutual inductance coefficient of the calculation device and perform apparent resistivity and time-depth conversion;

[0024] A5: Comprehensively analyze the characteristics of the apparent resistivity curves at each measuring point, determine the characteristics of resistivity variation with depth, and perform data inversion;

[0025] A6: Import terrain data and perform elevation correction on depth data;

[0026] A7: Summarize all process and result data, generate graphs and output them, and interpret the data.

[0027] The specific steps of the ground-penetrating radar data processing method are as follows:

[0028] B1: Collect data and transfer it to an editable file, making the data an editable data track;

[0029] B2: Sort the data files, delete erroneous data, and ensure data validity;

[0030] B3: Delete useless data channels and perform terrain correction;

[0031] B4: Reduce DC drift;

[0032] B5: Use frequency filtering to suppress or eliminate interference waves in a certain frequency band;

[0033] B6: By offsetting, the reflection point is moved to its original position, and the diffracted wave is converged to the initial diffracted point, reducing the effective layered reflection wave of the same tilt angle.

[0034] B7: Removes tilted, layered interference waves;

[0035] B7: Suppressing multiple reflected waves using the deconvolution method;

[0036] B8: Perform spatial filtering to make anomalies have better continuity or independence;

[0037] B9: Store the processing results, identify in-phase axes and divide wave groups on the output map.

[0038] The specific steps of the high-density resistivity method are as follows:

[0039] C1: Collect data and perform quality correction on the raw data;

[0040] C2: Performs data denoising and filtering operations;

[0041] C3: Converts the format of the preprocessed data to generate a *.dat file that can be recognized by the Res2Dinv inversion software;

[0042] C4: Import the converted data into the Res2Dinv inversion software;

[0043] C5: Obtain the inversion result map, and interpret and analyze the inversion results in conjunction with geological conditions and geophysical data.

[0044] The specific steps of the micro-motion data processing method are as follows:

[0045] D1: Select exploration methods based on the exploration task, site conditions, and instrument performance;

[0046] D2: The distribution of the detector array is determined by the choice of exploration method;

[0047] D3: After data acquisition, artificial source detection is processed using the SASW and MASW methods;

[0048] D4: After the initial data processing, the dispersion curve is further extracted to obtain the frequency-velocity curve;

[0049] D5: Invert the frequency-velocity curve to obtain the depth-velocity curve and form a wave velocity profile.

[0050] D6: Analyze and interpret the final results.

[0051] The specific steps of the comprehensive interpretation are as follows:

[0052] a: Transforming geological problems into geophysical problems;

[0053] b: Geophysical anomalies reflected in geological phenomena when using geophysical exploration methods;

[0054] c: Transform the obtained geophysical data into geological language or diagrams, give it geological meaning, and affirm its geological characteristics.

[0055] In A5, the data is inverted using the quasi-geomagnetic damped least squares method.

[0056] In B7, the tilt angle filter fk is used to remove tilted, layered interference waves.

[0057] The technical effects achieved by this invention are as follows:

[0058] This invention introduces the conical transient electromagnetic method, whose theory, data processing technology and integrated transmission and reception equipment are all independently developed. Research results show that this method has high sensitivity to faults and fracture zones and strong adaptability to urban environments.

[0059] This invention utilizes multiple geophysical methods to comprehensively process and interpret data, and cross-verifies multiple data, thus overcoming the problem of multiple solutions to a certain extent.

[0060] The data processing and interpretation method of the present invention is highly efficient, low-cost, and has higher accuracy and reliability. Attached Figure Description

[0061] Figure 1 This is a schematic diagram of the logging curve types of the present invention;

[0062] Figure 2 This is a flowchart of the data processing for the transient electromagnetic method using a conical field source in this invention.

[0063] Figure 3 This is a flowchart of the ground-penetrating radar data processing in this invention;

[0064] Figure 4 This is a flowchart of the high-density resistivity method data processing in this invention;

[0065] Figure 5 This is a flowchart of the micro-motion data processing in this invention;

[0066] Figure 6 This is a flowchart of the integrated interpretation process for geophysical exploration in this invention;

[0067] Figure 7 This is a diagram showing the location of the measuring line in this invention;

[0068] Figure 8 This is a comprehensive geophysical result map of the survey lines in this invention. Detailed Implementation

[0069] To make the objectives and advantages of this invention clearer, the invention will be specifically described below with reference to embodiments. It should be understood that the following text is merely used to describe one or more specific embodiments of the invention and does not strictly limit the scope of protection specifically claimed by the invention.

[0070] like Figure 1-8 As shown, a combined processing method for integrated geophysical data in the land-sea boundary zone includes data processing and interpretation of integrated geophysical methods such as the conical field source transient electromagnetic method, ground-penetrating radar, high-density resistivity method, and micro-motion. Through comprehensive data processing and interpretation, the geological structural characteristics and distribution of adverse geological bodies in the study area are identified, enabling detailed exploration of coastal cities. The specific content is divided into two parts:

[0071] The data processing and comprehensive interpretation methods are as follows:

[0072] The methods include:

[0073] Conical Field Source Transient Electromagnetic Method:

[0074] Most of the data processing flow for conical transient electromagnetics (CTEM) is consistent with that of traditional CTEM, including data processing, preprocessing, interference removal, and filtering. However, based on traditional CTEM data processing, the CTEM processing flow adds mutual inductance removal to correct the inductance effects at the transmitting and receiving ends caused by the mutual inductance between the transmitting and receiving coils. At the same time, in terms of apparent resistivity and depth conversion, the sloping step current response over the entire sampling delay range is considered, and a new apparent resistivity calculation formula is adopted to convert the severely low apparent resistivity value of small devices to the normal order of magnitude.

[0075] The conical transient electromagnetic data processing flow is as follows: Figure 2 As shown, the specific process is as follows:

[0076] First, the raw measurement data is organized and preprocessed, and then the organized data is converted into .TEM format;

[0077] Second, the STEM2.0 data processing system independently developed by the research group was used to remove or correct interfering data channels. Fraser filtering, Kalman filtering and Kalrous filtering were performed on the pre-processed data to determine the effective filtering method and select the corresponding filtering result.

[0078] Third, based on the full-delay theory of inclined step current excitation, mutual inductance influence correction is performed within the full-delay range according to the instrument-recorded turn-off time, and effective data of the middle and late stages are extracted.

[0079] IV. Based on the determined parameters of the transmitting and receiving devices, the mutual inductance coefficient of the devices is calculated, and the apparent resistivity and time-depth conversion is performed using the precise apparent resistivity conversion technology with inductive effect developed by the research group.

[0080] V. By comprehensively analyzing the characteristics of the apparent resistivity curves at each measuring point, the resistivity variation characteristics with depth were preliminarily determined, and the data inversion was carried out using the quasi-magnetic damped least squares method.

[0081] 6. Import terrain data and perform elevation correction on depth data;

[0082] 7. Summarize all process and result data, generate graphs and output them, and interpret the data.

[0083] Ground-penetrating radar data processing method:

[0084] Ground-penetrating radar (GPR) data processing involves processing digitized electromagnetic wave reflection signals to suppress random and regular interference. Its purpose is to suppress interfering waves, highlight effective reflected waves, enhance radar images, and extract parameters such as electromagnetic wave velocity, amplitude, and waveform from the image profile for geological interpretation.

[0085] The first step involves data acquisition, data transmission, and file editing. Following this, the raw data files are inspected and preprocessed. The purpose of file inspection is to ensure the data is in a suitable format for processing. Preprocessing typically includes data cleaning and outlier removal. Preprocessing removes errors, missing information, or duplicate data to ensure accuracy and consistency.

[0086] Data processing methods include deleting useless channels, gain adjustment, terrain correction, frequency filtering, FK tilt filtering, deconvolution, offset, spatial filtering, etc., as shown in the appendix. Figure 3 ;

[0087] The specific process is as follows:

[0088] First, collect the data and transfer it to an editable file, making the data an editable data track;

[0089] Second, sort the data files, delete erroneous data, and ensure data validity;

[0090] 3. Delete useless roads and perform terrain correction;

[0091] 4. DC drift reduction: If the radar data is all positive, all negative, or the positive and negative half-cycles are asymmetrical, it indicates that there is DC drift in the signal. The average value is obtained by summing the channel data and dividing by the number of sampling points. Then, the average value is subtracted from the channel data to obtain the DC drift reduction result.

[0092] 5. Use frequency filtering to suppress or eliminate interference waves in a certain frequency band. Filtering is used to eliminate the influence of environmental or system noise to improve data quality. Commonly used one-dimensional time-domain filtering methods include high-pass filtering, low-pass filtering, band-pass filtering, etc., but in practice, it is necessary to avoid "over-processing" the raw data to ensure the authenticity of the data interpretation.

[0093] VI. When there are effective layered reflected waves with the same tilt angle, the reflection point is moved to its original position by offsetting and the diffracted wave is converged to the initial diffracted point to improve the resolution of the profile.

[0094] 7. When it is determined that there are no effective layered reflected waves with the same tilt angle, fk tilt angle filtering is used to remove the interference waves of the tilted layer. If there are effective layered reflected waves with the same tilt angle, time offset or depth offset processing methods are used to return the interface of the tilted layer reflected waves to their original position and converge the diffracted waves.

[0095] 8. Suppress multiple reflected waves using deconvolution. However, deconvolution is not suitable when the reflected signal is weak or the data signal-to-noise ratio is low.

[0096] 9. Perform spatial filtering to make anomalies have better continuity or independence;

[0097] 10. Store the processing results, identify the in-phase axes and divide the wave groups on the result map.

[0098] High-density resistivity method:

[0099] High-density resistivity data processing technology is relatively mature. Res2Dinv software is commonly used for data stitching, bad sector removal, and inversion. Its main processing flow is as follows: Figure 4 As shown;

[0100] Its main processing flow is as follows:

[0101] 1. Collect data and perform quality correction on the raw data to eliminate bad data that may be caused by grounding problems or electrode polarization instability, so as to obtain high-quality data for subsequent processing.

[0102] Second, perform data denoising and filtering operations to remove environmental noise, improve data quality, and ensure the reliability of the inversion.

[0103] Third, the preprocessed data is converted to a format that generates a *.dat file that can be recognized by the Res2Dinv inversion software; during this process, terrain data also needs to be added to better describe the underground structure.

[0104] Fourth, import the converted data into Res2Dinv inversion software. During import, set the inversion parameters, including finite element mesh generation and model layer thickness factors. For each anomaly, select appropriate data constraints and model constraints based on the actual situation. This helps improve the accuracy and stability of the inversion results.

[0105] V. Through the above steps, high-quality inversion result maps are obtained, and the inversion results are interpreted and analyzed in conjunction with geological conditions and other geophysical data.

[0106] Micro-motion data processing method:

[0107] The data acquisition and processing steps of micro-motion technology are similar to those of active source seismic exploration methods, consisting of three steps: data acquisition, dispersion curve extraction, and dispersion curve inversion.

[0108] Due to the dispersion characteristics of surface waves, micromotion signals exhibit significant variations in amplitude and frequency over time and space, but they satisfy statistical stability within a certain spatiotemporal range and can be described by a stationary random process. Unlike traditional seismic exploration and seismology that uses ray theory to estimate seismic wave propagation velocity, due to the uncertainty of the micromotion source, the phase velocity of surface waves in micromotion signals can be obtained by calculating the spatial autocorrelation coefficients between stations in the observation system array.

[0109] Micro-motion data processing flow as follows Figure 5 As shown, the specific process is as follows:

[0110] First, the appropriate exploration method must be selected based on the exploration task, site conditions, and instrument performance. Natural sources are less expensive, environmentally friendly, and can be observed for extended periods, but their signal-to-noise ratio is low. Artificial sources have a higher signal-to-noise ratio and faster acquisition speed, but they are more expensive and may cause vibration and noise to the surrounding environment. Therefore, the exploration method must be selected based on the actual situation.

[0111] 2. The distribution of the detector array is determined by the choice of exploration method. Artificial sources are usually arranged in a straight line for data acquisition, while natural sources have more arrangement options, such as circular, straight, and L-shaped.

[0112] III. After data acquisition, artificial source detection is processed using SASW (Spectral Analysis of Surface Waves) and MASW (Multichannel Analysis of Surface Waves) methods. SASW typically uses a single seismic source and a single seismic receiver; MASW is suitable for multiple seismic receivers, usually arranged linearly, to record surface waves simultaneously. This allows for the capture of multiple frequency components at the same time. Natural sources are typically processed and analyzed using SPAC (Spatial Autocorrelation), fk (FK) methods, and NCF (Noise Cross-correlation Function).

[0113] IV. After the initial data processing, the dispersion curve is further extracted to obtain the frequency-velocity curve;

[0114] 5. Invert the frequency-velocity curve to obtain the depth-velocity curve and form a wave velocity profile.

[0115] VI. Finally, the results are analyzed and the stratigraphy is interpreted.

[0116] The comprehensive interpretation is as follows:

[0117] When using geophysical methods to solve a geological task, two transformations are required: first, the geological problem must be transformed into a geophysical problem; second, geophysical methods must be used to investigate the geophysical anomalies reflected in the geological phenomena; and finally, the obtained geophysical data must be transformed into geological language or diagrams, given geological meaning, and its geological effects affirmed. Each geophysical method has its application conditions and scope. Because the geological and geophysical characteristics of mineral deposits and natural geographical conditions often vary from place to place, the effectiveness of methods or the limitations of certain geophysical methods in solving geological problems are affected. Geophysical observation data contains rich and varied geological information, but it may also be affected by various interfering factors or contain human observation errors.

[0118] The comprehensive geophysical interpretation in this application is based on data processing using various geophysical methods, combined with existing geological and borehole data within the study area, to rationally transform the results of the data processing into geological bodies, such as... Figure 6 ;

[0119] Different geophysical methods measure different underground physical parameters, and the interfaces of physical parameters do not necessarily correspond to geological interfaces. This increases the difficulty of comprehensive comparison of multiple geophysical methods. To effectively combine multiple geophysical methods, the most reasonable geological interpretation results should be given, and the various geophysical methods should corroborate each other to accurately determine the stratigraphic distribution and geological structural characteristics, thus providing a theoretical basis for the exploration of adverse geological bodies.

[0120] The comprehensive interpretation method proposed in this application can be applied within a burial depth of 100m and has been verified. The exploration results are good, and the geological bodies after interpretation are in good agreement with the known geological data. The stratification is clear and the geological structural features are also obvious. In the process of comprehensive interpretation of multiple geophysical data, the characteristics and detection capabilities of each geophysical method are combined, and the method of depth stratification and the combination of primary and secondary interpretation targets are adopted to carry out comprehensive geophysical interpretation.

[0121] On the one hand, the ground-penetrating radar method has high shallow resolution, and its profile mainly reflects geological information at a depth of less than 5m, which is used to reveal the thickness limit of the surface concrete layer and the location of pipelines, etc.; the conical transient electromagnetic method has high lateral resolution, and its profile mainly reflects geological information in the depth range of 3-50m, focusing on revealing the location of faults and fracture zones and water-bearing properties from the perspective of lateral uniformity, while providing auxiliary evidence for stratigraphic interface identification and bedrock weathering degree judgment; the micro-motion exploration method has strong vertical detection capability, and its profile reflects a depth range of 5-100m.

[0122] It is mainly used to identify stratigraphic interfaces and to supplement information on the extension of anomalies reflected in transient electromagnetic profiles below 50m.

[0123] Finally, by combining the interpretation results of various geophysical methods, the stratigraphic distribution and geological structure characteristics of the study area can be obtained. In areas with abundant data, the accuracy of the detection results can also be verified by combining known geological and borehole data.

[0124] In view of the above, a specific embodiment of this application is as follows:

[0125] A geologically rich area within the Jiaodong Bay was selected as the study area, and a survey line with good exploration conditions was chosen within the study area. Geophysical exploration using multiple methods was carried out along this survey line, and the data was processed and interpreted comprehensively.

[0126] The location of the survey line is as follows Figure 7 As shown;

[0127] The survey line is laid out along Ruihu Road, with a length of 870m and a total of 174 survey points. Figure 8 The images show the measurement results of several geophysical methods used for this survey line. The radar imagery reveals a distinct reflective interface above the survey line (above the red line), presumably composed of Quaternary fill soil.

[0128] In the transient electromagnetic resistivity and fretting transverse wave velocity profile, the longitudinal difference between resistivity and wave velocity reflects a three-layer structure in depth. The first layer, with a bottom interface depth of about 15m, mainly consists of miscellaneous fill, dredged fill, and some silty clay.

[0129] The second stratum is unconformable, with a bottom interface depth of about 30m. The strata interface is clear and mainly consists of silty clay, serving as the bedrock surface that separates the Quaternary strata from the underlying conglomerate.

[0130] The third layer lies below the bedrock surface and consists mainly of strongly and moderately weathered conglomerate. Local areas below the bedrock surface (horizontal distance 480-550m) are presumed to be bedrock fracture zones.

[0131] After data processing and comprehensive interpretation according to the present invention, the stratigraphic distribution and structural characteristics of the survey line are clear and understandable, verifying the validity of this application.

[0132] The above description is merely a preferred embodiment of the present invention. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principles of the present invention, and these improvements and modifications should also be considered within the scope of protection of the present invention. Structures, devices, and operating methods not specifically described or explained in this invention are implemented according to conventional methods in the art unless otherwise specified or limited.

Claims

1. A method for joint processing of integrated geophysical data of the land-sea boundary zone, characterized in that: The method includes: data processing and comprehensive interpretation of each method; The specific steps of the comprehensive interpretation are as follows: a: Transforming geological problems into geophysical problems; b: Geophysical anomalies reflected in geological phenomena when using geophysical exploration methods; c: Transform the obtained geophysical data into geological language or diagrams, give it geological meaning, and affirm its geological characteristics; The data processing for each method includes the following steps: A: Transient electromagnetic method using a conical field source; The specific steps of the transient electromagnetic method using the conical field source are as follows: A1: Organize and preprocess the raw measurement data, and convert the organized data into .TEM format to make the data editable. A2: Remove or correct interference from the data channels, and test the pre-processed data using Fraser filtering, Kalman filtering, and Kalrous filtering to determine the effective filtering method and select the corresponding filtering result. A3: Based on the full-delay theory of inclined step current excitation, mutual inductance influence correction is performed within the full-delay range according to the instrument-recorded turn-off time, and effective data of the middle and late-stage data channels are extracted. A4: Calculate the mutual inductance coefficient of the calculation device and perform apparent resistivity and time-depth conversion; A5: Comprehensively analyze the characteristics of the apparent resistivity curves at each measuring point, determine the characteristics of resistivity variation with depth, and perform data inversion; A6: Import terrain data and perform elevation correction on depth data; A7: Summarize all process and result data, create graphs and output them, and interpret the data; B: Ground-penetrating radar data processing method; C: High-density resistivity method; D: Micro-motion data processing method.

2. The method for joint processing of integrated geophysical data of the land-sea boundary zone according to claim 1, characterized in that: The specific steps of the ground-penetrating radar data processing method are as follows: B1: Collect data and transfer it to an editable file, making the data an editable data track; B2: Sort the data files, delete erroneous data, and ensure data validity; B3: Delete useless data channels and perform terrain correction; B4: Reduce DC drift; B5: Use frequency filtering to suppress or eliminate interference waves in a certain frequency band; B6: Determine if there are effective layered reflected waves with the same tilt angle. If so, move the reflection point to its original position by offsetting the reflection and converge the diffracted wave to the initial diffraction point to reduce the number of effective layered reflected waves with the same tilt angle. If not, use fk tilt angle filtering to remove the interference waves of tilted layers. B7: Suppressing multiple reflected waves using the deconvolution method; B8: Perform spatial filtering; B9: Store the processing results, identify in-phase axes and divide wave groups on the output map.

3. The method for joint processing of integrated geophysical data of the land-sea boundary zone according to claim 1, characterized in that: The specific steps of the high-density resistivity method are as follows: C1: Collect data and perform quality correction on the raw data; C2: Performs data denoising and filtering operations; C3: Converts the format of the preprocessed data to generate a format recognizable by Res2Dinv inversion software. document; C4: Import the converted data into the Res2Dinv inversion software; C5: Obtain the inversion result map, and interpret and analyze the inversion results in conjunction with geological conditions and geophysical data.

4. The method for joint processing of integrated geophysical data of the land-sea boundary zone according to claim 1, characterized in that: The specific steps of the micro-motion data processing method are as follows: D1: Select exploration methods based on the exploration task, site conditions, and instrument performance; D2: The distribution of the detector array is determined by the choice of exploration method; D3: After data acquisition, artificial source detection is processed using the SASW and MASW methods; D4: After the initial data processing, the dispersion curve is extracted to obtain the frequency-velocity curve; D5: Invert the frequency-velocity curve to obtain the depth-velocity curve and form a wave velocity profile. D6: Analyze and interpret the final results.

5. The method for joint processing of integrated geophysical data of the land-sea boundary zone according to claim 1, characterized in that: In A5, the data is inverted using the quasi-geomagnetic damped least squares method.

6. The method for joint processing of integrated geophysical data of the land-sea boundary zone according to claim 2, characterized in that: In B6, if there are effective layered reflected waves with the same tilt angle, time offset or depth offset processing methods are used to reposition the tilted layer reflected wave interface and converge the diffracted wave.

Images