Methods for measuring and locating marine station structures applicable to silty coastal zones

Abstract

This invention provides a method for measuring and locating marine station structures suitable for silty coastal zones. The method includes acquiring multi-source monitoring data for the silty coastal zone; extracting deformation feature points of the marine station structure from the multi-source monitoring data and generating planar displacement maps and elevation displacement maps of the feature points relative to a reference period for each monitoring cycle; calculating spatial offsets based on the planar displacement maps and elevation displacement maps to generate geological deformation distribution maps for each monitoring cycle and mapping them to generate a foundation stability coefficient matrix; and fusing the foundation stability coefficient matrix for each monitoring cycle with historical positioning benchmarks to generate structural safety positioning coordinates. This invention can improve the accuracy and anti-interference capability of structural displacement monitoring in complex geological environments.

Description

Technical Field

[0001] This invention relates to the field of coastal engineering monitoring and geological disaster early warning technology, and more specifically, to a method for measuring and locating the structure of a marine station applicable to silty coastal zones. Background Technology

[0002] In the field of marine monitoring, especially for the structural measurement and positioning of marine stations in silty coastal zones, has always been a technical challenge. Existing measurement and positioning methods mostly rely on traditional techniques such as GPS positioning and leveling. These techniques have many limitations when faced with the complex geological conditions of silty coastal zones. For example, GPS signals are easily interfered with in the complex environment of coastal zones, leading to a decrease in positioning accuracy; leveling requires frequent manual operations, and on the soft soil foundation of silty coastal zones, the stability of the measurement benchmark is difficult to guarantee, affecting the accuracy of the measurement results. In addition, these traditional methods often fail to effectively consider the impact of geological deformation on the marine station structure, making it difficult to monitor and assess the safety status of the marine station structure in real time.

[0003] In implementing the embodiments of the present invention, the prior art has at least the following problems or defects: it cannot accurately reflect the dynamic deformation characteristics of the marine station structure in the silty coastal zone, it is difficult to adapt to complex geological environment changes, and it cannot effectively integrate multi-source monitoring data to improve positioning accuracy and reliability. Summary of the Invention

[0004] This invention provides a method for measuring and locating the structure of a marine station suitable for silty coastal zones, comprising:

[0005] Acquire multi-source monitoring data for silty coastal zones;

[0006] Extract deformation feature points of the marine station structure from the multi-source monitoring data and generate planar displacement maps and elevation displacement maps of the feature points relative to the reference cycle for each monitoring period;

[0007] Based on the plane displacement map and the elevation displacement map, spatial offset is calculated to generate geological deformation distribution maps for each monitoring period and mapped to generate a foundation stability coefficient matrix;

[0008] The foundation stability coefficient matrix of each monitoring period is fused with historical positioning benchmarks to generate structural safety positioning coordinates.

[0009] As a further improvement to this application, the calculation of spatial offset based on the plane displacement map and the elevation displacement map includes:

[0010] The displacement values ​​of each measuring point in the plane displacement diagram of each monitoring cycle are corrected by tidal compensation with the plane offset of the corresponding measuring point in the reference cycle to obtain the plane offset of each measuring point in the current cycle.

[0011] The settlement value of each measuring point in the elevation displacement map of each monitoring period is corrected by sedimentation rate compensation with the elevation offset of the corresponding measuring point in the reference period to obtain the elevation offset of each measuring point in the current period.

[0012] As a further improvement to this application, the generation of geological deformation distribution maps for each monitoring period includes:

[0013] Based on the plane offset and elevation offset of each measuring point, the measuring points are divided into dynamic deformation units and static units through the geological activity threshold, and plane movement coefficient and elevation movement coefficient are set.

[0014] The composite deformation index is calculated based on the plane movement coefficient and elevation movement coefficient of each measuring point.

[0015] The geological time-series filtering of the composite deformation index is used to obtain the deformation distribution values ​​of each measuring point, which constitutes the geological deformation distribution map.

[0016] As a further improvement to this application, the method of dividing the measuring points into dynamic deformation units and static units based on geological activity thresholds includes:

[0017] When the planar offset of the measuring point is greater than or equal to the critical value of the planar displacement, it is determined to be a planar dynamic unit, and its planar movement coefficient is set to 1.2 times the critical value of the planar displacement.

[0018] When the plane offset of the measuring point is less than the critical value of the plane displacement, it is determined to be a static plane element, and its plane movement coefficient is set to the ratio of the measured displacement to the critical value.

[0019] When the elevation offset of the measuring point is greater than or equal to the elevation settlement threshold, it is determined to be a vertical dynamic unit, and its elevation movement coefficient is set to 1.5 times the elevation settlement threshold.

[0020] When the elevation offset of a measuring point is less than the elevation settlement threshold, it is determined to be a vertical static unit, and its elevation movement coefficient is set to the logarithmic transformation value of the measured settlement and the threshold.

[0021] As a further improvement to this application, the composite deformation index is calculated using the following formula based on the plane movement coefficient and elevation movement coefficient of each measuring point:

[0022]

