A method, device, equipment and medium for determining optimal earth radial parameters

Abstract

The present application provides a method, device, equipment and medium for determining the optimal earth radial parameters. Among them, the original observation data of a number of seafloor stations in the epicentral area of submarine earthquakes are obtained; the seafloor station coordinate position solution data and the observed values of the displacement time series are determined based on the original observation data of each seafloor station; a Burgers viscoelastic earth model is constructed, and the reset dislocation Love number and the complex Green's function are determined based on the Burgers viscoelastic earth model; based on the reset dislocation Love number and the complex Green's function, the calculated values of the displacement time series are determined according to the observed values of the displacement time series; the optimal transient viscosity parameter and the optimal steady-state viscosity parameter are determined based on the observed values and the calculated values of the displacement time series, and the transient viscosity parameter and the optimal steady-state viscosity parameter are jointly composed of the optimal earth radial parameters. By adopting the above scheme, the earth radial parameters can be accurately determined to ensure the correct understanding and accurate prediction of the earthquake process.

Description

Technical Field

[0001] The present invention relates to the technical field of GNSS marine applications, and more particularly, to a method, apparatus, device and medium for determining optimal earth radial parameters. Background Art

[0002] Currently, the research on the viscoelastic structure of the earth model mainly focuses on the radial stratification of the land, because there is usually a relatively dense observation network in the land area, which can provide data with a high time sampling rate. These data are crucial for constraining the inversion process. GNSS stations are widely distributed on land and can provide continuous and high-precision position information, enabling researchers to monitor the movement and deformation of the crust more accurately. This rich data resource provides a solid foundation for the research on the earth model in the land area, enabling the radial stratification model of the land to better reflect the viscoelastic characteristics of the crust.

[0003] However, when the research perspective turns to the seabed, the situation is quite different. The complexity of the seabed environment far exceeds that of the land. Its structure not only has stratification characteristics in the vertical direction but also has significant differences in the horizontal direction. This lateral difference is mainly due to factors such as the diversity of the seabed topography, the variation of the lithosphere thickness, and the interaction between different plates. For example, at the plate boundary in the subduction zone, the crustal structure and rock properties change drastically, resulting in obvious spatial inhomogeneity of the viscoelastic characteristics. In addition, seabed volcanic activities, sediment distribution, and marine geological structures also have important impacts on the seabed structure, making the viscoelastic characteristics of the seabed more complex and diverse.

[0004] Due to the existence of these complex factors, it is no longer possible to accurately describe the viscoelastic characteristics of the seabed structure using only a simple radial stratification model. The traditional radial stratification method ignores the structural changes in the horizontal direction and cannot fully consider the influence of the seabed topography and geological structure on the viscoelastic characteristics. Therefore, when conducting seabed earthquake research, if this simplified model continues to be used, it will lead to deviations in the simulation and prediction of the earthquake process. Specifically, key processes such as the earthquake source mechanism, the propagation path of seismic waves, and the crustal deformation caused by the earthquake may not be correctly explained and predicted due to the inaccuracy of the model. And the earth radial parameters, as the core components of the earth model, are of great significance for the research on the earthquake process. These parameters not only affect the earthquake source mechanism and the propagation characteristics of seismic waves but are also closely related to the crustal deformation and stress distribution caused by the earthquake. If the earth radial parameters cannot be accurately determined, it will directly affect the understanding and prediction of the earthquake process. Summary of the Invention

[0005] In view of this, the purpose of the present invention is to provide a method, device, equipment and medium for determining the optimal earth radial parameters, so as to accurately determine the earth radial parameters and ensure the correct understanding and accurate prediction of the earthquake process.

[0006] In a first aspect, an embodiment of the present application provides a method for determining the optimal earth radial parameters, and the method includes:

[0007] Obtain the original observation data of several undersea stations in the epicenter area of the undersea earthquake;

[0008] Based on the original observation data of each undersea station, determine the undersea station coordinate position solution data and the observed values of the displacement time series;

[0009] Construct a Burgers viscoelastic earth model, and based on the Burgers viscoelastic earth model, determine the reset dislocation Love number and the complex Green's function, where the Burgers viscoelastic earth model includes transient viscosity parameters and steady-state viscosity parameters, and the transient viscosity parameters and the steady-state viscosity parameters together constitute the earth radial parameters;

[0010] Based on the reset dislocation Love number and the complex Green's function, determine the calculated values of the displacement time series according to the observed values of the displacement time series;

[0011] Based on the observed values and calculated values of the displacement time series, determine the optimal transient viscosity parameter and the optimal steady-state viscosity parameter, and jointly form the optimal earth radial parameters by the transient viscosity parameter and the optimal steady-state viscosity parameter.

[0012] Optionally, the obtaining the original observation data of several undersea stations in the epicenter area of the undersea earthquake includes:

[0013] Determine the actual undersea station information of the research area according to the location information of the research area;

[0014] Screen the undersea station information in the epicenter area from the actual undersea station information of the research area;

[0015] Extract the original observation data of each undersea station from the undersea station information in the epicenter area, where the original observation data includes GNSS data, attitude data, and acoustic ranging data of each undersea station in the epicenter area of the undersea earthquake.

[0016] Optionally, the determining the undersea station coordinate position solution data and the observed values of the displacement time series based on the original observation data of each undersea station includes:

[0017] According to the original observation data of each undersea station, use the sound speed spatio-temporal variation error compensation positioning model to solve the undersea station coordinate position to obtain the undersea station coordinate position solution data;

[0018] Screen and extract the displacement time series of each subsea station according to the subsea station coordinate position solution data, and determine the observed values of the displacement time series.

[0019] Optionally, the determining the complex dislocation Love number and the complex Green's function based on the Burgers viscoelastic Earth model includes:

[0020] Introduce a complex variable parameter s into the stress-strain constitutive relation of the Burgers viscoelastic Earth model, and perform Laplace transform to obtain the constitutive relation in the complex domain;

[0021] Solve the target equations in the complex domain according to the constitutive relation in the complex domain to obtain the complex dislocation Love number and the complex Green's function.

[0022] Optionally, the target equations include:

[0023] The spherical center initial condition function, the inner core-outer core boundary condition function, the outer core-mantle boundary condition function, the free surface boundary condition function on the Earth's surface, and the source function.

