Ocean station structure measuring and positioning method suitable for muddy coastal zone

Abstract

The invention provides an ocean station structure measuring and positioning method suitable for a muddy coastal zone. The ocean station structure measuring and positioning method comprises the following steps: acquiring multi-source monitoring data of a muddy coastal zone area; extracting deformation feature points of the ocean station structure in the multi-source monitoring data, and generating a plane displacement diagram and an elevation displacement diagram of each monitoring period relative to the feature points of the reference period; performing space offset calculation based on the plane displacement diagram and the elevation displacement diagram to generate a geological deformation distribution diagram of each monitoring period, and mapping to generate a foundation stability coefficient matrix; and fusing the foundation stability coefficient matrix of each monitoring period with a historical positioning reference to generate a structure safety positioning coordinate. According to the invention, the accuracy and anti-interference capability of structure displacement monitoring in a complex geological environment can be improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of coastal engineering monitoring and geological disaster early warning, more particularly, the present application relates to a marine station structure measurement positioning method suitable for a muddy coastal zone. BACKGROUND

[0002] In the field of marine monitoring, especially for the measurement positioning of marine station structures in muddy coastal zones, it has always been a technical challenge. Existing measurement positioning methods rely heavily on traditional GPS positioning, leveling and other technologies. These technologies have many limitations when faced with the complex geological conditions of muddy coastal zones. For example, GPS signals are easily disturbed in the complex environment of the coastal zone, resulting in a decrease in positioning accuracy; leveling requires frequent manual operation, and on the soft soil foundation of the muddy coastal zone, the stability of the measurement reference point 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 marine station structures, making it difficult to monitor and evaluate the safety status of marine station structures in real time.

[0003] In the process of implementing the embodiments of the present application, there are at least the following problems or defects in the prior art: it is difficult to accurately reflect the dynamic deformation characteristics of the marine station structure in the muddy coastal zone, it is difficult to adapt to complex changes in the geological environment, and it is difficult to effectively fuse multi-source monitoring data to improve positioning accuracy and reliability. SUMMARY

[0004] The present application provides a marine station structure measurement positioning method suitable for a muddy coastal zone, comprising: obtaining multi-source monitoring data of the muddy coastal zone region; extracting deformation feature points of the marine station structure in the multi-source monitoring data and generating a plane displacement map and an elevation displacement map of each monitoring period relative to the reference period feature points; generating a geological deformation distribution map of each monitoring period based on the plane displacement map and the elevation displacement map and mapping to generate a foundation stability coefficient matrix; fusing the foundation stability coefficient matrix of each monitoring period with the historical positioning reference to generate a structure safety positioning coordinate.

[0005] As a further improvement of the present application, the spatial offset calculation based on the plane displacement map and the elevation displacement map comprises: performing tidal compensation correction on the displacement values of each measuring point in the plane displacement map of each monitoring period and the plane offset of the corresponding measuring point in the reference period to obtain the plane offset of each measuring point in the current period; performing deposition rate compensation correction on the settlement values of each measuring point in the elevation displacement map of each monitoring period and 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.

[0006] As a further improvement of the present application, generating a geological deformation distribution map 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 the 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 composite deformation index is subjected to geological time series filtering to obtain deformation distribution values ​​of each measuring point, thereby forming the geological deformation distribution map.

[0007] As a further improvement of the present application, the method of dividing the measuring points into dynamic deformation units and static units by using geological activity thresholds includes: When the plane offset of the measuring point is greater than or equal to the plane displacement critical value, it is determined to be a plane dynamic unit, and its plane movement coefficient is set to 1.2 times the plane displacement critical value; When the plane offset of the measuring point is less than the critical value of the plane displacement, it is determined to be a plane static unit, and its plane movement coefficient is set as the ratio of the measured displacement to the critical value; 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 conversion value of the measured settlement and the threshold.

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

[0009] in: For measuring points Composite deformation index; is the plane movement coefficient; is the elevation movement coefficient; It is the characteristic parameter of the silty soil layer (determined through on-site soil mechanics tests, reflecting the response characteristics of the soil layer to plane displacement and elevation settlement).

[0010] As a further improvement of the present application, the mapping generates a foundation stability coefficient matrix, including: The deformation distribution value of each measuring point in the geological deformation distribution map is processed with the following formula for anti-interference:

[0011] in: is the foundation stability coefficient; is the deformation distribution value after filtering; is the tidal cycle influencing factor; is the sediment density coefficient; is the geological memory decay constant; is the local soil density at the measuring point.

[0012] As a further improvement of this application, the following formula is used to perform tidal compensation correction of the plane offset:

[0013] in: The measurement point of the nth monitoring cycle The plane offset of is the silt rheological compensation coefficient; is the measured plane displacement value in the nth period; is the duration of the nth monitoring cycle; For measuring points Tidal velocity field data; is the plane offset of the n-1th period.

[0014] As a further improvement of this application, the following formula is used to perform sedimentation compensation correction for elevation offset:

[0015] in: The measurement point of the nth monitoring cycle The elevation offset of is the sediment consolidation coefficient; is the measured elevation settlement value in the nth period; is the compression modulus of the i-th type of sediment; is the thickness change of the corresponding sedimentary layer of type i; is the elevation offset of the n-1th cycle.

[0016] As a further improvement of the present application, the formula for fusion-generated structural safety positioning coordinates is:

[0017] in: represents the Hadamard product; represents the weighted fusion operator considering the creep characteristics of soil; Securely locate the coordinates of the fused structure; is the measured coordinate data of the current monitoring period; It is the optimal reference value in the historical safety coordinate database.

[0018] As a further improvement of the present application, generating the plane displacement map and the elevation displacement map includes: Adopt four-dimensional sonar interference measurement technology to obtain seabed topography change data in continuous monitoring period; The strain field data of the structure foundation are collected by a distributed optical fiber sensing system, and the strain field data are calculated into plane displacement components and elevation displacement components; Adaptive Kalman filtering is used for noise suppression on the measurement area disturbed by waves to generate an anti-interference displacement atlas.