[0023] in: For measuring points The composite deformation index; This is the planar translation coefficient; This is the elevation shift factor; These are the characteristic parameters of the silty soil layer (determined through field soil mechanics tests, reflecting the soil layer's response characteristics to planar displacement and elevation settlement).

[0024] As a further improvement to this application, the mapping to generate the foundation stability coefficient matrix includes:

[0025] The deformation distribution values ​​at each measuring point on the geological deformation distribution map are processed using the following formula to resist interference:

[0026]

[0027] in: This is the foundation stability coefficient; The value is the deformation distribution after filtering; Tidal cycle influencing factors; The density coefficient of the sediment; It is the geological memory decay constant; This represents the local soil density at the measuring point.

[0028] As a further improvement to this application, the tidal compensation correction for the plane offset is performed using the following formula:

[0029]

[0030] in: For the nth monitoring cycle, the measuring point Planar offset; This is the rheological compensation coefficient for silt; This represents the measured planar displacement value during the nth period. The duration of the nth monitoring cycle; For measuring points Tidal current field data; This is the plane offset for the (n-1)th period.

[0031] As a further improvement to this application, the depositional compensation correction for elevation offset is performed using the following formula:

[0032]

[0033] in: For the nth monitoring cycle, the measuring point Elevation offset; The degree of consolidation of the sediment; This represents the measured elevation settlement value for the nth period. Let be the compressibility modulus of the i-th type of sediment; This corresponds to the thickness variation of the i-th type of sediment layer; This represents the elevation offset for the (n-1)th cycle.

[0034] As a further improvement to this application, the formula for fusing and generating structural safety positioning coordinates is as follows:

[0035]

[0036] in: It represents the Hadamardi (or Hadama) stack; This represents a weighted fusion operator that considers the creep characteristics of soil. The coordinates are used to determine the safety positioning of the merged structure; This refers to the measured coordinate data for the current monitoring period; It is the optimal reference value in the historical safety coordinate database.

[0037] As a further improvement to this application, the generation of the planar displacement map and the elevation displacement map includes:

[0038] Four-dimensional sonar interferometry technology was used to obtain seabed topographic change data for a continuous monitoring period;

[0039] Strain field data of the structural foundation is collected through a distributed optical fiber sensing system; the strain field data is then calculated into planar displacement components and elevation displacement components.

[0040] An adaptive Kalman filter is used to suppress noise in the measurement area affected by wave interference, generating an anti-interference displacement map.

[0041] The embodiments of the present invention have at least the following beneficial effects:

[0042] 1. This invention, by acquiring multi-source monitoring data from silty coastal zones and extracting deformation characteristic points of marine station structures, can comprehensively reflect structural displacement under complex geological environments. Tidal compensation is used to correct planar offsets, and sedimentation rate compensation is used to correct elevation offsets, effectively eliminating the interference of tidal forces and sedimentation processes on the measurement results and significantly improving the accuracy of displacement monitoring. This method solves the technical challenges of traditional measurement techniques in silty coastal zones, where they are greatly affected by environmental factors and accumulate measurement errors, providing a reliable data foundation for subsequent structural safety assessments.

[0043] 2. This invention innovatively divides measuring points into dynamic deformation units and static units based on geological activity thresholds, and calculates a composite deformation index based on plane movement coefficients and elevation movement coefficients, achieving a refined assessment of foundation stability. By applying geological time-series filtering to the composite deformation index, interference from short-term geological activities is eliminated, accurately reflecting the long-term deformation trend of the foundation. This method solves the problem that traditional monitoring methods struggle to distinguish between dynamic deformation and static stability areas, providing a scientific basis for engineering safety early warning, and is particularly suitable for geological environments with significant rheological characteristics, such as silty soil layers.

[0044] 3. This invention intelligently fuses the foundation stability coefficient matrix from each monitoring period with historical positioning benchmarks, employing a weighted fusion operator that considers soil creep characteristics to dynamically generate structural safety positioning coordinates. This method effectively solves the problem of positioning benchmark drift caused by soil consolidation and creep during long-term monitoring, ensuring the long-term stability of marine station structures under complex geological conditions. Compared to traditional static positioning methods, this invention can adapt to changes in the geological environment, providing reliable technical support for the safe operation and maintenance of marine engineering projects in silty coastal zones. Attached Figure Description

[0045] The above and other objects, features, and advantages of exemplary embodiments of the present invention will become readily apparent from the following detailed description taken in conjunction with the accompanying drawings. Several embodiments of the invention are illustrated in the drawings by way of example, not limitation, in which:

[0046] Figure 1 This is a flowchart illustrating a method for measuring and locating the structure of a marine station in a silty coastal zone, as provided in an embodiment of the present invention. Detailed Implementation

[0047] The principles and spirit of the invention will now be described with reference to several exemplary embodiments. It should be understood that these embodiments are provided merely to enable those skilled in the art to better understand and implement the invention, and are not intended to limit the scope of the invention in any way. Rather, these embodiments are provided to make the invention more thorough and complete, and to fully convey the scope of the invention to those skilled in the art.

[0048] Those skilled in the art will recognize that embodiments of the present invention can be implemented as a system, apparatus, device, method, or computer program product. Therefore, the present invention can be specifically implemented in the following forms: entirely hardware, entirely software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.

[0049] It should be noted that the number of any elements in the accompanying drawings is for illustrative purposes only and not as a limitation, and any naming is for distinction only and has no limiting meaning.

[0050] The following is for reference. Figure 1 , Figure 1 This is a schematic flowchart illustrating a method for measuring and locating the structure of a marine station in a silty coastal zone, as provided in an embodiment of the present invention. Figure 1 As shown, a method for measuring and locating the structure of a marine station suitable for silty coastal zones includes:

[0051] S1. Obtain multi-source monitoring data for silty coastal areas;

[0052] S2. Extract the deformation feature points of the marine station structure from the multi-source monitoring data and generate a planar displacement map and an elevation displacement map of the feature points relative to the reference cycle for each monitoring cycle.

[0053] S3. Based on the plane displacement map and the elevation displacement map, calculate the spatial offset to generate a geological deformation distribution map for each monitoring period and map it to generate a foundation stability coefficient matrix;

[0054] S4. The foundation stability coefficient matrix of each monitoring period is fused with the historical positioning benchmark to generate structural safety positioning coordinates.

[0055] It should be noted that this invention proposes a method for measuring and locating marine station structures applicable to silty coastal zones, aiming to address the problems of insufficient positioning accuracy and inability to effectively assess structural safety under complex geological conditions in existing technologies. Multi-source monitoring data includes, but is not limited to, seabed topographic change data and structural foundation strain field data, which can comprehensively reflect the dynamic changes of the marine station structure. Multi-source monitoring data is a collection of data from different monitoring methods and sensors, with diverse sources including, but not limited to, sonar and fiber optic sensors, used to provide deformation information of the marine station structure in different dimensions. By extracting deformation feature points from these data, a plane displacement map and an elevation displacement map are generated, which respectively show the horizontal and vertical displacement of the marine station structure. The geological deformation distribution map is generated based on the plane displacement map and the elevation displacement map. By calculating the spatial offset of each measuring point, the spatial distribution of geological deformation can be more accurately reflected. The foundation stability coefficient matrix is ​​obtained by mapping the geological deformation distribution map; it can quantify the stability of the foundation and provide key data for the final generation of structural safety positioning coordinates. The structural safety positioning coordinates are the result of fusing the foundation stability coefficient matrix with historical positioning benchmarks. They integrate current geological conditions and historical data, and can accurately reflect the safe position of the marine station structure.

[0056] Furthermore, the multi-source monitoring data includes seabed topographic change data acquired by four-dimensional sonar interferometry and structural foundation strain field data collected by a distributed fiber optic sensing system. Four-dimensional sonar interferometry is an advanced seabed topographic monitoring method that provides high-precision seabed topographic change data, while the distributed fiber optic sensing system can monitor the strain of the structural foundation in real time, providing important data for structural safety assessment. Deformation feature points refer to points in the monitoring data that reflect the main deformation of the structure, extracted from the multi-source monitoring data using specific algorithms. Planar displacement maps and elevation displacement maps are obtained by analyzing the displacements of deformation feature points, representing the horizontal and vertical displacements of the structure, respectively. When calculating spatial offsets, tidal compensation corrections are applied to the displacement values ​​in the planar displacement map, and sedimentation rate compensation corrections are applied to the settlement values ​​in the elevation displacement map. These corrections improve the accuracy of the calculations. The geological deformation distribution map is generated based on the horizontal and vertical offsets of each measuring point. The measuring points are divided into dynamic and static deformation units using a geological activity threshold, and corresponding movement coefficients are set. These coefficients reflect the degree of deformation at each measuring point. Finally, the deformation distribution values ​​for each measuring point are obtained through geological time-series filtering, thus constructing the geological deformation distribution map.

[0057] Preferably, in the tidal compensation correction of the horizontal offset, the silt rheological compensation coefficient needs to be considered. This coefficient reflects the weight of the rheological properties of silty soil on the displacement correction, and its value ranges from 0 to 1. The correction formula also involves parameters such as measured horizontal displacement values, tidal velocity field data, and historical correction values. The accurate acquisition and reasonable setting of these parameters are crucial to the accuracy of the correction results. The sedimentary compensation correction of the elevation offset needs to consider the sediment consolidation coefficient, which reflects the weight of the degree of sediment consolidation on the settlement correction, and its value also ranges from 0 to 1. The correction process also needs to consider parameters such as the compressibility modulus of the sediment and the corresponding change in sedimentary layer thickness. These parameters reflect the compressibility characteristics of the soil layer and are very important for accurately calculating the elevation offset. When generating the geological deformation distribution map, the calculation of the composite deformation index involves the horizontal displacement coefficient and the elevation displacement coefficient, as well as the characteristic parameters of the silty soil layer. These parameters are determined through field soil mechanics tests and can reflect the response characteristics of the soil layer to horizontal displacement and elevation settlement. Through these detailed steps and parameter settings, spatial offset calculations and geological deformation distribution maps can be performed more accurately, thereby providing more reliable data support for the safe positioning of marine station structures.

[0058] In some embodiments, the calculation of spatial offset based on the planar displacement map and the elevation displacement map includes:

[0059] The displacement values ​​of each measuring point in the plane displacement diagram of each monitoring cycle are corrected by tidal compensation with the plane offset of the corresponding measuring point in the reference cycle to obtain the plane offset of each measuring point in the current cycle.

[0060] The settlement value of each measuring point in the elevation displacement map of each monitoring period is corrected by sedimentation rate compensation with the elevation offset of the corresponding measuring point in the reference period to obtain the elevation offset of each measuring point in the current period.

[0061] It should be noted that this invention pays particular attention to the impact of tides and sedimentation rates on the measurement results when calculating spatial offsets. Specifically, for displacement values ​​in the planar displacement map, tidal compensation correction is required because tidal forces cause periodic displacement changes in the marine station structure; without correction, the accuracy of displacement measurements would be affected. Similarly, for settlement values ​​in the elevation displacement map, sedimentation rate compensation correction is required because sediment consolidation and settlement affect elevation measurements. Here, tidal compensation and sedimentation rate compensation correction refer to adjusting the measured displacement and settlement values ​​using specific algorithms to eliminate errors caused by tides and sedimentation rates, thereby obtaining more accurate planar and elevation offsets.

[0062] Specifically, tidal compensation correction involves parameters including measured horizontal displacement values, tidal velocity field data, and historical correction values. Measured horizontal displacement values ​​refer to the horizontal displacement data of the marine station structure directly measured by monitoring equipment; tidal velocity field data refers to the tidal velocity information obtained from tidal models or on-site current meters during the monitoring period, reflecting the degree of tidal influence on the displacement of the marine station structure; historical correction values ​​refer to the horizontal offset that has been corrected in the previous monitoring period, used for further correction based on the measurement data of the current period. Sedimentation rate compensation correction involves parameters such as measured elevation settlement values, sediment compressibility modulus, and sediment thickness variation. Measured elevation settlement values ​​refer to the vertical settlement data of the marine station structure directly measured by monitoring equipment; sediment compressibility modulus refers to a physical parameter reflecting the compressibility characteristics of sediments, which is related to the type and properties of the sediments; sediment thickness variation refers to the change in sediment thickness during the monitoring period, which can be obtained through stratigraphic monitoring equipment. Accurate measurement and reasonable setting of these parameters are key to effective correction.

[0063] Preferably, when performing tidal compensation correction for planar offsets, a time-series-based algorithm can be used. This algorithm takes into account the measured planar displacement value of the current period, tidal velocity field data, and historical correction values ​​from the previous period. The algorithm first calculates the tidal influence of the current period, then adjusts the measured displacement value according to the silt rheological compensation coefficient, and finally combines the historical correction values ​​to obtain the corrected planar offset. For sedimentary compensation correction of elevation offsets, an algorithm based on sedimentary layer characteristics can be used. This algorithm takes into account the measured elevation settlement value of the current period, the compressibility modulus of the sediment, and the change in sedimentary layer thickness. The algorithm first calculates the compressibility of the sedimentary layer, then adjusts the measured settlement value according to the sediment consolidation coefficient, and finally combines the historical correction values ​​to obtain the corrected elevation offset. Through these specific algorithms and parameter settings, the accuracy of spatial offset calculations can be effectively improved, providing more reliable data support for subsequent geological deformation analysis and structural safety positioning.

[0064] In some embodiments, generating geological deformation distribution maps for each monitoring period includes:

[0065] Based on the plane offset and elevation offset of each measuring point, the measuring points are divided into dynamic deformation units and static units through the geological activity threshold, and plane movement coefficient and elevation movement coefficient are set.

[0066] The composite deformation index is calculated based on the plane movement coefficient and elevation movement coefficient of each measuring point.

[0067] The geological time-series filtering of the composite deformation index is used to obtain the deformation distribution values ​​of each measuring point, which constitutes the geological deformation distribution map.

[0068] It should be noted that this invention employs an analysis method based on the horizontal and vertical offsets of measuring points when generating the geological deformation distribution map. This method divides the measuring points into dynamic and static deformation units using a geological activity threshold, and sets horizontal and vertical movement coefficients. The geological activity threshold is a standard used to distinguish the degree of deformation activity of measuring points. Based on this standard, measuring points can be divided into dynamic and static deformation units, thus more accurately reflecting the geological deformation situation. The horizontal and vertical movement coefficients are parameters used to quantify the degree of deformation of measuring points in the horizontal and vertical directions. These parameters can be used to calculate a composite deformation index, thereby obtaining a geological deformation distribution map, providing an important basis for assessing foundation stability.

[0069] Specifically, geological activity thresholds include critical values ​​for horizontal displacement and elevation settlement. The critical value for horizontal displacement refers to the threshold value at which a measuring point is classified as a dynamic horizontal unit if its horizontal displacement exceeds this value; otherwise, it is classified as a static horizontal unit. The elevation settlement threshold refers to the threshold value at which a measuring point is classified as a dynamic vertical unit if its elevation displacement exceeds this value; otherwise, it is classified as a static vertical unit. The horizontal displacement coefficient and elevation displacement coefficient are parameters set based on the deformation of the measuring point, reflecting the degree of deformation. For dynamic horizontal units, the horizontal displacement coefficient is set to 1.2 times the critical value for horizontal displacement; for static horizontal units, it is set to the ratio of the measured displacement to the critical value. For dynamic vertical units, the elevation displacement coefficient is set to 1.5 times the elevation settlement threshold; for static vertical units, it is set to the logarithmic transformation of the measured settlement to the threshold. These parameter settings more accurately reflect the deformation of the measuring point, providing a basis for subsequent geological deformation analysis.

[0070] Preferably, a model based on a composite deformation index can be used when generating the geological deformation distribution map. This model takes into account the horizontal displacement coefficient and the vertical displacement coefficient of each measuring point, and calculates the composite deformation index using a specific algorithm. This algorithm considers the characteristic parameters of the silty soil layer, which are determined through in-situ soil mechanics tests and reflect the soil layer's response characteristics to horizontal displacement and vertical settlement. Specifically, the composite deformation index is calculated by multiplying the horizontal displacement coefficient and the vertical displacement coefficient by the corresponding characteristic parameters of the silty soil layer, then summing the results and taking the square root. The resulting composite deformation index undergoes geological time-series filtering to eliminate noise and outliers, ultimately yielding the deformation distribution values ​​for each measuring point, thus forming the geological deformation distribution map. This processing effectively improves the accuracy and reliability of the geological deformation distribution map, providing stronger support for foundation stability assessment.

[0071] In some embodiments, dividing the measuring points into dynamic deformation units and static units based on geological activity thresholds includes:

[0072] When the planar offset of the measuring point is greater than or equal to the critical value of the planar displacement, it is determined to be a planar dynamic unit, and its planar movement coefficient is set to 1.2 times the critical value of the planar displacement.

[0073] When the plane offset of the measuring point is less than the critical value of the plane displacement, it is determined to be a static plane element, and its plane movement coefficient is set to the ratio of the measured displacement to the critical value.

[0074] When the elevation offset of the measuring point is greater than or equal to the elevation settlement threshold, it is determined to be a vertical dynamic unit, and its elevation movement coefficient is set to 1.5 times the elevation settlement threshold.

[0075] When the elevation offset of a measuring point is less than the elevation settlement threshold, it is determined to be a vertical static unit, and its elevation movement coefficient is set to the logarithmic transformation value of the measured settlement and the threshold.

[0076] It should be noted that this invention employs specific judgment criteria and coefficient setting methods when classifying measuring points into dynamic deformation units and static units using geological activity thresholds. This method can more accurately identify the deformation state of measuring points, thus providing more reliable data support for subsequent geological deformation analysis. Specifically, when the plane offset of a measuring point is greater than or equal to the plane displacement critical value, it is determined to be a dynamic plane unit, and its plane movement coefficient is set to 1.2 times the plane displacement critical value; when the plane offset of a measuring point is less than the plane displacement critical value, it is determined to be a static plane unit, and its plane movement coefficient is set to the ratio of the measured displacement to the critical value. Similarly, when the elevation offset of a measuring point is greater than or equal to the elevation settlement threshold, it is determined to be a vertical dynamic unit, and its elevation movement coefficient is set to 1.5 times the elevation settlement threshold; when the elevation offset of a measuring point is less than the elevation settlement threshold, it is determined to be a vertical static unit, and its elevation movement coefficient is set to the logarithmic transformation value of the measured settlement and the threshold. These judgment criteria and coefficient settings can effectively distinguish the deformation state of measuring points, providing a basis for subsequent calculation of composite deformation index.

[0077] Specifically, the critical value for horizontal displacement and the elevation settlement threshold are key parameters used to distinguish the deformation state of a measuring point. The critical value for horizontal displacement refers to the threshold value in the horizontal direction; if the horizontal displacement of a measuring point exceeds this value, it is considered to be in a dynamic deformation state; otherwise, it is considered to be in a static deformation state. The elevation settlement threshold refers to the threshold value in the vertical direction; if the elevation displacement of a measuring point exceeds this value, it is considered to be in a dynamic deformation state; otherwise, it is considered to be in a static deformation state. The horizontal displacement coefficient and the elevation displacement coefficient are parameters used to quantify the degree of deformation of a measuring point. For dynamic horizontal elements, the horizontal displacement coefficient is set to 1.2 times the critical value for horizontal displacement, indicating that the deformation degree of the measuring point is relatively significant; for static horizontal elements, the horizontal displacement coefficient is set to the ratio of the measured displacement to the critical value, indicating that the deformation degree of the measuring point is relatively weak. For vertical dynamic units, the elevation movement coefficient is set to 1.5 times the elevation settlement threshold, indicating that the settlement at that measuring point is relatively significant. For vertical static units, the elevation movement coefficient is set to the logarithmic transformation value of the measured settlement and the threshold, indicating that the settlement at that measuring point is relatively weak. These parameter settings can more accurately reflect the deformation state of the measuring points, providing important basis for subsequent geological deformation analysis.

[0078] Preferably, a threshold-based determination model can be used when determining the deformation state of the measuring point. This model takes the plane offset and elevation offset of the measuring point as input and compares them with a critical value for plane displacement and a threshold value for elevation settlement to determine the deformation state of the measuring point. The specific steps are as follows: First, obtain the plane offset and elevation offset of the measuring point; then, compare the plane offset with the critical value for plane displacement. If it is greater than or equal to the critical value, it is determined to be a dynamic plane element; otherwise, it is determined to be a static plane element. Next, compare the elevation offset with the elevation settlement threshold. If it is greater than or equal to the threshold, it is determined to be a dynamic vertical element; otherwise, it is determined to be a static vertical element. Based on the determination results, set the corresponding plane movement coefficient and elevation movement coefficient. For example, for a dynamic plane element, the plane movement coefficient is set to 1.2 times the critical value for plane displacement; for a static plane element, the plane movement coefficient is set to the ratio of the measured displacement to the critical value. For a dynamic vertical element, the elevation movement coefficient is set to 1.5 times the elevation settlement threshold; for a static vertical element, the elevation movement coefficient is set to the logarithmic transformation value of the measured settlement and the threshold. By using this specific judgment model and parameter settings, the deformation state of the measuring points can be identified more accurately, providing more reliable data support for subsequent geological deformation analysis.

[0079] In some embodiments, the composite deformation index is calculated using the following formula based on the plane movement coefficient and elevation movement coefficient of each measuring point:

[0080]

[0081] in: For measuring points The composite deformation index; This is the planar translation coefficient; This is the elevation shift factor; These are the characteristic parameters of the silty soil layer (determined through field soil mechanics tests, reflecting the soil layer's response characteristics to planar displacement and elevation settlement).

[0082] It should be noted that this invention employs a method that comprehensively considers both horizontal and vertical displacement coefficients when calculating the composite deformation index. This method combines the horizontal and vertical displacement coefficients through a specific calculation formula to generate a composite deformation index that comprehensively reflects the deformation of the measuring point. The composite deformation index is a quantitative indicator used to assess the combined degree of deformation of the measuring point in both horizontal and vertical directions. This approach more accurately reflects the distribution of geological deformation, providing important evidence for subsequent foundation stability assessments.

[0083] Specifically, the calculation of the composite deformation index involves two key parameters: the horizontal displacement coefficient and the vertical displacement coefficient. The horizontal displacement coefficient reflects the degree of deformation of the measuring point in the horizontal direction, while the vertical displacement coefficient reflects the degree of deformation in the vertical direction. When calculating the composite deformation index, the characteristic parameters of the silty soil layer also need to be considered. These parameters are determined through field soil mechanics tests and reflect the soil layer's response characteristics to horizontal displacement and vertical settlement. These parameters include α, β, and γ, where α and β represent the weights of the horizontal displacement coefficient and the vertical displacement coefficient, respectively, and γ is used to adjust the scale of the final composite deformation index. The specific values ​​of these parameters need to be adjusted according to the actual geological conditions and monitoring data to ensure the accuracy and reliability of the calculation results.

[0084] Preferably, a weighted comprehensive model can be used to calculate the composite deformation index. The input parameters of this model include the horizontal and vertical movement coefficients of the measuring point, as well as the characteristic parameters α, β, and γ of the silty soil layer determined through field tests. The specific calculation steps are as follows: First, multiply the horizontal movement coefficient by the weight α, and multiply the vertical movement coefficient by the weight β; then, add the two results together; finally, take the γ-th root of this sum to obtain the composite deformation index. This calculation method can comprehensively consider the deformation of the measuring point in both the horizontal and vertical directions, and adjust the contribution of deformation in different directions through weight parameters, thereby obtaining a more comprehensive and accurate deformation assessment index. Through this detailed calculation process, it can be ensured that the composite deformation index can effectively reflect the actual deformation of the measuring point, providing reliable data support for the generation of geological deformation distribution maps.

[0085] In some embodiments, the mapping to generate the foundation stability coefficient matrix includes:

[0086] The deformation distribution values ​​at each measuring point on the geological deformation distribution map are processed using the following formula to resist interference:

[0087]

[0088] in: This is the foundation stability coefficient; The value is the deformation distribution after filtering; This refers to the tidal cycle influence factor (reflecting the periodic disturbance of tidal action on foundation deformation). This is the sediment density coefficient (a correction parameter related to soil compaction). It is the geological memory decay constant (characterizing the decay characteristics of soil deformation over time). The density of the local soil at the measuring point (obtained through borehole sampling or geophysical exploration).

[0089] It should be noted that this invention employs an anti-interference processing method when generating the foundation stability coefficient matrix. This method obtains the foundation stability coefficient by performing specific calculations on the deformation distribution values ​​at each measuring point of the geological deformation distribution map. This process considers various influencing factors, such as tidal cycle influence factors, sediment density coefficients, geological memory decay constants, and local soil density at the measuring points, to ensure that the calculation results accurately reflect the stability of the foundation. The foundation stability coefficient is an important quantitative indicator used to assess the stability of the foundation under current geological conditions, providing crucial data support for the safe positioning of marine station structures.

[0090] Specifically, the calculation of the foundation stability coefficient involves several key parameters. The tidal cycle influence factor reflects the periodic disturbance of foundation deformation caused by tidal forces; its value is typically adjusted based on historical tidal data and foundation response. The sediment density coefficient, related to soil compaction, is used to correct for deformation effects caused by changes in sediment density. The geological memory decay constant characterizes the decay characteristics of soil deformation over time, reflecting the long-term stability of foundation deformation. The local soil density at the measuring point is obtained through borehole sampling or geophysical methods, reflecting the actual density of the soil at that point. The proper setting and accurate acquisition of these parameters are crucial for the calculation of the foundation stability coefficient. In the calculation process, the deformation distribution values ​​at each measuring point on the geological deformation distribution map are first filtered to obtain the filtered deformation distribution values. Then, combined with the above parameters, the foundation stability coefficient is obtained through a specific calculation formula.

[0091] Preferably, a model that comprehensively considers multiple influencing factors can be used when generating the foundation stability coefficient matrix. The input parameters of this model include filtered deformation distribution values, tidal cycle influence factors, sediment density coefficient, geological memory decay constant, and local soil density at the measuring point. The specific calculation steps are as follows: First, obtain the deformation distribution values ​​of each measuring point on the geological deformation distribution map and filter them to eliminate the influence of noise and outliers. Then, divide the filtered deformation distribution values ​​by a comprehensive correction factor, which consists of the tidal cycle influence factor, sediment density coefficient, geological memory decay constant, and local soil density at the measuring point. In this way, a quantitative index that accurately reflects foundation stability—the foundation stability coefficient—can be obtained. This calculation method can effectively consider the influence of multiple geological and environmental factors on foundation stability, providing more reliable data support for the safe positioning of marine station structures.

[0092] In some embodiments, the tidal compensation correction for the plane offset is performed using the following formula:

[0093]

[0094] in: For the nth monitoring cycle, the measuring point Planar offset; This is the silt rheological compensation coefficient (with a value ranging from 0 to 1, reflecting the weight of silty soil rheological properties on displacement correction). This represents the measured planar displacement value during the nth period. The duration of the nth monitoring cycle; For measuring points Tidal current field data (obtained through tidal models or on-site current meters); This is the plane offset for the (n-1)th cycle (historical correction value).

[0095] It should be noted that this invention employs a comprehensive correction method that considers multiple factors when performing tidal compensation correction for planar offset. This method uses a specific algorithm to adjust the measured planar displacement value to eliminate the influence of tidal velocity field data on the measurement results. The correction process also considers the silt rheological compensation coefficient, an important parameter reflecting the weight of silty soil rheological properties on displacement correction. This correction method yields a more accurate planar offset, thereby improving the accuracy of marine station structural measurement and positioning.

[0096] Specifically, the tidal compensation correction for planar offset involves several key parameters. The measured planar displacement value refers to the horizontal displacement data of the marine station structure obtained directly from monitoring equipment. Tidal velocity field data refers to the tidal velocity information acquired by the tidal model or on-site current meter during the monitoring period, reflecting the degree of tidal influence on the marine station structure displacement. The silt rheological compensation coefficient is a value between 0 and 1, used to adjust the degree of influence of tidal velocity on displacement; its value is determined based on the rheological properties of silty soil. During the correction process, the tidal influence of the current period is first calculated, then the measured displacement value is adjusted according to the silt rheological compensation coefficient, and finally, the corrected planar offset is obtained by combining the historical correction value from the previous period. This correction method can effectively eliminate the interference of tides on the measurement results and improve the accuracy of the planar offset.

[0097] Preferably, a time-series-based correction model can be used to correct the plane offset using tidal compensation. The model's input parameters include the measured plane displacement value for the current period, tidal velocity field data, historical correction values ​​from the previous period, and a silt rheological compensation coefficient. The specific correction steps are as follows: First, calculate the tidal influence for the current period, which is the product of the tidal velocity field data and the monitoring period duration. Then, adjust the measured plane displacement value according to the silt rheological compensation coefficient to obtain the tidal-corrected displacement value. Finally, perform a weighted average of the tidal-corrected displacement value and the historical correction values ​​from the previous period to obtain the corrected plane offset. This correction model effectively considers the periodic influence of tides and the rheological characteristics of silty soil, thus obtaining a more accurate plane offset and providing more reliable data support for the safe positioning of marine station structures.

[0098] In some embodiments, the following formula is used to perform deposition compensation correction for elevation offset:

[0099]

[0100] in: For the nth monitoring cycle, the measuring point Elevation offset; This is the sediment consolidation coefficient (ranging from 0 to 1, reflecting the weight of sediment consolidation degree on settlement correction). This represents the measured elevation settlement value for the nth period. The compression modulus of the i-th type of sediment (a physical parameter reflecting the compressibility of the soil layer); This corresponds to the thickness variation of the i-th type of sedimentary layer (obtained through stratigraphic monitoring). This is the elevation offset for the (n-1)th cycle (historical correction value).

[0101] It should be noted that this invention employs a correction method that comprehensively considers both the degree of sediment consolidation and compressibility characteristics when performing sedimentary compensation correction for elevation offsets. This method uses a specific algorithm to adjust the measured elevation settlement value to eliminate the influence of sediment consolidation and compression on the measurement results. The correction process also considers the sediment consolidation coefficient, an important parameter reflecting the weight of sediment consolidation on settlement correction. This correction method yields more accurate elevation offsets, thereby improving the accuracy of marine station structural measurement and positioning.

[0102] Specifically, sedimentary compensation correction for elevation offset involves several key parameters. The measured elevation settlement value refers to the vertical settlement data of the marine station structure obtained directly from monitoring equipment. The sediment consolidation coefficient, a value between 0 and 1, is used to adjust the influence of sediment consolidation on settlement; its value is determined based on the consolidation characteristics of the sediment. The sediment compressibility modulus is a physical parameter reflecting the compressibility characteristics of sediments, which is related to the type and properties of the sediment. In the correction process, the compression of the sedimentary layer is first calculated, which is the product of the sediment compressibility modulus and the change in sedimentary layer thickness. Then, the measured settlement value is adjusted according to the sediment consolidation coefficient. Finally, the corrected elevation offset is obtained by combining the historical correction value from the previous period. This correction method effectively eliminates the interference of sediment consolidation and compression on the measurement results, improving the accuracy of the elevation offset.

[0103] Preferably, a correction model based on sedimentary layer characteristics can be used when performing sedimentary compensation correction of elevation offset. The input parameters of this model include the measured elevation settlement value for the current period, the compressibility modulus of the sediment, the change in sedimentary layer thickness, the historical correction value from the previous period, and the sediment consolidation coefficient. The specific correction steps are as follows: First, calculate the compressibility of the sedimentary layer, which is the product of the sedimentary compressibility modulus and the change in sedimentary layer thickness. Then, adjust the measured settlement value according to the sediment consolidation coefficient to obtain the sediment-corrected settlement value. Finally, perform a weighted average of the sediment-corrected settlement value and the historical correction value from the previous period to obtain the corrected elevation offset. This correction model can effectively consider the consolidation and compressibility characteristics of sediments, thereby obtaining a more accurate elevation offset and providing more reliable data support for the safe positioning of marine station structures.

[0104] In some embodiments, the formula for fusing and generating structural safety positioning coordinates is:

[0105]

[0106] in: This represents the Hadamard product (element-wise multiplication of matrices). This represents a weighted fusion operator that considers the creep characteristics of soil (the weight allocation is determined based on the soil creep model). The coordinates are used to determine the safety positioning of the merged structure; This refers to the measured coordinate data for the current monitoring period; It is the optimal reference value in the historical safe coordinate database (based on the coordinate benchmark of historical stable state).

[0107] It should be noted that this invention employs a method that comprehensively considers both measured coordinate data from the current monitoring period and historical safety coordinate data when fusing and generating structural safety positioning coordinates. This method uses a specific fusion algorithm to combine the foundation stability coefficient matrix with historical positioning benchmarks to generate structural safety positioning coordinates. These structural safety positioning coordinates are a crucial output, accurately reflecting the safe position of the marine station structure under current geological conditions, providing critical data support for the operation and maintenance of the marine station.

[0108] Specifically, the fusion of structural safety positioning coordinates involves several key parameters. The foundation stability coefficient matrix, calculated in the previous steps, reflects the foundation stability at each measuring point. The measured coordinate data for the current monitoring period refers to the coordinate data of the marine station structure acquired by measuring equipment within the current monitoring period. The optimal reference values ​​in the historical safety coordinate database refer to coordinate benchmarks based on historical stable states; these data are typically selected from stable coordinate data obtained from long-term monitoring. During the fusion process, the foundation stability coefficient matrix is ​​first processed, and then fused with the measured coordinate data for the current monitoring period and the historical safety coordinate data. The fusion algorithm involves the Hadamard product and a weighted fusion operator considering soil creep characteristics. These operations ensure that the fusion result comprehensively considers current and historical data, while also taking into account the soil creep characteristics.

[0109] Preferably, a weighted fusion model can be used when generating structural safety positioning coordinates. The input parameters of this model include the foundation stability coefficient matrix, measured coordinate data from the current monitoring period, and historical safety coordinate data. The specific fusion steps are as follows: First, the foundation stability coefficient matrix is ​​processed to obtain the stability coefficient for each measuring point. Then, the stability coefficient of each measuring point is multiplied by the measured coordinate data from the current monitoring period using a Hadamard product operation, i.e., corresponding elements are multiplied. Next, the result is weighted and fused with the historical safety coordinate data, with the weights determined based on the soil creep characteristics. This fusion model effectively considers both current and historical data, while also taking into account the soil creep characteristics, thereby obtaining more accurate structural safety positioning coordinates and providing more reliable data support for the safety assessment and positioning of marine station structures.

[0110] In some embodiments, generating the planar displacement map and the elevation displacement map includes:

[0111] Four-dimensional sonar interferometry technology was used to obtain seabed topographic change data for a continuous monitoring period;

[0112] Strain field data of the structural foundation is collected through a distributed optical fiber sensing system; the strain field data is then calculated into planar displacement components and elevation displacement components.

[0113] An adaptive Kalman filter is used to suppress noise in the measurement area affected by wave interference, generating an anti-interference displacement map.

[0114] It should be noted that this invention employs an implementation method that integrates multiple monitoring technologies and data processing methods when generating planar displacement maps and elevation displacement maps. This method acquires seabed topographic change data over a continuous monitoring period using four-dimensional sonar interferometry and collects structural foundation strain field data through a distributed fiber optic sensing system. After processing, this data generates interference-resistant displacement maps, thus providing accurate foundational data for subsequent geological deformation analysis. This method is particularly suitable for complex geological environments such as silty coastal zones, effectively improving the accuracy and reliability of measurements.

[0115] Specifically, four-dimensional sonar interferometry is a method for monitoring seabed topography, providing high-precision data on seabed topographic changes. This data reflects the topographic changes of the marine station foundation at different points in time, providing important information for generating planar displacement maps. Distributed fiber optic sensing systems are used to collect strain field data of the structural foundation. This data reflects the stress and strain within the structure and is crucial for assessing structural safety. The process of converting strain field data into planar and elevation displacement components involves algorithms that transform strain data into displacement data. For measurement areas affected by wave interference, adaptive Kalman filtering is used for noise suppression. This algorithm dynamically adjusts filtering parameters based on the characteristics of the measurement data, effectively reducing the interference of environmental factors such as waves on the measurement results and generating anti-interference displacement maps.

[0116] Preferably, a comprehensive data processing model can be used to generate the planar displacement map and the elevation displacement map. The input parameters of this model include seabed topographic change data acquired by four-dimensional sonar interferometry and structural foundation strain field data collected by a distributed fiber optic sensing system. The specific processing steps are as follows: First, the four-dimensional sonar interferometry data is preprocessed to extract key features of seabed topographic changes. Then, the strain field data is solved to obtain the planar displacement component and the elevation displacement component. Next, for measurement areas affected by wave interference, an adaptive Kalman filter algorithm is applied for noise suppression. The adaptive Kalman filter algorithm dynamically adjusts the filtering parameters according to the characteristics of the measurement data to achieve the best noise suppression effect. Finally, the processed data is used to generate an anti-interference displacement map, including the planar displacement map and the elevation displacement map. This comprehensive data processing model can effectively process multiple monitoring data, improve the accuracy and reliability of the displacement map, and provide a solid data foundation for subsequent geological deformation analysis and structural safety assessment.

[0117] The above embodiments of the present invention have the following beneficial effects:

[0118] 1. This invention, by acquiring multi-source monitoring data from silty coastal zones and extracting deformation characteristic points of marine station structures, can comprehensively reflect structural displacement under complex geological environments. Tidal compensation is used to correct planar offsets, and sedimentation rate compensation is used to correct elevation offsets, effectively eliminating the interference of tidal forces and sedimentation processes on the measurement results and significantly improving the accuracy of displacement monitoring. This method solves the technical challenges of traditional measurement techniques in silty coastal zones, where they are greatly affected by environmental factors and accumulate measurement errors, providing a reliable data foundation for subsequent structural safety assessments.

[0119] 2. This invention innovatively divides measuring points into dynamic deformation units and static units based on geological activity thresholds, and calculates a composite deformation index based on plane movement coefficients and elevation movement coefficients, achieving a refined assessment of foundation stability. By applying geological time-series filtering to the composite deformation index, interference from short-term geological activities is eliminated, accurately reflecting the long-term deformation trend of the foundation. This method solves the problem that traditional monitoring methods struggle to distinguish between dynamic deformation and static stability areas, providing a scientific basis for engineering safety early warning, and is particularly suitable for geological environments with significant rheological characteristics, such as silty soil layers.

[0120] 3. This invention intelligently fuses the foundation stability coefficient matrix from each monitoring period with historical positioning benchmarks, employing a weighted fusion operator that considers soil creep characteristics to dynamically generate structural safety positioning coordinates. This method effectively solves the problem of positioning benchmark drift caused by soil consolidation and creep during long-term monitoring, ensuring the long-term stability of marine station structures under complex geological conditions. Compared to traditional static positioning methods, this invention can adapt to changes in the geological environment, providing reliable technical support for the safe operation and maintenance of marine engineering projects in silty coastal zones.

[0121] Furthermore, the storage medium in the embodiments of this application stores program instructions capable of implementing all the above methods. These program instructions can be stored in the storage medium in the form of a software product, including several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) or processor to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks, or terminal devices such as computers, servers, mobile phones, and tablets.

[0122] The above description is merely an explanation of some preferred embodiments of the present invention and the technical principles employed. Those skilled in the art should understand that the scope of the invention as described in the embodiments of the present invention is not limited to technical solutions formed by specific combinations of the above-described technical features, but should also cover other technical solutions formed by arbitrary combinations of the above-described technical features or their equivalents without departing from the above-described inventive concept. For example, technical solutions formed by substituting the above-described features with (but not limited to) technical features with similar functions disclosed in the embodiments of the present invention.

Claims

1. A method for measuring and locating the structure of a marine station suitable for silty coastal zones, characterized in that, Includes the following steps: Acquire multi-source monitoring data for silty coastal zones; Extract deformation feature points of the marine station structure from the multi-source monitoring data and generate planar displacement maps and elevation displacement maps of the feature points relative to the reference cycle for each monitoring period; Based on the plane displacement map and the elevation displacement map, spatial offset is calculated to generate geological deformation distribution maps for each monitoring period and mapped to generate a foundation stability coefficient matrix; The foundation stability coefficient matrix of each monitoring period is fused with historical positioning benchmarks to generate structural safety positioning coordinates; The calculation of spatial offset based on the plane displacement map and the elevation displacement map includes: The displacement values ​​of each measuring point in the plane displacement diagram of each monitoring cycle are corrected by tidal compensation with the plane offset of the corresponding measuring point in the reference cycle to obtain the plane offset of each measuring point in the current cycle. The settlement values ​​of each measuring point in the elevation displacement map of each monitoring period are compared with the elevation offset of the corresponding measuring point in the reference period by sedimentation rate compensation correction to obtain the elevation offset of each measuring point in the current period. The mapping generates the foundation stability coefficient matrix, including: The deformation distribution values ​​at each measuring point on the geological deformation distribution map are processed using the following formula to resist interference: in: This is the foundation stability coefficient; The value is the deformation distribution after filtering; Tidal cycle influencing factors; The density coefficient of the sediment; It is the geological memory decay constant; This refers to the local soil density at the measuring point. Use the following formula to perform tidal compensation correction for plane offset: in: For the nth monitoring cycle, the measuring point Planar offset; This is the rheological compensation coefficient for silt; This represents the measured planar displacement value during the nth period. The duration of the nth monitoring cycle; For measuring points Tidal current field data; This is the plane offset for the (n-1)th period; The deposition rate compensation correction for elevation offset is performed using the following formula: in: For the nth monitoring cycle, the measuring point Elevation offset; The degree of consolidation of the sediment; This represents the measured elevation settlement value for the nth period. Let be the compressibility modulus of the i-th type of sediment; This corresponds to the thickness variation of the i-th type of sediment layer; This represents the elevation offset for the (n-1)th cycle. The formula for generating the structural safety positioning coordinates through fusion is: in: It represents the Hadamardi (or Hadama) stack; This represents a weighted fusion operator that considers the creep characteristics of soil. The coordinates are used to determine the safety positioning of the merged structure; This refers to the measured coordinate data for the current monitoring period; It is the optimal reference value in the historical safety coordinate database.

2. The method according to claim 1, characterized in that, The generation of geological deformation distribution maps for each monitoring period includes: Based on the plane offset and elevation offset of each measuring point, the measuring points are divided into dynamic deformation units and static units through the geological activity threshold, and plane movement coefficient and elevation movement coefficient are set. The composite deformation index is calculated based on the plane movement coefficient and elevation movement coefficient of each measuring point. The geological time-series filtering of the composite deformation index is used to obtain the deformation distribution values ​​of each measuring point, which constitutes the geological deformation distribution map.

3. The method according to claim 2, characterized in that, The method of dividing the measuring points into dynamic deformation units and static units based on geological activity thresholds includes: When the planar offset of the measuring point is greater than or equal to the critical value of the planar displacement, it is determined to be a planar dynamic unit, and its planar movement coefficient is set to 1.2 times the critical value of the planar displacement. When the plane offset of the measuring point is less than the critical value of the plane displacement, it is determined to be a static plane element, and its plane movement coefficient is set to the ratio of the measured displacement to the critical value. When the elevation offset of the measuring point is greater than or equal to the elevation settlement threshold, it is determined to be a vertical dynamic unit, and its elevation movement coefficient is set to 1.5 times the elevation settlement threshold. When the elevation offset of a measuring point is less than the elevation settlement threshold, it is determined to be a vertical static unit, and its elevation movement coefficient is set to the logarithmic transformation value of the measured settlement and the threshold.

4. The method according to claim 2 or 3, characterized in that, The composite deformation index is calculated using the following formula based on the plane displacement coefficient and elevation displacement coefficient of each measuring point: in: For measuring points The composite deformation index; This is the planar translation coefficient; This is the elevation shift factor; These are characteristic parameters of the silty soil layer, determined through field soil mechanics tests, reflecting the soil layer's response characteristics to planar displacement and elevation settlement.

5. The method according to claim 1, characterized in that, The generation of the planar displacement map and the elevation displacement map includes: Four-dimensional sonar interferometry technology was used to obtain seabed topographic change data for a continuous monitoring period; Strain field data of the structural foundation is collected through a distributed optical fiber sensing system; the strain field data is then calculated into planar displacement components and elevation displacement components. An adaptive Kalman filter is used to suppress noise in the measurement area affected by wave interference, generating an anti-interference displacement map.

Images