[0024] Optionally, the determining the calculated values of the displacement time series based on the complex dislocation Love number and the complex Green's function according to the observed values of the displacement time series specifically includes:

[0025] Use the Riemann-Mellin inversion formula to perform inverse transform on the complex dislocation Love number and the complex Green's function to obtain the dislocation Love number and the Green's function in the time domain;

[0026] Use the subsea station coordinate position solution data to solve the dislocation Love number and the Green's function to obtain the calculated values of the displacement time series.

[0027] Optionally, the determining the optimal transient viscosity parameter and the optimal steady-state viscosity parameter based on the observed values and the calculated values of the displacement time series includes:

[0028] Constrain and adjust the observed values and the calculated values of the displacement time series to determine the optimal calculated values of the displacement time series that minimize the root mean square error between the observed values and the calculated values of the displacement time series;

[0029] Respectively determine the transient viscosity parameter and the steady-state viscosity parameter corresponding to the optimal calculated values of the displacement time series in the Burgers viscoelastic Earth model as the optimal transient viscosity parameter and the optimal steady-state viscosity parameter.

[0030] In a second aspect, an embodiment of the present application provides an optimal Earth radial parameter determination device, and the device includes:

[0031] An observation data acquisition module for acquiring original observation data of a number of undersea stations in the epicentral area of an undersea earthquake;

[0032] A displacement time series observation value determination module for determining undersea station coordinate position calculation data and observation values of the displacement time series based on the original observation data of each undersea station;

[0033] A viscoelastic earth model construction module for constructing a Burgers viscoelastic earth model, and determining a reset dislocation Love number and a complex Green's function based on the Burgers viscoelastic earth model, wherein the Burgers viscoelastic earth model includes a transient viscosity parameter and a steady-state viscosity parameter, and the transient viscosity parameter and the steady-state viscosity parameter together constitute earth radial parameters;

[0034] A displacement time series calculated value determination module for determining calculated values of the displacement time series based on the reset dislocation Love number and the complex Green's function according to the observation values of the displacement time series;

[0035] An optimal earth radial parameter determination module for determining an optimal transient viscosity parameter and an optimal steady-state viscosity parameter based on the observation values and calculated values of the displacement time series, and combining the transient viscosity parameter and the optimal steady-state viscosity parameter to form optimal earth radial parameters.

[0036] Optionally, the acquiring of the original observation data of a number of undersea stations in the epicentral area of the undersea earthquake includes:

[0037] Determining actual undersea station information of the research area according to the position information of the research area;

[0038] Screening out undersea station information in the epicentral area from the actual undersea station information of the research area;

[0039] Extracting the original observation data of each undersea station from the undersea station information in the epicentral area, wherein the original observation data includes GNSS data, attitude data, and acoustic ranging data of each undersea station in the epicentral area of the undersea earthquake.

[0040] Optionally, the determining of the undersea station coordinate position calculation data and the observation values of the displacement time series based on the original observation data of each undersea station includes:

[0041] Calculating the undersea station coordinate position to obtain the undersea station coordinate position calculation data by using a sound speed spatio-temporal variation error compensation positioning model according to the original observation data of each undersea station;

[0042] Screen and extract the displacement time series of each seabed station according to the seabed station coordinate position solution data, and determine the observed values of the displacement time series.

[0043] Optionally, determining the complex dislocation Love number and the complex Green's function based on the Burgers viscoelastic Earth model includes:

[0044] Introduce a complex variable parameter s into the stress-strain constitutive relation of the Burgers viscoelastic Earth model, and perform Laplace transform to obtain the constitutive relation in the complex domain;

[0045] Solve the target equations in the complex domain according to the constitutive relation in the complex domain to obtain the complex dislocation Love number and the complex Green's function.

[0046] Optionally, the target equations include:

[0047] The spherical center initial condition function, the inner core-outer core boundary condition function, the outer core-mantle boundary condition function, the free surface boundary condition function on the Earth's surface, and the source function.

[0048] Optionally, based on the complex dislocation Love number and the complex Green's function, determining the calculated values of the displacement time series according to the observed values of the displacement time series includes:

[0049] Use the Riemann-Mellin inversion formula to perform inverse transformation on the complex dislocation Love number and the complex Green's function to obtain the dislocation Love number and the Green's function in the time domain;

[0050] Use the seabed station coordinate position solution data to solve the dislocation Love number and the Green's function to obtain the calculated values of the displacement time series.

[0051] Optionally, determining the optimal transient viscosity parameter and the optimal steady-state viscosity parameter based on the observed values and the calculated values of the displacement time series includes:

[0052] Constrain and adjust the observed values and the calculated values of the displacement time series to determine the optimal calculated values of the displacement time series that minimize the root mean square error between the observed values and the calculated values of the displacement time series;

[0053] Respectively determine the transient viscosity parameter and the steady-state viscosity parameter corresponding to the optimal calculated values of the displacement time series in the Burgers viscoelastic Earth model as the optimal transient viscosity parameter and the optimal steady-state viscosity parameter.

[0054] In a third aspect, an embodiment of the present application provides a computer device, including: a processor, a memory, and a bus. The memory stores machine-readable instructions executable by the processor. When the computer device runs, the processor communicates with the memory through the bus. When the machine-readable instructions are executed by the processor, the steps of the optimal earth radial parameter determination method described in any optional implementation manner of the first aspect are executed.

[0055] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium. A computer program is stored on the computer-readable storage medium. When the computer program is run by a processor, the steps of the optimal earth radial parameter determination method described in any optional implementation manner of the first aspect are executed.

[0056] The technical solutions provided by the present application include, but are not limited to, the following beneficial effects:

[0057] The present application first obtains the original observation data of the seafloor stations in the epicenter area of the submarine earthquake, which can directly supplement the deficiencies of land observation data and provide more comprehensive and accurate earthquake monitoring information. These data include GNSS data, attitude data, acoustic ranging data, etc., which can reflect the precise position and displacement changes of the seafloor stations. By obtaining these data, the source mechanism and earthquake process of the submarine earthquake can be understood more accurately, providing basic data support for subsequent model construction and parameter inversion.