[0019] The above embodiments of the present application have at least the following beneficial effects: 1. The present application can comprehensively reflect the structure displacement under complex geological environment by obtaining multi-source monitoring data of the muddy coastal zone and extracting the deformation feature points of the ocean station structure. The tidal compensation is used to correct the plane offset, and the deposition rate compensation is used to correct the elevation offset, which effectively eliminates the interference of tidal action and deposition process on the measurement results, and significantly improves the accuracy of displacement monitoring. The method solves the technical problems of large environmental factor influence and measurement error accumulation of traditional measurement technology in the muddy coastal zone, and provides a reliable data basis for subsequent structure safety evaluation.

[0020] 2. The present application innovatively divides the measuring points into dynamic deformation units and static units by the geological activity threshold, and calculates the composite deformation index based on the plane movement coefficient and the elevation movement coefficient, which realizes the fine evaluation of the foundation stability. The short-term geological activity interference is eliminated by the geological time sequence filtering processing of the composite deformation index, and the long-term deformation trend of the foundation is accurately reflected. The method solves the problem that the traditional monitoring means cannot distinguish between dynamic deformation and static stable area, and provides a scientific basis for engineering safety warning, especially suitable for geological environment such as muddy soil layer with significant rheological properties.

[0021] 3. The present application dynamically generates the structure safety positioning coordinates by intelligently fusing the foundation stability coefficient matrix of each monitoring period with the historical positioning reference, and using the weighted fusion operator considering the creep characteristics of soil. This method effectively solves the problem of positioning reference drift caused by soil consolidation, creep and other factors in the long-term monitoring process, and ensures the long-term stability of the ocean station structure under complex geological conditions. Compared with the traditional static positioning method, the present application can adapt to the change of geological environment, and provides reliable technical support for the safe operation and maintenance of the ocean engineering in the muddy coastal zone. BRIEF DESCRIPTION OF DRAWINGS

[0022] The above and other objects, features and advantages of the exemplary embodiments of the present application will be more apparent from the following detailed description read in conjunction with the accompanying drawings, in which: Figure 1A flowchart of a method for measuring and positioning a marine station structure suitable for a muddy coastal zone according to an embodiment of the present application is shown in FIG. 1. DETAILED DESCRIPTION

[0023] The principles and spirits of the present application will be described below with reference to several exemplary embodiments. It should be understood that these embodiments are given only so that those skilled in the art can better understand and implement the present application, and do not limit the scope of the present application in any way. On the contrary, these embodiments are provided so that the present application is more thorough and complete, and the scope of the present application is fully conveyed to those skilled in the art.

[0024] Those skilled in the art understand that the embodiments of the present application can be implemented as a system, a device, an apparatus, a method or a computer program product. Therefore, the present application can be embodied as a complete hardware, a complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.

[0025] It should be noted that any number of elements in the drawings are used for illustration only and not limitation, and any naming is only for distinction and does not have any limiting meaning.

[0026] Reference will be made to the following Figure 1 , Figure 1 A flowchart of a method for measuring and positioning a marine station structure suitable for a muddy coastal zone according to an embodiment of the present application is shown in FIG. 1. Figure 1 As shown in FIG. 1, a method for measuring and positioning a marine station structure suitable for a muddy coastal zone includes: S1, acquiring multi-source monitoring data of a muddy coastal zone area; S2, extracting deformation feature points of a marine station structure in the multi-source monitoring data and generating a planar displacement map and an elevation displacement map of feature points of each monitoring period relative to a reference period; S3, calculating spatial offset based on the planar displacement map and the elevation displacement map to generate a geological deformation distribution map of each monitoring period and map to generate a foundation stability coefficient matrix; S4, fusing the foundation stability coefficient matrix of each monitoring period with a historical positioning reference to generate a structure safety positioning coordinate.

[0027] It should be noted that the present application proposes a method for measuring and positioning the structure of a marine station suitable for a muddy coastal zone, aiming to solve the problems of insufficient positioning accuracy and inability to effectively evaluate the safety of the structure under complex geological conditions in the prior art. The multi-source monitoring data, including but not limited to seabed topography change data and structure foundation strain field data, can comprehensively reflect the dynamic changes of the marine station structure. Multi-source monitoring data is a collection of data from different monitoring means and sensors, which are diverse, including but not limited to sonar and optical fiber sensors, to provide deformation information of the marine station structure in different dimensions. By extracting the deformation feature points in these data, plane displacement maps and elevation displacement maps are generated, which respectively show the displacement of the marine station structure in the horizontal and vertical directions. The generation of the geological deformation distribution map is based on the plane displacement map and the elevation displacement map, and by calculating the spatial offset of each measuring point, it can more accurately reflect the spatial distribution of geological deformation. The foundation stability coefficient matrix is obtained by mapping the geological deformation distribution map, which can quantify the stability of the foundation and provide key data for the final structure safety positioning coordinate generation. The structure safety positioning coordinate is the result of fusing the foundation stability coefficient matrix with the historical positioning benchmark, which integrates the current geological conditions and historical data, and can accurately reflect the safe position of the marine station structure.

[0028] Further, the multi-source monitoring data includes seabed topography change data obtained by four-dimensional sonar interferometric measurement technology, and structure foundation strain field data collected by distributed optical fiber sensing system. Four-dimensional sonar interferometric measurement technology is an advanced seabed topography monitoring method that can provide high-precision seabed topography change data, while the distributed optical fiber sensing system can monitor the strain of the structure foundation in real time, providing important basis for structure safety evaluation. Deformation feature points are points in the monitoring data that can reflect the main deformation of the structure, which are extracted from multi-source monitoring data by a specific algorithm. Plane displacement maps and elevation displacement maps are obtained by analyzing the displacement of deformation feature points, which respectively represent the displacement of the structure in the horizontal and vertical directions. When calculating the spatial offset, the displacement values in the plane displacement map need to be tide compensation corrected, and the settlement values in the elevation displacement map need to be sedimentation rate compensation corrected, which can improve the accuracy of the calculation. The generation of the geological deformation distribution map is based on the plane offset and elevation offset of each measuring point, and by setting the corresponding movement coefficients, these coefficients reflect the deformation degree of the measuring point. Finally, through the geological time series filtering processing, the deformation distribution values of each measuring point are obtained, forming the geological deformation distribution map.

[0029] Preferably, in the process of tidal compensation correction of the plane offset, the silt rheological compensation coefficient needs to be considered, which reflects the weight of the silt rheological properties on displacement correction, and its value range is 0 to 1. The correction formula also involves parameters such as measured plane displacement value, tidal flow velocity field data and historical correction value, and the accurate acquisition and reasonable setting of these parameters are crucial to the accuracy of the correction result. The deposition compensation correction of the elevation offset needs to consider the sediment consolidation degree coefficient, which reflects the weight of the sediment consolidation degree on settlement correction, and also has a value range of 0 to 1. In the correction process, parameters such as the compression modulus of the sediment and the corresponding sediment layer thickness change need to be considered, which can reflect the compression characteristics of the soil layer and are very important for accurate calculation of the elevation offset. In generating the geological deformation distribution map, the calculation of the composite deformation index involves the plane movement coefficient and the elevation movement coefficient, as well as the silt layer characteristic parameters, which are determined through field soil mechanics test and can reflect the response characteristics of the soil layer to plane displacement and elevation settlement. Through these detailed steps and parameter settings, more accurate spatial offset calculation and geological deformation distribution map generation can be achieved, thereby providing more reliable data support for the safe positioning of the ocean station structure.

[0030] In some embodiments, the spatial offset calculation based on the plane displacement map and the elevation displacement map comprises: The displacement values of each monitoring period in the plane displacement map are subjected to tidal compensation correction with the plane offset of the corresponding measuring point in the reference period to obtain the plane offset of each measuring point in the current period. The settlement values of each monitoring period in the elevation displacement map are subjected to deposition rate compensation correction 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.

[0031] It should be noted that the present application pays special attention to the influence of tides and deposition rate on measurement results when calculating spatial offset. Specifically, for the displacement values in the plane displacement map, tidal compensation correction is needed, because tidal action will cause periodic displacement changes of the ocean station structure, which will affect the accuracy of displacement measurement if not corrected. Similarly, for the settlement values in the elevation displacement map, deposition rate compensation correction is needed, because the consolidation and settlement of sediments will affect the elevation measurement. The tidal compensation correction and deposition rate compensation correction here refer to adjusting the measured displacement and settlement values through a specific algorithm to eliminate the errors caused by tides and deposition rate, so as to obtain more accurate plane offset and elevation offset.

[0032] Specifically, the parameters involved in the tidal compensation correction include the measured plane displacement value, the tidal current velocity field data, and the historical correction value, etc. The measured plane displacement value refers to the displacement data of the ocean station structure in the horizontal direction directly measured by the monitoring equipment; the tidal current velocity field data refers to the tidal current velocity information obtained by the tidal model or the field current meter in the monitoring period, which reflects the influence degree of the tide on the displacement of the ocean station structure; the historical correction value refers to the plane displacement that has been corrected in the last monitoring period, which is used to further correct the measurement data of the current period. The deposition rate compensation correction involves the measured elevation subsidence value, the compression modulus of the sediment, and the thickness change amount of the sediment layer, etc. The measured elevation subsidence value refers to the subsidence data of the ocean station structure in the vertical direction directly measured by the monitoring equipment; the compression modulus of the sediment refers to the physical parameter reflecting the compression characteristics of the sediment, which is related to the type and nature of the sediment; the thickness change amount of the sediment layer refers to the thickness change of the sediment layer in the monitoring period, which can be obtained by the stratum monitoring equipment. The accurate measurement and reasonable setting of these parameters are the key to effective correction.

[0033] Preferably, when performing tidal compensation correction of the plane displacement, a time series-based algorithm can be used, which inputs the measured plane displacement value of the current period, the tidal current velocity field data, and the historical correction value of the last period. The algorithm first calculates the tidal influence amount of the current period, then adjusts the measured displacement value according to the silt rheological compensation coefficient, and finally obtains the corrected plane displacement combined with the historical correction value. For the deposition compensation correction of the elevation displacement, a deposition layer characteristic-based algorithm can be used, which inputs the measured elevation subsidence value of the current period, the compression modulus of the sediment, and the thickness change amount of the sediment layer. The algorithm first calculates the compression amount of the sediment layer, then adjusts the measured subsidence value according to the sediment consolidation degree coefficient, and finally obtains the corrected elevation displacement combined with the historical correction value. Through these specific algorithms and parameter settings, the accuracy of the spatial displacement calculation can be effectively improved, providing more reliable data support for subsequent geological deformation analysis and structure safety positioning.

[0034] In some embodiments, the generating the geological deformation distribution map of each monitoring period comprises: dividing the measuring points into dynamic deformation units and static units based on the plane displacement and the elevation displacement of each measuring point and setting the plane movement coefficient and the elevation movement coefficient; calculating the composite deformation index based on the plane movement coefficient and the elevation movement coefficient of each measuring point; performing geological time series filtering processing on the composite deformation index to obtain the deformation distribution value of each measuring point, constituting the geological deformation distribution map.

[0035] It should be noted that the present application adopts an analysis method based on the plane displacement and the elevation displacement of the measuring points when generating the geological deformation distribution map. This method divides the measuring points into dynamic deformation units and static units through a geological activity threshold, and sets the plane movement coefficient and the elevation movement coefficient. The geological activity threshold is a standard for distinguishing the activity degree of the deformation of the measuring points. Based on this standard, the measuring points can be divided into dynamic deformation units and static units, so as to more accurately reflect the geological deformation. The plane movement coefficient and the elevation movement coefficient are parameters for quantifying the degree of deformation of the measuring points in the horizontal and vertical directions. Through these parameters, the composite deformation index can be calculated, and then the geological deformation distribution map is obtained, which provides an important basis for evaluating the stability of the foundation.