[0058] Then, based on the original observation data of each seafloor station, the coordinate position solution data of the seafloor station and the observed values of the displacement time series are determined, which can accurately determine the position changes of the seafloor stations and reflect the displacement changes of the seafloor stations during the earthquake process. These observed values are important input data for subsequent model construction and parameter inversion, and can help describe the earthquake process and crustal deformation more accurately.

[0059] Then, a Burgers viscoelastic earth model is constructed. Based on the Burgers viscoelastic earth model, the reset dislocation Love number and the complex Green's function are determined, which can more accurately describe the viscoelastic characteristics of the seafloor structure. This model includes transient viscosity parameters and steady-state viscosity parameters, and can simultaneously consider the transient and steady-state deformations during the earthquake process. By calculating the reset dislocation Love number and the complex Green's function, the theoretical solutions of the displacement field caused by the earthquake can be obtained. These theoretical solutions provide a basis for subsequent displacement time series calculations and can help simulate the earthquake process and crustal deformation more accurately.

[0060] Next, based on the dislocation Love numbers and complex Green's functions, the calculated values of the displacement time series are determined from the observed values of the displacement time series, enabling a comparison between the theoretical solution of the model and the actual observed data to verify the accuracy and reliability of the model. By comparing the observed values and the calculated values, the performance of the model can be evaluated, providing a basis for subsequent parameter inversion.

[0061] Finally, based on the observed values and calculated values of the displacement time series, the optimal transient viscosity parameter and the optimal steady-state viscosity parameter are determined. Combining the transient viscosity parameter and the optimal steady-state viscosity parameter to form the best Earth radial parameter can minimize the error between the calculated value and the observed value of the model, improving the accuracy and reliability of the model. Combining the optimal transient viscosity parameter and the steady-state viscosity parameter to form the best Earth radial parameter can more accurately describe the viscoelastic properties of the seabed structure, providing important parameter support for the study of seismic processes. These parameters help to more accurately simulate seismic processes, predict seismic impacts, and formulate more effective disaster prevention and mitigation measures.

[0062] To make the above objects, features, and advantages of the present invention more obvious and understandable, the following specifically enumerates preferred embodiments and, in conjunction with the accompanying drawings, provides the following detailed description. Brief Description of the Drawings

[0063] To more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for use in the embodiments. It should be understood that the following drawings only show some embodiments of the present invention and should not be regarded as limiting the scope. For those of ordinary skill in the art, without creative efforts, other related drawings can also be obtained based on these drawings.

[0064] Figure 1 Shows the flowchart of a method for determining the best Earth radial parameter provided in Embodiment 1 of the present invention;

[0065] Figure 2 Shows the flowchart of a method for obtaining original observed data provided in Embodiment 1 of the present invention;

[0066] Figure 3 Shows a schematic diagram of a partial position of a seabed station provided in Embodiment 1 of the present invention;

[0067] Figure 4 Shows the flowchart of a method for determining the observed values of the displacement time series provided in Embodiment 1 of the present invention;

[0068] Figure 5 Shows the flowchart of a method for constructing a viscoelastic Earth model provided in Embodiment 1 of the present invention;

[0069] Figure 6Shows a schematic diagram of a viscoelastic earth model provided by Embodiment 1 of the present invention;

[0070] Figure 7 Shows a flowchart of calculated values of a displacement time series provided by Embodiment 1 of the present invention;

[0071] Figure 8 Shows a flowchart of a method for determining optimal parameters provided by Embodiment 1 of the present invention;

[0072] Figure 9 Shows a schematic structural diagram of an apparatus for determining optimal earth radial parameters provided by Embodiment 2 of the present invention;

[0073] Figure 10 Shows a schematic structural diagram of a computer device provided by Embodiment 3 of the present invention. Detailed implementation manners

[0074] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the described embodiments are only some of the embodiments of the present invention, rather than all of the embodiments. Usually, the components of the embodiments of the present invention described and illustrated herein can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely represents selected embodiments of the present invention. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0075] Embodiment 1

[0076] For the convenience of understanding the present application, the following Figure 1 describes Embodiment 1 of the present application in detail with reference to the content described in the flowchart of a method for determining optimal earth radial parameters provided by Embodiment 1 of the present invention shown.

[0077] See Figure 1 as shown, Figure 1 shows a flowchart of a method for determining optimal earth radial parameters provided by Embodiment 1 of the present invention, wherein the method includes steps S101 to S105:

[0078] S101: Obtain the original observation data of several undersea stations in the epicentral area of the undersea earthquake.

[0079] Specifically, based on the location information of the research area, the actual seabed station information is determined, including the coordinate position, the plate it is located in, and the quantity of the seabed stations. The seabed station information in the epicenter area is screened out from the actual seabed station information of the research area, and the original GNSS-A (GNSS-Acoustic, Global Navigation Satellite System - Acoustic Observation Model) observation data of each seabed station are extracted, such as GNSS (Global Navigation Satellite System) data, attitude data, and acoustic ranging data.

[0080] S102: Determine the calculated data of the seabed station coordinate position and the observed values of the displacement time series based on the original observation data of each seabed station.

[0081] Specifically, based on the original observation data of each seabed station, the coordinate position of the seabed station is calculated using the sound speed spatio-temporal variation error compensation positioning model. This process can compensate for the variation of the sound speed under different time and space conditions, thereby improving the positioning accuracy. Through this calculation process, the calculated data of the seabed station coordinate position can be obtained. Then, based on the calculated data of the seabed station coordinate position, the displacement time series of each seabed station are screened and extracted, and the observed values of the displacement time series are determined. These observed values reflect the displacement changes of the seabed station after the earthquake and provide important input data for subsequent model construction and parameter inversion.

[0082] S103: Construct a Burgers viscoelastic Earth model, and determine the reduced dislocation Love number and the complex Green's function based on the Burgers viscoelastic Earth model. Among them, the Burgers viscoelastic Earth model includes transient viscosity parameters and steady-state viscosity parameters, and the transient viscosity parameters and the steady-state viscosity parameters together constitute the Earth's radial parameters.

[0083] Specifically, construct a Burgers viscoelastic Earth model. This model combines a Kelvin body and a Maxwell body and can describe the viscoelastic behavior of the Earth's medium. The model includes transient viscosity parameters and steady-state viscosity parameters, and these two parameters together constitute the Earth's radial parameters, which are the targets of subsequent parameter inversion. Next, a complex variable parameter s is introduced into the stress-strain constitutive relationship of the model, and a Laplace transform is performed to obtain the constitutive relationship in the complex domain. This process transforms the viscoelastic problem into an equivalent elastic problem in the complex domain, simplifying the calculation process. Finally, according to the constitutive relationship in the complex domain, the target equations in the complex domain are solved to obtain the reduced dislocation Love number and the complex Green's function. These functions describe the displacement field caused by the earthquake and provide a basis for subsequent calculation of the displacement time series.

[0084] S104: Based on the reduced dislocation Love number and the complex Green's function, determine the calculated values of the displacement time series according to the observed values of the displacement time series.

[0085] Specifically, by using the Riemann-Mellin inversion formula, the dislocation Love numbers and complex Green's functions are inversely transformed to obtain the dislocation Love numbers and Green's functions in the time domain. Then, using the data calculated from the coordinate positions of the seafloor stations, the dislocation Love numbers and Green's functions are solved to obtain the calculated values of the displacement time series, which reflect the displacement changes of the seafloor stations after the earthquake.

[0086] S105: Determine the optimal transient viscosity parameter and the optimal steady-state viscosity parameter based on the observed values and calculated values of the displacement time series, and jointly form the best earth radial parameter with the transient viscosity parameter and the optimal steady-state viscosity parameter.

[0087] Specifically, the least squares principle is used to constrain and adjust the observed values and calculated values of the displacement time series to determine the optimal calculated values of the displacement time series that minimize the root mean square error between the observed values and calculated values. This process minimizes the difference between the observed values and calculated values by adjusting the model parameters. Next, the transient viscosity parameter and the steady-state viscosity parameter corresponding to the optimal calculated values of the displacement time series in the Burgers viscoelastic earth model are respectively determined as the optimal transient viscosity parameter and the optimal steady-state viscosity parameter. These two optimal parameters jointly form the best earth radial parameter, which can accurately describe the viscoelastic properties of the earth's medium and provide important parameter support for earthquake research and geophysical modeling.

[0088] In an alternative embodiment, refer to Figure 2 as shown Figure 2 shows a flowchart of a method for obtaining original observation data provided in Embodiment 1 of the present invention. Among them, obtaining the original observation data of several seafloor stations in the epicentral area of the submarine earthquake includes steps S201 to S203:

[0089] S201: Determine the actual seafloor station information of the research area according to the location information of the research area.

[0090] Specifically, according to the location information of the research area, the information of all seafloor stations in the area is determined, that is, the actual seafloor station information of the research area, including at least: the coordinate position information of the seafloor stations, the tectonic plates where each seafloor station is located, and the quantity information of the seafloor stations. For example, when studying a major earthquake that occurred in a certain area in 2011, the research area is the area where the earthquake occurred, and it is necessary to determine all the seafloor stations in the area and their relevant information. This information can be obtained by referring to relevant materials, databases or cooperating with research institutions. Determining the actual seafloor station information of the research area is a prerequisite for subsequent steps, ensuring the integrity and accuracy of the data.

[0091] S202: Screen out the seafloor station information in the epicentral area from the actual seafloor station information of the research area.

[0092] Specifically, according to the actual seafloor station information in the research area, the seafloor stations located in the epicentral area are screened out. The epicentral area generally refers to a certain range around the earthquake epicenter, and this range can be determined according to research needs and the scale of the earthquake. The information of the seafloor stations in the epicentral area that is screened out includes at least: the name of the seafloor station, the coordinate position information of the seafloor station, the tectonic plates where each seafloor station is located, and the quantity information of the seafloor stations.

[0093] S203: Extract the original observation data of each seafloor station from the information of the seafloor stations in the epicentral area. Among them, the original observation data includes GNSS data, attitude data, and acoustic ranging data of each seafloor station in the seafloor earthquake epicentral area.

[0094] Specifically, according to the information of the seafloor stations in the epicentral area, extract the original observation data of each seafloor station, that is, the GNSS-A data of each seafloor station, including GNSS data, attitude data, and acoustic ranging data. The GNSS data is used to determine the position and movement of the seafloor station, the attitude data is used to describe the attitude changes of the seafloor station, and the acoustic ranging data is used to measure the distance between the seafloor station and the surrounding environment. These original observation data provide a basis for subsequent coordinate position calculation and displacement time series extraction, ensuring the accuracy and reliability of the data. The extracted original observation data needs to be preprocessed, such as data cleaning, noise removal, and data correction, to improve the quality and usability of the data.

[0095] For example, see Figure 3 as shown Figure 3 FIG. shows a schematic diagram of the positions of some seafloor stations provided in the first embodiment of the present invention. Among them, the original observation data of each seafloor station obtained is from the original seafloor GNSS-A data of 6 seafloor points publicly released by the JCG (the Group of Japan Coast Guard) team in the past decade (2011.3 - 2020.6). In the figure, the positions of each seafloor station are marked with latitude (°N) and longitude (°E), the five-pointed star represents the epicenter position, the light-colored arrow represents the calculated value, the dark-colored arrow represents the observed value, the displacement rate of the seafloor station on the left of the figure is 3 cm / a, and the displacement rate of the seafloor station on the right of the figure is 1.5 cm / a. The data of the FUKU, MYGI, and MYGW stations start from March 28, 2011, and end in mid-June 2020, with 35 - 36 groups of data for each station; the data of the KAMN and KAMS stations start on April 4, 2011, and end in mid-June 2020, with 26 and 29 groups of data respectively; the observation time of the CHOS station is from April 17, 2011, to late June 2020, with a total of 28 groups of data; the observation period involved in the 6 stations is from March 28, 2011, to June 30, 2020, during which a total of 4 earthquakes with a moment magnitude Mw of more than 7.0 occurred.

[0096] In an alternative embodiment, see Figure 4 as shown, Figure 4 FIG. shows a flowchart of a method for determining an observation value of a displacement time series provided in the first embodiment of the present invention. Among them, determining the seafloor station coordinate position solution data and the observation value of the displacement time series based on the original observation data of each seafloor station includes steps S401 to S402:

[0097] S401: According to the original observation data of each seafloor station, use the sound speed spatio-temporal variation error compensation positioning model to calculate the coordinate position of the seafloor station to obtain the seafloor station coordinate position solution data.

[0098] Specifically, according to the original observation data of each seafloor station, based on the ray tracing technology, use the sound speed spatio-temporal variation error compensation positioning model to calculate the coordinate position of each seafloor station in each period to obtain the seafloor station coordinate position solution data. The sound speed spatio-temporal variation error compensation positioning model can compensate for the variation of the sound speed under different time and space conditions, thereby improving the positioning accuracy. Through this calculation process, the seafloor station coordinate position solution data can be obtained, providing a basis for the subsequent extraction of the displacement time series.

[0099] S402: According to the seafloor station coordinate position solution data, screen and extract the displacement time series of each seafloor station to determine the observation value of the displacement time series.

[0100] Specifically, collect the coordinate position solution data of the seafloor station. These data usually include the three-dimensional coordinate information of each seafloor station at different time points. Determine the research time range, for example, data from a certain time point before the earthquake to several years after the earthquake. For each seafloor station, based on its coordinate position solution data, construct its displacement time series within the research time range. This usually involves comparing the coordinates at each time point with the initial coordinates (coordinates before the earthquake) to calculate the displacement vector. Calculate the displacement of each seafloor station at each time point. The displacement can be calculated by the following formula: According to the research requirements, screen out the displacement components of interest. For example, if the research is on vertical displacement, only extract the vertical displacement components. Ensure the integrity and consistency of the data, and process missing data or outliers. Store the screened and extracted displacement time series as observation values, usually in the form of a table or a time series database for subsequent analysis and modeling.

[0101] In an alternative embodiment, see Figure 5 as shown, Figure 5 FIG. shows a flowchart of a method for constructing a viscoelastic earth model provided in the first embodiment of the present invention. Among them, determining the reset dislocation Love number and the complex Green's function based on the Burgers viscoelastic earth model includes steps S501 to S502:

[0102] S501: Introduce a complex variable parameter s into the stress-strain constitutive relation of the Burgers viscoelastic earth model, and perform Laplace transform to obtain the constitutive relation in the complex domain.

[0103] Specifically, before performing step S501, construct a Burgers viscoelastic earth model by combining a Kelvin body and a Maxwell body according to the dislocation theory of viscoelastic semi-infinite space, as shown in Figure 6 shown. Figure 6 FIG. shows a schematic diagram of a viscoelastic earth model provided in the first embodiment of the present invention. Fault rupture (shown as the rupture area in the figure) is generated by the subduction and extrusion of the land plate and the seabed plate (when two plates meet, one plate is inserted under the other relatively passive plate, and this inserted plate is the subduction plate shown in the figure). Among them, the Burgers viscoelastic earth model refers to the existing continental crust model and the oceanic crust model of CRUST 1.0, and takes the thickness of the elastic lithosphere in the area where the seabed reference station is located as 25 km deep, that is, 0-25 km is the pure elastic layer, and 25 km to the core-mantle boundary is the viscous structure layer; the Burgers viscoelastic model includes two viscosity parameters, transient and steady state, and takes them as the earth's radial parameters. is the viscous parameter of the viscoelastic land mantle wedge (i.e., the viscoelastic mantle wedge) in the Maxwell body, is the viscous parameter of the viscoelastic ocean mantle wedge in the Maxwell body, is the viscous parameter of the viscoelastic land mantle wedge in the Kelvin body, is the viscous parameter of the viscoelastic ocean mantle wedge in the Kelvin body. is the viscous parameter of the Kelvin body, is the viscous parameter of the Maxwell body, is the Lame constant of the Kelvin body, is the Lame constant of the Maxwell body.

[0104] In step S501, introduce a complex variable parameter s into the stress-strain constitutive relation of the viscoelastic model; where is the stress-strain constitutive relation formula of the viscoelastic model, is the stress, and the dot above represents the derivative with respect to time t, and are the Lame constants, is the viscous parameter, I is the identity matrix, is the trace of the matrix, and the subscript represents the viscoelastic model; It is a function obtained by differentiating displacement with respect to time; perform Laplace transform on the stress-strain constitutive relation of the viscoelastic model with the introduced complex parameter s to obtain the constitutive relation in the complex domain; among them, a new complex parameter s is introduced, and the expression in the complex domain obtained by performing Laplace transform is .

[0105] S502: Solve the target equations in the complex domain according to the constitutive relation in the complex domain to obtain the reduced dislocation Love number and the complex Green's function.

[0106] Specifically, solve the equations in the complex domain composed of the source function, the initial conditions at the center of the sphere, the inner core-outer core boundary conditions, the outer core-mantle boundary conditions, the free surface boundary conditions of the earth's surface, and the source function, etc., to obtain the reduced dislocation Love number and the complex Green's function generated by four independent point sources in the complex domain; among them, is the final solution formula for the reduced dislocation Love number and the complex Green's function; among them, the reduced dislocation Love number and the complex Green's function are very similar to the elastic constitutive relation, and the Lame parameters λ(s) and μ(s) of the reduced dislocation Love number and the complex Green's function are complex functions, denoted by .

[0107] Among them, r is the distance from the center of the sphere, θ is the co-latitude, is the longitude, t is the time, y1 and y3 are the radial and horizontal components of the displacement, the superscript T represents toroidal deformation, n is the order, m is the degree, , , are all surface vector spherical harmonics, and any vector on the unit sphere can be expressed by these three spherical harmonics.

[0108] In an alternative embodiment, the target equations include:

[0109] The initial conditions function at the center of the sphere, the inner core-outer core boundary conditions function, the outer core-mantle boundary conditions function, the free surface boundary conditions function of the earth's surface, and the source function.

[0110] In an alternative embodiment, as shown in Figure 7 shown, Figure 7 shows a flowchart of the calculated value of a displacement time series provided in the first embodiment of the present invention. Among them, based on the reduced dislocation Love number and the complex Green's function, the calculated value of the displacement time series is determined according to the observed value of the displacement time series, including steps S701~S702:

[0111] S701: Use the Riemann-Mellin inversion formula to perform an inverse transform on the said dislocation Love number and the said complex Green's function to obtain the dislocation Love number and the Green's function in the time domain.

[0112] Specifically, use the Riemann-Mellin inversion formula to perform an inverse transform on the dislocation Love number and the complex Green's function to obtain the dislocation Love number and the Green's function in the time domain. The Riemann-Mellin inversion formula is a mathematical tool that can transform a function in the complex domain into a function in the time domain, enabling complex functions in the complex domain to be compared and analyzed with actual observed data. Through this inverse transform process, the dislocation Love number and the Green's function in the time domain can be obtained, providing a basis for subsequent calculation of the displacement time series.

[0113] S702: Use the data of the calculated seabed station coordinate positions to solve the said dislocation Love number and the said Green's function to obtain the calculated value of the said displacement time series.

[0114] Specifically, according to the data of the calculated seabed station coordinate positions, solve the dislocation Love number and the Green's function to obtain the calculated value of the displacement time series. These calculated values reflect the displacement changes of the seabed station at different time points and are theoretical values calculated based on the Burgers viscoelastic earth model, the dislocation Love number, and the complex Green's function. By comparing with the observed values, the accuracy of the model and the rationality of the parameters can be evaluated, providing an important basis for subsequent parameter inversion and model optimization.

[0115] Specifically, according to the said dislocation Love number and the complex Green's function, perform an inverse transform through the Riemann-Mellin inversion formula to calculate the displacement time series of each seabed station and use it as the calculated value of the displacement time series; among them, the said dislocation Love number and the complex Green's function can be uniformly expressed as y(s), which represent the deformation under the viscoelastic model. To further obtain y(t) in the time domain, a further Laplace inverse transform is required; the Riemann-Mellin inversion formula is .

[0116] where i is the imaginary unit, c needs to be greater than all the singularities of the integrand, and s is a complex variable parameter.

[0117] In an alternative implementation, refer to Figure 8 as shown Figure 8The flowchart of an optimal parameter determination method provided by Embodiment 1 of the present invention is shown. Among them, the optimal transient viscosity parameter and the optimal steady-state viscosity parameter are determined based on the observed values and calculated values of the displacement time series, including steps S801 to S802:

[0118] S801: Constrain and adjust the observed values and calculated values of the displacement time series to determine the optimal calculated value of the displacement time series that minimizes the root mean square error between the observed values and calculated values of the displacement time series.

[0119] Specifically, first define an error function, usually the root mean square error (RMSE), to quantify the difference between the observed value and the calculated value. By adjusting the transient viscosity parameter and the steady-state viscosity parameter in the Burgers viscoelastic earth model, calculate the displacement time series under different parameter combinations, and use optimization algorithms (such as gradient descent, genetic algorithm, or simulated annealing, etc.) to find the parameter combination that minimizes the RMSE. In this process, starting from the initial parameters, gradually adjust the parameter values, and recalculate the RMSE after each adjustment until the optimal parameter combination that minimizes the error is found.

[0120] S802: Respectively determine the transient viscosity parameter and the steady-state viscosity parameter corresponding to the optimal calculated value of the displacement time series in the Burgers viscoelastic earth model as the optimal transient viscosity parameter and the optimal steady-state viscosity parameter.

[0121] Specifically, once the optimal parameter combination is determined, the corresponding transient viscosity parameter and steady-state viscosity parameter are respectively determined as the optimal transient viscosity parameter and the optimal steady-state viscosity parameter. These parameters not only make the calculation results of the model closest to the actual observed data, but also can accurately reflect the viscoelastic characteristics of the earth medium. Through this process, the best earth radial parameters can be obtained, providing important data support for seismic research and geophysical modeling.

[0122] Specifically, according to the observed value of the displacement time series and the calculated value of the displacement time series, use the least square principle to constrain the observed value and the calculated value of the displacement time series, and adjust them to minimize the root mean square error between the theoretical calculated value and the observed value during 2 to 10 years after the earthquake, so as to obtain the transient and steady-state viscosity parameters of the Burgers viscoelastic model on the ocean side, and use them as the best earth radial parameters; among them, in the least square principle, the functional relationship between the calculated value and the observed value is defined as , is the observed value, is the parameter to be estimated, is the calculated value, is the number of parameters to be estimated, represents the parameters to be estimated, denotes the th calculated value, where . Substitute displacement time series to obtain the solution of the least - squares equation; among them, substituting the displacement time series gives

[0123]

[0124] , the solution of the least - squares equation is ; X is the design matrix, denotes the th calculated value in the th observation equation, and H is the observation value matrix.

[0125] A method for determining the optimal earth radial parameters provided by the present application uses the GNSS - A data on the seabed to calculate the displacement changes of the seabed reference points, that is, the observed values of the displacement time series, and through the complex - domain transformation of the visco - elastic seismic dislocation theory, constructs a visco - elastic earth model, and calculates the displacement changes of the seabed reference points after the earthquake, that is, the calculated values of the displacement time series, through inversion and inverse transformation, and adjusts the visco - elastic parameters of the seabed according to the least - squares constraint principle to obtain the optimal earth radial parameters, which is beneficial to the refined description of the seabed structure near the earthquake source area in the earth model, the inversion of the visco - elastic structure of the earth, and more accurate characterization of the earthquake source mechanism of the seabed earthquake. At the same time, it is also beneficial to the refined management of the spatio - temporal changes after the earthquake.

[0126] Embodiment 2

[0127] Refer to Figure 9 as shown, Figure 9 shows a schematic structural diagram of an apparatus for determining the optimal earth radial parameters provided by Embodiment 2 of the present invention, where the apparatus includes:

[0128] An observation data acquisition module 901, configured to acquire the original observation data of a plurality of seabed stations in the epicentral area of the seabed earthquake;

[0129] A displacement time series observed value determination module 902, configured to determine the seabed station coordinate position solution data and the observed values of the displacement time series based on the original observation data of each seabed station;

[0130] A visco - elastic earth model construction module 903, configured to construct a Burgers visco - elastic earth model, and determine the reset dislocation Love number and the complex Green's function based on the Burgers visco - elastic earth model, where the Burgers visco - elastic earth model includes a transient viscosity parameter and a steady - state viscosity parameter, and the transient viscosity parameter and the steady - state viscosity parameter together form the earth radial parameters;