[0036] Specifically, the geological activity threshold includes a plane displacement critical value and an elevation subsidence threshold. The plane displacement critical value refers to, in the horizontal direction, when the plane displacement of the measuring point exceeds this value, the measuring point is determined to be a plane dynamic unit, otherwise it is a plane static unit. The elevation subsidence threshold refers to, in the vertical direction, when the elevation displacement of the measuring point exceeds this value, the measuring point is determined to be a vertical dynamic unit, otherwise it is a vertical static unit. The plane movement coefficient and the elevation movement coefficient are parameters set according to the deformation of the measuring points, which are used to reflect the degree of deformation of the measuring points. For the plane dynamic unit, the plane movement coefficient is set to 1.2 times the plane displacement critical value; for the plane static unit, the plane movement coefficient is set to the ratio of the measured displacement to the critical value. For the vertical dynamic unit, the elevation movement coefficient is set to 1.5 times the elevation subsidence threshold; for the vertical static unit, the elevation movement coefficient is set to the logarithmic conversion value of the measured subsidence to the threshold. The setting of these parameters can more accurately reflect the deformation of the measuring points, and provide a basis for subsequent geological deformation analysis.

[0037] Preferably, when generating the geological deformation distribution map, a model based on the composite deformation index can be used. This model inputs the plane movement coefficient and the elevation movement coefficient of each measuring point, and calculates the composite deformation index through a specific algorithm. This algorithm takes into account the characteristic parameters of the silt soil layer, which are determined through field soil mechanics tests and can reflect the response characteristics of the soil layer to plane displacement and elevation subsidence. Specifically, the calculation method of the composite deformation index is to multiply the plane movement coefficient and the elevation movement coefficient by the corresponding characteristic parameters of the silt soil layer respectively, then add the results and take the square root. The obtained composite deformation index is subjected to geological time sequence filtering processing to eliminate noise and outliers, and finally the deformation distribution value of each measuring point is obtained to form the geological deformation distribution map. This processing process can effectively improve the accuracy and reliability of the geological deformation distribution map, and provide stronger support for the evaluation of the stability of the foundation.

[0038] In some embodiments, the dividing the measuring points into dynamic deformation units and static units through the geological activity threshold comprises: When the planar offset of the measuring point is greater than or equal to the planar displacement threshold, it is determined as a planar dynamic unit, and the planar movement coefficient is set to 1.2 times the planar displacement threshold; When the planar offset of the measuring point is less than the planar displacement threshold, it is determined as a planar static unit, and the planar movement coefficient is set to the ratio of the measured displacement to the threshold value; When the elevation offset of the measuring point is greater than or equal to the elevation subsidence threshold, it is determined as a vertical dynamic unit, and the elevation movement coefficient is set to 1.5 times the elevation subsidence threshold; When the elevation offset of the measuring point is less than the elevation subsidence threshold, it is determined as a vertical static unit, and the elevation movement coefficient is set to the logarithmic conversion value of the measured subsidence and the threshold value.

[0039] It should be noted that the present application adopts specific determination criteria and coefficient setting methods when dividing the measuring points into dynamic deformation units and static units through geological activity thresholds. This method can more accurately identify the deformation state of the measuring points, thereby providing more reliable data support for subsequent geological deformation analysis. Specifically, when the planar offset of the measuring point is greater than or equal to the planar displacement threshold, it is determined as a planar dynamic unit, and the planar movement coefficient is set to 1.2 times the planar displacement threshold; when the planar offset of the measuring point is less than the planar displacement threshold, it is determined as a planar static unit, and the planar movement coefficient is set to the ratio of the measured displacement to the threshold value. Similarly, when the elevation offset of the measuring point is greater than or equal to the elevation subsidence threshold, it is determined as a vertical dynamic unit, and the elevation movement coefficient is set to 1.5 times the elevation subsidence threshold; when the elevation offset of the measuring point is less than the elevation subsidence threshold, it is determined as a vertical static unit, and the elevation movement coefficient is set to the logarithmic conversion value of the measured subsidence and the threshold value. These determination criteria and coefficient settings can effectively distinguish the deformation state of the measuring points and provide a basis for subsequent composite deformation index calculation.

[0040] Specifically, the plane displacement threshold and the elevation subsidence threshold are key parameters for distinguishing the deformation state of a measuring point. The plane displacement threshold refers to that in the horizontal direction, when the plane offset of a measuring point exceeds this value, the measuring point is considered to be in a dynamic deformation state; otherwise, the measuring point is considered to be in a static deformation state. The elevation subsidence threshold refers to that in the vertical direction, when the elevation offset of a measuring point exceeds this value, the measuring point is considered to be in a dynamic deformation state; otherwise, the measuring point is considered to be in a static deformation state. The plane movement coefficient and the elevation movement coefficient are parameters for quantifying the deformation degree of a measuring point. For a plane dynamic unit, the plane movement coefficient is set to 1.2 times the plane displacement threshold, indicating that the deformation degree of the measuring point is relatively significant; for a plane static unit, the plane movement coefficient is set to the ratio of the measured displacement to the threshold, indicating that the deformation degree of the measuring point is relatively weak. For a vertical dynamic unit, the elevation movement coefficient is set to 1.5 times the elevation subsidence threshold, indicating that the subsidence degree of the measuring point is relatively significant; for a vertical static unit, the elevation movement coefficient is set to the logarithmic conversion value of the measured subsidence to the threshold, indicating that the subsidence degree of the measuring point is relatively weak. The setting of these parameters can more accurately reflect the deformation state of the measuring point, providing an important basis for subsequent geological deformation analysis.