[0131] The displacement time series calculated value determination module 904 is configured to determine the calculated value of the displacement time series based on the reset dislocation Love number and the complex Green's function according to the observed value of the displacement time series;

[0132] The optimal earth radial parameter determination module 905 is configured to determine the optimal transient viscosity parameter and the optimal steady-state viscosity parameter based on the observed value and the calculated value of the displacement time series, and jointly form the optimal earth radial parameter by combining the transient viscosity parameter and the optimal steady-state viscosity parameter.

[0133] In an optional implementation manner, the obtaining of the original observation data of a plurality of undersea stations in the epicenter area of the undersea earthquake includes:

[0134] Determining the actual undersea station information of the research area according to the position information of the research area;

[0135] Screening the undersea station information in the epicenter area from the actual undersea station information of the research area;

[0136] Extracting the original observation data of each undersea station from the undersea station information in the epicenter area, wherein the original observation data includes GNSS data, attitude data, and acoustic ranging data of each undersea station in the epicenter area of the undersea earthquake.

[0137] In an optional implementation manner, the determining of the undersea station coordinate position calculation data and the observed value of the displacement time series based on the original observation data of each undersea station includes:

[0138] Calculating the undersea station coordinate position by using the sound speed spatio-temporal variation error compensation positioning model according to the original observation data of each undersea station to obtain the undersea station coordinate position calculation data;

[0139] Screening and extracting the displacement time series of each undersea station according to the undersea station coordinate position calculation data to determine the observed value of the displacement time series.

[0140] In an optional implementation manner, the determining of the reset dislocation Love number and the complex Green's function based on the Burgers viscoelastic earth model includes:

[0141] Introducing a complex variable parameter s into the stress-strain constitutive relation of the Burgers viscoelastic earth model and performing Laplace transform to obtain the constitutive relation in the complex domain;

[0142] Solving the target equation set in the complex domain according to the constitutive relation in the complex domain to obtain the reset dislocation Love number and the complex Green's function.

[0143] In an alternative embodiment, the target system of equations includes:

[0144] The initial condition function of the center of the sphere, the inner core - outer core boundary condition function, the outer core - mantle boundary condition function, the free surface boundary condition function of the earth's surface, and the source function.

[0145] Optionally, determining the calculated value of the displacement time series based on the observed value of the displacement time series according to the reset dislocation Love number and the complex Green's function includes:

[0146] Using the Riemann - Mellin inversion formula, performing an inverse transform on the reset dislocation Love number and the complex Green's function to obtain the dislocation Love number and the Green's function in the time domain;

[0147] Using the seabed station coordinate position solution data to solve the dislocation Love number and the Green's function to obtain the calculated value of the displacement time series.

[0148] In an alternative embodiment, determining the optimal transient viscosity parameter and the optimal steady - state viscosity parameter based on the observed value and the calculated value of the displacement time series includes:

[0149] Constraining and adjusting the observed value and the calculated value of the displacement time series to determine the optimal calculated value of the displacement time series that minimizes the root - mean - square error between the observed value and the calculated value of the displacement time series;

[0150] Determining the transient viscosity parameter and the steady - state viscosity parameter corresponding to the optimal calculated value of the displacement time series in the Burgers visco - elastic earth model as the optimal transient viscosity parameter and the optimal steady - state viscosity parameter respectively.

[0151] Embodiment III

[0152] Based on the same application concept, refer to Figure 10 as shown, Figure 10 shows a schematic structural diagram of a computer device provided in Embodiment III of the present invention. Among them, as Figure 10 shown, a computer device 1000 provided in Embodiment III of the present application includes:

[0153] A processor 1001, a memory 1002, and a bus 1003. The memory 1002 stores machine - readable instructions executable by the processor 1001. When the computer device 1000 runs, communication is carried out between the processor 1001 and the memory 1002 through the bus 1003. When the machine - readable instructions are run by the processor 1001, the steps of the optimal earth radial parameter determination method shown in Embodiment I above are executed.

[0154] Example 4

[0155] Based on the same application concept, an embodiment of the present application further provides a computer-readable storage medium, on which a computer program is stored. When the computer program is run by a processor, it executes the steps of the optimal earth radial parameter determination method described in any one of the above embodiments.

[0156] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the above-described system and device can refer to the corresponding processes in the foregoing method embodiments, and will not be elaborated herein.

[0157] The computer program product for determining the optimal earth radial parameters provided by the embodiments of the present invention includes a computer-readable storage medium storing program codes. The instructions included in the program codes can be used to execute the methods described in the foregoing method embodiments. For the specific implementation, reference can be made to the method embodiments, and will not be elaborated herein.

[0158] The optimal earth radial parameter determination device provided by the embodiments of the present invention can be specific hardware on a device or software or firmware installed on the device, etc. For the system provided by the embodiments of the present invention, its implementation principle and the technical effects produced are the same as those of the foregoing method embodiments. For the sake of brief description, for the parts not mentioned in the system embodiments, reference can be made to the corresponding content in the foregoing method embodiments. Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the foregoing-described system, device, and unit can all refer to the corresponding processes in the above method embodiments, and will not be elaborated herein.

[0159] In the embodiments provided by the present invention, it should be understood that the disclosed system and method can be implemented in other ways. The system embodiments described above are only illustrative. For example, the division of the units is only a logical function division, and there can be other division methods in actual implementation. For another example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed coupling or direct coupling or communication connection between each other can be through some communication interfaces. The indirect coupling or communication connection of the device or unit can be in an electrical, mechanical or other form.

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

[0161] In addition, each functional unit in the embodiments provided by the present invention may be integrated into one processing unit, may exist separately as individual physical units, or two or more units may be integrated into one unit.

[0162] If the described function is implemented in the form of a software functional unit and sold or used as an independent product, it may be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of this technical solution, may be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0163] It should be noted that: similar reference numerals and letters represent similar items in the following drawings. Therefore, once an item is defined in one drawing, it does not need to be further defined and explained in subsequent drawings. In addition, the terms "first", "second", "third", etc. are only used for distinguishing descriptions and cannot be understood as indicating or implying relative importance.

[0164] Finally, it should be noted that: the above-described embodiments are only specific implementation manners of the present invention, used to illustrate the technical solutions of the present invention, rather than limiting it. The protection scope of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: any person skilled in the art within the technical scope disclosed by the present invention can still modify the technical solutions recorded in the foregoing embodiments, or can easily think of changes, or perform equivalent replacements for some of the technical features; and these modifications, changes, or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention. All should be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims.

Claims

1. A method for determining the optimal earth radial parameters, characterized in that, The method includes: Obtaining the original observation data of several undersea stations in the epicentral area of the undersea earthquake; Determining the undersea station coordinate position solution data and the observed values of the displacement time series based on the original observation data of each undersea station; Constructing a Burgers viscoelastic earth model, and determining the reduced dislocation Love number and the complex Green's function based on the Burgers viscoelastic earth model, wherein the Burgers viscoelastic earth model includes a transient viscosity parameter and a steady-state viscosity parameter, and the transient viscosity parameter and the steady-state viscosity parameter together constitute the earth's radial parameter; Based on the reduced dislocation Love number and the complex Green's function, determining the calculated values of the displacement time series according to the observed values of the displacement time series; Determining the optimal transient viscosity parameter and the optimal steady-state viscosity parameter based on the observed values and the calculated values of the displacement time series, and jointly forming the optimal earth's radial parameter with the optimal transient viscosity parameter and the optimal steady-state viscosity parameter; The determining the calculated values of the displacement time series according to the observed values of the displacement time series based on the reduced dislocation Love number and the complex Green's function includes: Using the Riemann-Mellin inversion formula to perform an inverse transformation on the reduced dislocation Love number and the complex Green's function to obtain the dislocation Love number and the Green's function in the time domain; Using the undersea station coordinate position solution data to solve the dislocation Love number and the Green's function to obtain the calculated values of the displacement time series.

2. The method according to claim 1, characterized in that The obtaining the original observation data of several undersea stations in the epicentral area of the undersea earthquake includes: Determining the actual undersea station information of the research area according to the position information of the research area; Screening and obtaining the undersea station information in the epicentral area from the actual undersea station information of the research area; Extracting the original observation data of each undersea station from the undersea station information in the epicentral area, wherein the original observation data includes GNSS data, attitude data, and acoustic ranging data of each undersea station in the epicentral area of the undersea earthquake.

3. The method according to claim 1, wherein The determining the undersea station coordinate position solution data and the observed values of the displacement time series based on the original observation data of each undersea station includes: According to the original observation data of each undersea station, using the sound speed spatio-temporal variation error compensation positioning model to solve the undersea station coordinate position to obtain the undersea station coordinate position solution data; According to the undersea station coordinate position solution data, screening and extracting the displacement time series of each undersea station to determine the observed values of the displacement time series.

4. The method according to claim 1, wherein The determining the reduced dislocation Love number and the complex Green's function based on the Burgers viscoelastic earth model includes: Introducing a complex variable parameter s into the stress-strain constitutive relationship of the Burgers viscoelastic earth model and performing a Laplace transform to obtain the constitutive relationship in the complex domain; Solving the target equations in the complex domain according to the constitutive relationship in the complex domain to obtain the reduced dislocation Love number and the complex Green's function.

5. The method according to claim 4, wherein The target equations include: The initial condition function of the sphere center, the inner core - outer core boundary condition function, the outer core - mantle boundary condition function, the free surface boundary condition function of the earth's surface, and the source function.

6. The method according to claim 1, wherein Determining the optimal transient viscosity parameter and the optimal steady - state viscosity parameter based on the observed values and calculated values of the displacement time series, including: Constraining and adjusting the observed values and calculated values of the displacement time series to determine the optimal calculated value of the displacement time series that minimizes the root - mean - square error between the observed values and calculated values of the displacement time series; Respectively determining the transient viscosity parameter and the steady - state viscosity parameter corresponding to the optimal calculated value of the displacement time series in the Burgers visco - elastic earth model as the optimal transient viscosity parameter and the optimal steady - state viscosity parameter.

7. An apparatus for determining optimal earth radial parameters, characterized in that, The device includes: An observation data acquisition module for acquiring the original observation data of several seafloor stations in the epicentral area of the submarine earthquake; A displacement time series observed value determination module for determining the seafloor station coordinate position solution data and the observed values of the displacement time series based on the original observation data of each seafloor station; A visco - elastic earth model construction module for constructing a Burgers visco - elastic earth model, and determining the reduced dislocation Love number and the complex Green's function based on the Burgers visco - elastic earth model, where the Burgers visco - elastic earth model includes a transient viscosity parameter and a steady - state viscosity parameter, and the transient viscosity parameter and the steady - state viscosity parameter jointly form the earth's radial parameters; A displacement time series calculated value determination module for determining the calculated values of the displacement time series based on the reduced dislocation Love number and the complex Green's function according to the observed values of the displacement time series; An optimal earth radial parameter determination module for determining the optimal transient viscosity parameter and the optimal steady - state viscosity parameter based on the observed values and calculated values of the displacement time series, and jointly forming the optimal transient viscosity parameter and the optimal steady - state viscosity parameter into the best earth radial parameters; Determining the calculated values of the displacement time series based on the reduced dislocation Love number and the complex Green's function according to the observed values of the displacement time series, including: Using the Riemann - Mellin inversion formula to perform an inverse transformation on the reduced dislocation Love number and the complex Green's function to obtain the dislocation Love number and the Green's function in the time domain; Using the seafloor station coordinate position solution data to solve the dislocation Love number and the Green's function to obtain the calculated values of the displacement time series.

8. A computer device, characterized in that, Including: A processor, a memory, and a bus. The memory stores machine - readable instructions executable by the processor. When the computer device runs, the processor communicates with the memory through the bus. When the machine - readable instructions are executed by the processor, the steps of the best earth radial parameter determination method as described in any one of claims 1 to 6 are executed.

9. A computer-readable storage medium, characterized in that, A computer program is stored on the computer - readable storage medium. When the computer program is run by the processor, the steps of the best earth radial parameter determination method as described in any one of claims 1 to 6 are executed.

Images