[0041] Preferably, when determining the deformation state of a measuring point, a threshold-based determination model can be used. The model inputs the plane offset and the elevation offset of the measuring point, and determines the deformation state of the measuring point by comparing with the plane displacement threshold and the elevation subsidence threshold. The specific steps are as follows: first, obtain the plane offset and the elevation offset of the measuring point; then, compare the plane offset with the plane displacement threshold, if greater than or equal to the threshold, determine it as a plane dynamic unit, otherwise determine it as a plane static unit; then, compare the elevation offset with the elevation subsidence threshold, if greater than or equal to the threshold, determine it as a vertical dynamic unit, otherwise determine it as a vertical static unit. According to the determination result, set the corresponding plane movement coefficient and elevation movement coefficient. For example, for a plane dynamic unit, the plane movement coefficient is set to 1.2 times the plane displacement threshold; for a plane static unit, the plane movement coefficient is set to the ratio of the measured displacement to the threshold. For a vertical dynamic unit, the elevation movement coefficient is set to 1.5 times the elevation subsidence threshold; for a vertical static unit, the elevation movement coefficient is set to the logarithmic conversion value of the measured subsidence to the threshold. Through this specific determination model and parameter setting, the deformation state of the measuring point can be more accurately identified, providing more reliable data support for subsequent geological deformation analysis.

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

[0043] wherein: a composite deformation index for the measuring point a plane movement coefficient; an elevation movement coefficient; an elevation movement coefficient; a silt layer characteristic parameter (determined through on-site soil mechanics test, reflecting the response characteristics of the soil layer to plane displacement and elevation settlement).

[0044] It should be noted that the present application adopts a method of comprehensively considering the plane movement coefficient and the elevation movement coefficient when calculating the composite deformation index. This method combines the plane movement coefficient and the elevation movement coefficient through a specific calculation formula to generate a composite deformation index that can comprehensively reflect the deformation of the measuring point. The composite deformation index is a quantitative index used to evaluate the comprehensive deformation degree of the measuring point in the horizontal and vertical directions. In this way, the distribution of geological deformation can be more accurately reflected, providing important basis for subsequent evaluation of foundation stability.

[0045] Specifically, the calculation of the composite deformation index involves two key parameters: the plane movement coefficient and the elevation movement coefficient. The plane movement coefficient reflects the deformation degree of the measuring point in the horizontal direction, while the elevation movement coefficient reflects the deformation degree of the measuring point in the vertical direction. When calculating the composite deformation index, the silt layer characteristic parameters also need to be considered. These parameters are determined through on-site soil mechanics test and can reflect the response characteristics of the soil layer to plane displacement and elevation settlement. These parameters include α, β and γ, where α and β represent the weights of the plane movement coefficient and the elevation movement 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.

[0046] Preferably, in the calculation of the composite deformation index, a comprehensive model based on weights can be used. The input parameters of this model include the plane movement coefficient and the elevation movement coefficient of the measuring point, as well as the silt layer characteristic parameters α, β and γ determined through on-site test. The specific calculation steps are as follows: first, multiply the plane movement coefficient by the weight α and the elevation movement coefficient by the weight β; then, add the two results; finally, take the γth root of this sum to get the composite deformation index. This calculation method can comprehensively consider the deformation of the measuring point in the horizontal and vertical directions, and adjust the contribution of deformation in different directions through weight parameters, so as to obtain a more comprehensive and accurate deformation evaluation 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 map.

[0047] In some embodiments, the mapping generates a foundation stability coefficient matrix, comprising: The anti-interference processing of the deformation distribution value of each measuring point of the geological deformation distribution map is performed using the following formula:

[0048] Wherein: is the foundation stability coefficient; is the filtered deformation distribution value; is the tidal period influence factor (reflecting the periodic interference of tidal action on the foundation deformation); is the sediment density coefficient (a correction parameter related to the soil density); is the geological memory decay constant (characterizing the decay characteristics of soil deformation over time); is the local soil density of the measuring point (obtained by drilling sampling or geophysical prospecting method).

[0049] It should be noted that the present application adopts 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 of each measuring point of the geological deformation distribution map. This process takes into account various influencing factors, such as the tidal period influence factor, the sediment density coefficient, the geological memory decay constant, and the local soil density of the measuring point, to ensure that the calculation results accurately reflect the stability of the foundation. The foundation stability coefficient is an important quantitative index for evaluating the stability of the foundation under current geological conditions, providing key data support for the safe positioning of the ocean station structure.

[0050] Specifically, the calculation of the foundation stability coefficient involves multiple key parameters. The tidal period influence factor reflects the periodic interference of tidal action on the foundation deformation, and its value is usually adjusted according to historical tidal data and foundation response. The sediment density coefficient is related to the soil density and is used to correct the deformation influence 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 the foundation deformation. The local soil density of the measuring point is obtained by drilling sampling or geophysical prospecting method, reflecting the actual density of the soil at the measuring point. The reasonable 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 of each measuring point of the geological deformation distribution map are first filtered to obtain the filtered deformation distribution values, and then combined with the above parameters, the foundation stability coefficient is obtained through specific calculation formula.

[0051] Preferably, when generating the foundation stability coefficient matrix, a model that comprehensively considers multiple influencing factors can be used. The input parameters of the model include the filtered deformation distribution value, the tidal cycle influence factor, the sediment density coefficient, the geological memory decay constant, and the local soil density of the measuring point. The specific calculation steps are as follows: First, obtain the deformation distribution value of each measuring point in the geological deformation distribution map, and filter it to eliminate the influence of noise and outliers. Then, divide the filtered deformation distribution value by a comprehensive correction factor, which is composed of the tidal cycle influence factor, the sediment density coefficient, the geological memory decay constant, and the local soil density of the measuring point. In this way, a quantitative indicator that can accurately reflect the stability of the foundation can be obtained - the foundation stability coefficient. This calculation method can effectively consider the influence of various geological and environmental factors on the stability of the foundation, and provide more reliable data support for the safe positioning of the ocean station structure.

[0052] In some embodiments, the following equation is used to perform tidal compensation correction of the plane offset:

[0053] in: The measurement point of the nth monitoring cycle The plane offset of is the silt rheology compensation coefficient (the value range is 0~1, reflecting the weight of the rheological characteristics of silt soil on displacement correction); is the measured plane displacement value in the nth period; is the duration of the nth monitoring cycle; For measuring points Tidal velocity field data (obtained through tidal models or on-site current meters); is the plane offset of the n-1th cycle (historical correction value).

[0054] It should be noted that the present invention utilizes a correction method that comprehensively considers multiple factors when performing tidal compensation for planar offsets. This method uses a specific algorithm to adjust the measured planar displacement values ​​to eliminate the influence of tidal velocity field data on the measurement results. The correction process also considers the silt rheology compensation coefficient, a key parameter that reflects the weight of the rheological properties of silty soil on displacement correction. This correction method can produce more accurate planar offsets, thereby improving the accuracy of ocean station structure measurement and positioning.

[0055] Specifically, the tidal compensation correction of the plane offset involves several key parameters. The measured plane displacement value refers to the displacement data of the ocean station structure in the horizontal direction directly measured by the monitoring equipment. The tidal current velocity field data refers to the tidal current velocity information obtained by the tidal model or the field current meter during the monitoring period, which reflects the influence degree of the tide on the displacement of the ocean station structure. The silt rheological compensation coefficient is a value between 0 and 1, which is used to adjust the influence degree of the tidal current on the displacement, and its value is determined according to the rheological properties of the silt soil. In the correction process, first, the tidal influence amount of the current period is calculated, then the measured displacement value is adjusted according to the silt rheological compensation coefficient, and finally, the corrected plane offset is obtained by combining the historical correction value of the last period. This correction method can effectively eliminate the interference of the tide on the measurement results and improve the accuracy of the plane offset.

[0056] Preferably, when performing tidal compensation correction of the plane offset, a correction model based on time series can be used. The input parameters of the model include the measured plane displacement value of the current period, the tidal current velocity field data, the historical correction value of the last period, and the silt rheological compensation coefficient. The specific correction steps are as follows: first, calculate the tidal influence amount of the current period, which is the product of the tidal current velocity field data and the monitoring period length. Then, adjust the measured plane displacement value according to the silt rheological compensation coefficient to obtain the tidal corrected displacement value. Finally, the tidal corrected displacement value and the historical correction value of the last period are weighted and averaged to obtain the corrected plane offset. This correction model can effectively consider the periodic influence of the tide and the rheological properties of the silt soil, so as to obtain more accurate plane offset and provide more reliable data support for the safe positioning of the ocean station structure.

[0057] In some embodiments, the deposition compensation correction of the elevation offset is performed using the following formula:

[0058] wherein: is the elevation offset of the monitoring point in the nth monitoring period; is the elevation offset of the monitoring point in the nth monitoring period; is the deposition consolidation degree coefficient (value range 0~1, reflecting the weight of the deposition consolidation degree on the settlement correction); is the measured elevation settlement value in the nth period; is the compression modulus of the ith type of sediment (reflecting the physical parameter of the soil layer compression characteristics); is the thickness change amount of the corresponding ith type of sediment layer (obtained by stratum monitoring); is the elevation offset of the nth-1 period (historical correction value).

[0059] It should be noted that the present application adopts a correction method considering the consolidation degree and compression characteristics of sediments when performing the deposition compensation correction of the elevation offset. This method adjusts the measured elevation settlement value through a specific algorithm to eliminate the influence of sediment consolidation and compression on the measurement results. The sediment consolidation degree coefficient is also considered in the correction process, which is an important parameter reflecting the influence of sediment consolidation degree on the settlement correction weight. Through this correction method, more accurate elevation offset can be obtained, thereby improving the positioning accuracy of the marine station structure measurement.

[0060] Specifically, the deposition compensation correction of the elevation offset involves several key parameters. The measured elevation settlement value refers to the vertical settlement data of the marine station structure measured directly by monitoring equipment. The sediment consolidation degree coefficient is a value between 0 and 1, used to adjust the influence of sediment consolidation on settlement, and its value is determined according to the consolidation characteristics of the sediment. The compression modulus of the sediment refers to a physical parameter reflecting the compression characteristics of the sediment, which is related to the type and nature of the sediment. In the correction process, first, the compression amount of the sediment layer is calculated, which is the product of the compression modulus of the sediment and the thickness change amount of the sediment layer. Then, the measured settlement value is adjusted according to the sediment consolidation degree coefficient, and finally, the historical correction value of the last period is combined to obtain the corrected elevation offset. This correction method can effectively eliminate the interference of sediment consolidation and compression on the measurement results, and improve the accuracy of the elevation offset.

[0061] Preferably, when performing the deposition compensation correction of the elevation offset, a correction model based on the characteristics of the sediment layer can be used. The input parameters of this model include the measured elevation settlement value of the current period, the compression modulus of the sediment, the thickness change amount of the sediment layer, the historical correction value of the last period, and the sediment consolidation degree coefficient. The specific correction steps are as follows: first, calculate the compression amount of the sediment layer, which is the product of the compression modulus of the sediment and the thickness change amount of the sediment layer. Then, adjust the measured settlement value according to the sediment consolidation degree coefficient to obtain the settlement value after deposition correction. Finally, the settlement value after deposition correction and the historical correction value of the last period are weighted and averaged to obtain the corrected elevation offset. This correction model can effectively consider the consolidation and compression characteristics of the sediment, thereby obtaining more accurate elevation offset and providing more reliable data support for the safe positioning of the marine station structure.

[0062] In some embodiments, the formula for generating the structure safety positioning coordinates by fusion is:

[0063] wherein: represents Hadamard product (multiplication of corresponding elements of matrices); represents a weighted fusion operator considering the creep characteristics of the soil body (weight distribution is determined according to the soil creep model). for the fused structure safety positioning coordinates; for the current monitoring period measured coordinate data; for the optimal reference value in the historical safety coordinate database (based on the coordinate reference of the historical stable state).

[0064] It should be noted that the present application adopts a method of comprehensively considering the current monitoring period measured coordinate data and the historical safety coordinate data when generating the fused structure safety positioning coordinates. This method generates the structure safety positioning coordinates by combining the foundation stability coefficient matrix with the historical positioning reference through a specific fusion algorithm. The structure safety positioning coordinates are an important output result, which can accurately reflect the safe position of the ocean station structure under the current geological conditions, and provide key data support for the operation and maintenance of the ocean station.

[0065] Specifically, the generation of the fused structure safety positioning coordinates involves several key parameters. The foundation stability coefficient matrix is calculated through the previous steps, which reflects the stability of the foundation of each measuring point. The current monitoring period measured coordinate data refers to the coordinate data of the ocean station structure obtained by measuring equipment in the current monitoring period. The optimal reference value in the historical safety coordinate database refers to the coordinate reference based on the historical stable state, which is usually selected from the stable coordinate data obtained from long-term monitoring. In the fusion process, the foundation stability coefficient matrix is first processed, and then it is fused with the current monitoring period measured coordinate data and the historical safety coordinate data. The Hadamard product and the weighted fusion operator considering the creep characteristics of the soil are involved in the fusion algorithm, which ensures that the fusion result can comprehensively consider the current and historical data, and also considers the creep characteristics of the soil.

[0066] Preferably, when generating the fused structure safety positioning coordinates, a weight-based fusion model can be used. The input parameters of this model include the foundation stability coefficient matrix, the current monitoring period measured coordinate data and the historical safety coordinate data. The specific fusion steps are as follows: first, process the foundation stability coefficient matrix to obtain the stability coefficient of each measuring point. Then, perform Hadamard product operation on the stability coefficient of each measuring point and the current monitoring period measured coordinate data, i.e. multiply the corresponding elements. Next, perform weighted fusion on the obtained results and the historical safety coordinate data, and the weight is determined according to the creep characteristics of the soil. This fusion model can effectively consider the current monitoring data and the historical data, and also consider the creep characteristics of the soil, so as to obtain more accurate structure safety positioning coordinates, and provide more reliable data support for the safety evaluation and positioning of the ocean station structure.

[0067] In some embodiments, the generation of the planar displacement map and the elevation displacement map comprises: The four-dimensional sonar interferometric measurement technology is adopted to obtain seabed topography change data in continuous monitoring periods; The strain field data of the structure foundation are collected by a distributed optical fiber sensing system, and the strain field data are calculated into plane displacement components and elevation displacement components; Adaptive Kalman filtering is used for noise suppression for the measurement area disturbed by waves to generate an anti-interference displacement atlas.

[0068] It should be noted that the present application adopts an implementation mode of integrating multiple monitoring technologies and data processing methods when generating the plane displacement atlas and the elevation displacement atlas. This method obtains seabed topography change data in continuous monitoring periods by the four-dimensional sonar interferometric measurement technology, and collects the strain field data of the structure foundation by the distributed optical fiber sensing system. After processing, these data can generate an anti-interference displacement atlas, thereby providing accurate basic data for subsequent geological deformation analysis. This method is particularly suitable for complex geological environments such as muddy coastal zones, and can effectively improve the accuracy and reliability of measurement.

[0069] Specifically, the four-dimensional sonar interferometric measurement technology is a seabed topography monitoring means, which can provide high-precision seabed topography change data. These data reflect the topography change of the ocean station foundation at different time points, and provide an important basis for generating the plane displacement atlas. The distributed optical fiber sensing system is used to collect the strain field data of the structure foundation, which can reflect the stress and strain conditions inside the structure and is an important basis for evaluating the safety of the structure. The process of calculating the strain field data into plane displacement components and elevation displacement components involves an algorithm for converting strain data into displacement data. For the measurement area disturbed by waves, adaptive Kalman filtering is used for noise suppression, which is an algorithm that can dynamically adjust the filtering parameters according to the characteristics of the measurement data, and can effectively reduce the interference of environmental factors such as waves on the measurement results to generate an anti-interference displacement atlas.

[0070] Preferably, in generating the plane displacement map and the elevation displacement map, a comprehensive data processing model can be adopted. The input parameters of the model include the seafloor topography change data obtained by the four-dimensional interferometric bathymetry technology and the structural foundation strain field data collected by the distributed optical fiber sensing system. The specific processing steps are as follows: first, the four-dimensional interferometric bathymetry data is pre-processed to extract the key features of the seafloor topography change. Then, the strain field data is calculated to obtain the plane displacement component and the elevation displacement component. Next, for the measurement area disturbed by waves, an adaptive Kalman filtering algorithm is applied for noise suppression. The adaptive Kalman filtering 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 generated into an anti-interference displacement map, including the plane displacement map and the elevation displacement map. This comprehensive data processing model can effectively process various 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 evaluation.

[0071] The above-mentioned various embodiments of the present application have the following beneficial effects: 1、The present application can comprehensively reflect the structural displacement under complex geological environment by acquiring multi-source monitoring data of the muddy coastal zone and extracting the deformation feature points of the ocean station structure. The plane displacement is corrected by tidal compensation, and the elevation displacement is corrected by sedimentation rate compensation, effectively eliminating the interference of tidal action and sedimentation process on the measurement results, and significantly improving the accuracy of displacement monitoring. This method solves the technical problems of large environmental factor influence and measurement error accumulation of traditional measurement technology in the muddy coastal zone, and provides a reliable data foundation for subsequent structural safety evaluation.

[0072] 2、The present application innovatively divides the measuring points into dynamic deformation units and static units by the geological activity threshold, and calculates the composite deformation index based on the plane movement coefficient and the elevation movement coefficient, realizing the fine evaluation of the foundation stability. Through the geological time sequence filtering processing of the composite deformation index, the interference of short-term geological activity is eliminated, and the long-term deformation trend of the foundation is accurately reflected. This method solves the problem that traditional monitoring methods cannot distinguish between dynamic deformation and static stable regions, and provides a scientific basis for engineering safety warning, especially suitable for muddy soil layers with significant rheological properties.

[0073] 3、The application intelligently fuses the ground stability coefficient matrix of each monitoring period with the historical positioning reference, adopts a weighted fusion operator considering the creep characteristics of the soil body, and dynamically generates the structural safety positioning coordinates. This method effectively solves the problem of positioning reference drift caused by factors such as soil consolidation and creep in the long-term monitoring process, and ensures the long-term stability of the ocean station structure under complex geological conditions. Compared with the traditional static positioning method, the application can adapt to the changes in the geological environment, and provides reliable technical support for the safe operation and maintenance of marine engineering in the muddy coastal area.

[0074] Further, the storage medium of the embodiments of the present application stores program instructions capable of implementing all the methods described above, wherein the program instructions can be stored in the storage medium in the form of a software product, including a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) or a processor to execute all or part of the steps of the methods described in the embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, and various media capable of storing program codes, or a terminal device such as a computer, a server, a mobile phone, a tablet, etc.

[0075] The above description is only some of the preferred embodiments of the present application and an explanation of the technical principles applied. Those skilled in the art should understand that the scope of the application involved in the embodiments of the present application is not limited to the technical solutions formed by the specific combinations of the above technical features, and should also cover other technical solutions formed by any combination of the above technical features or equivalent features without departing from the above inventive concept. For example, the above features are replaced with technical features having similar functions disclosed in the embodiments of the present application (but not limited to) to form technical solutions.

Claims

1. A method for measuring and positioning ocean station structures in muddy coastal zones, characterized in that: The steps include: Obtain multi-source monitoring data for muddy coastal areas; Extracting deformation feature points of the ocean station structure from the multi-source monitoring data and generating a plane displacement map and an elevation displacement map of the feature points of each monitoring period relative to a reference period; Based on the plane displacement map and the elevation displacement map, a spatial offset calculation is performed to generate a geological deformation distribution map of each monitoring period and a foundation stability coefficient matrix is ​​generated by mapping; The foundation stability coefficient matrix of each monitoring period is integrated with the historical positioning benchmark to generate the structural safety positioning coordinates.

2. The method according to claim 1, characterized in that The spatial offset calculation based on the plane displacement map and the elevation displacement map includes: The displacement value of each measuring point in the plane displacement diagram of each monitoring period is corrected by tidal compensation with the plane offset of the corresponding measuring point in the reference period to obtain the plane offset of each measuring point in the current period; 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.

3. The method according to claim 2, characterized in that The generating of the geological deformation distribution map of 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 the 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 composite deformation index is subjected to geological time series filtering to obtain deformation distribution values ​​of each measuring point, thereby forming the geological deformation distribution map.

4. The method according to claim 3, characterized in that The method of dividing the measuring points into dynamic deformation units and static units by using the geological activity threshold includes: When the plane offset of the measuring point is greater than or equal to the plane displacement critical value, it is determined to be a plane dynamic unit, and its plane movement coefficient is set to 1.2 times the plane displacement critical value; When the plane offset of the measuring point is less than the critical value of the plane displacement, it is determined to be a plane static unit, and its plane movement coefficient is set as the ratio of the measured displacement to the critical value; 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 conversion value of the measured settlement and the threshold.

5. The method according to claim 3 or 4, characterized in that The composite deformation index is calculated based on the plane movement coefficient and elevation movement coefficient of each measuring point using the following formula: in: For measuring points Composite deformation index; is the plane movement coefficient; is the elevation movement coefficient; It is the characteristic parameter of the silty soil layer (determined through on-site soil mechanics tests, reflecting the response characteristics of the soil layer to plane displacement and elevation settlement).

6. The method according to claim 5, characterized in that The mapping generates a foundation stability coefficient matrix, including: The deformation distribution value of each measuring point in the geological deformation distribution map is processed with the following formula for anti-interference: in: is the foundation stability coefficient; is the deformation distribution value after filtering; is the tidal cycle influencing factor; is the sediment density coefficient; is the geological memory decay constant; is the local soil density at the measuring point.

7. The method according to claim 2, characterized in that Use the following formula to make tidal compensation correction for plane offset: in: The measurement point of the nth monitoring cycle The plane offset of is the silt rheological compensation coefficient; is the measured plane displacement value in the nth period; is the duration of the nth monitoring cycle; For measuring points Tidal velocity field data; is the plane offset of the n-1th period.

8. The method according to claim 7, characterized in that Use the following formula to make sedimentation compensation correction for elevation offset: in: The measurement point of the nth monitoring cycle The elevation offset of is the sediment consolidation coefficient; is the measured elevation settlement value in the nth period; is the compression modulus of the i-th type of sediment; is the thickness change of the corresponding sedimentary layer of type i; is the elevation offset of the n-1th cycle.

9. The method according to claim 1, characterized in that The formula for the fusion-generated structure safety positioning coordinates is: in: represents the Hadamard product; represents the weighted fusion operator considering the creep characteristics of soil; Securely locate the coordinates of the fused structure; is the measured coordinate data of the current monitoring period; It is the optimal reference value in the historical safety coordinate database.

10. The method according to claim 1, characterized in that The generating of the plane displacement map and the elevation displacement map includes: Use four-dimensional sonar interferometry to obtain seabed topography change data over a continuous monitoring period; Collecting structural foundation strain field data through a distributed optical fiber sensing system; and resolving the strain field data into a plane displacement component and an elevation displacement component; Adaptive Kalman filtering is used to suppress noise in the measurement area disturbed by waves to generate an anti-interference displacement map.

Images