Four-component borehole strain observation data small window analysis method and related device

Abstract

The application discloses a small window analysis method for four-component borehole strain observation data and a related device, relates to the technical field of four-component borehole strain observation data analysis, and comprises the following steps: segmenting four-component borehole strain observation data to be analyzed to obtain four-component borehole strain observation data in multiple unit time periods; performing data analysis on the four-component borehole strain observation data in the unit time period; determining observation errors of observed relative tectonic ground stresses and errors of non-tectonic ground stresses caused by each influencing factor; further determining errors between observed stresses and relative tectonic ground stresses; further determining a value range of a first slope of linear change of the observed relative tectonic ground stress with time; and selecting a slope in the value range of the first slope as a second slope of linear change of the relative tectonic ground stress with time by using a linear elasticity assumption, so that the relative tectonic ground stress in the four-component borehole strain observation data can be extracted.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of four-component borehole strain observation data analysis, in particular to a small window analysis method for four-component borehole strain observation data and a related device. BACKGROUND

[0002] Scientific processing and interpretation of four-component borehole strain observation data is a key scientific problem that needs to be urgently solved in the fields of geodynamics and seismology. At present, most scholars focus on the data processing method and application method of four-component borehole strain observation data, that is, the research on the observation phenomenon in the four-component borehole strain observation data still stays in the stage of interference exclusion for a specific observation phenomenon, and there is no method to extract relative tectonic stress. SUMMARY

[0003] The purpose of the present application is to provide a small window analysis method for four-component borehole strain observation data and a related device, which can extract relative tectonic stress in four-component borehole strain observation data.

[0004] To achieve the above purpose, the present application provides the following solutions:

[0005] In a first aspect, the present application provides a small window analysis method for four-component borehole strain observation data, which comprises:

[0006] Obtaining four-component borehole strain observation data to be analyzed; wherein the four-component borehole strain observation data to be analyzed comprises observation stresses at multiple time points, the observation stress is the sum of an observed relative tectonic stress and a non-tectonic stress caused by each influencing factor, and the influencing factor is a factor that affects the observation of the four-component borehole strain observation data;

[0007] Segmenting the four-component borehole strain observation data to be analyzed to obtain four-component borehole strain observation data in multiple unit time periods; wherein the observed relative tectonic stress in the four-component borehole strain observation data in the unit time period changes linearly with time, and the non-tectonic stress caused by each influencing factor changes nonlinearly with time;

[0008] For each unit time period of four-component borehole strain observation data, data analysis is performed on the four-component borehole strain observation data in the unit time period to determine the observation error of the observed relative tectonic stress and the error of the non-tectonic stress caused by each influencing factor, and based on the observation error of the observed relative tectonic stress and the error of the non-tectonic stress caused by each influencing factor, the error between the observation stress and the relative tectonic stress is determined;

[0009] For each four-component borehole strain observation data of a unit time period, a range of a first slope of linear change of the observed relative tectonic crustal stress with time is determined based on an error between the observed stress and the relative tectonic crustal stress, a second slope of linear change of the relative tectonic crustal stress with time is selected as one of the range of the first slope by using a linear elasticity assumption, and the relative tectonic crustal stress at each time of the unit time period is determined based on the second slope.

[0010] Optionally, the influence factors include temperature, water level, earth tide and borehole stress relaxation.

[0011] Optionally, the unit time period is 1 hour or 2 hours.

[0012] Optionally, a calculation formula of the error between the observed stress and the relative tectonic crustal stress is:

[0013]

[0014] wherein, δσ1 is the error between the observed stress and the relative tectonic crustal stress; δσ0 is an observation error of the observed relative tectonic crustal stress; m is a total number of the influence factors; δσ 1k is an error of the non-tectonic crustal stress caused by the kth influence factor.

[0015] Optionally, the range of the first slope of linear change of the observed relative tectonic crustal stress with time is determined based on the error between the observed stress and the relative tectonic crustal stress, and specifically includes:

[0016] An upper limit and a lower limit of the observed stress at an end time of the unit time period are determined based on the error between the observed stress and the relative tectonic crustal stress; the upper limit is a sum of the observed stress at the end time and the error between the observed stress and the relative tectonic crustal stress, and the lower limit is a difference between the observed stress at the end time and the error between the observed stress and the relative tectonic crustal stress;

[0017] A first linear expression of linear change of the observed stress with time is established; a slope distance of the first linear expression on a y-axis is the observed stress at an initial time of the unit time period;

[0018] The upper limit and the lower limit are taken as inputs, the range of the slope of linear change of the observed stress with time is determined by using the first linear expression, and the range of the slope of linear change of the observed stress with time is taken as the range of the first slope of linear change of the observed relative tectonic crustal stress with time.

[0019] Optionally, a calculation formula of the second slope of linear change of the relative tectonic crustal stress with time is:

[0020] η=σ x σ y-(2C-σ xy )(2D-σ xy ;

[0021] wherein, η is a slope, min|η| is a second slope of the relative tectonic stress linearly changing with time; σ x is a normal stress in the x direction, determined based on a value range of the first slope; σ y is a normal stress in the y direction, determined based on a value range of the first slope; C is a constant, C = σ x0 + σ xy0 , σ x0 is a normal stress in the x direction at an initial moment in a unit time period, σ xy0 is a shear stress at the initial moment in the unit time period; σ xy is a shear stress, determined based on a value range of the first slope; D = σ y0 + σ xy0 , σ y0 is a normal stress in the y direction at the initial moment in the unit time period.

[0022] Optionally, based on the second slope, the relative tectonic stress at each moment in a unit time period is determined, and specifically includes:

[0023] a second linear expression of the relative tectonic stress linearly changing with time is established; a slope of the second linear expression is the second slope, and an ordinate of the second linear expression is the relative tectonic stress at an initial moment in a unit time period;

[0024] the relative tectonic stress at each moment in the unit time period is determined by using the second linear expression.

[0025] In a second aspect, a computer device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor executes the computer program to implement the four-component borehole strain observation data small window analysis method.

[0026] In a third aspect, a computer readable storage medium is provided, and the computer readable storage medium stores a computer program, and the computer program is executable on a processor to implement the four-component borehole strain observation data small window analysis method.

[0027] In a fourth aspect, a computer program product is provided, and the computer program product includes a computer program, and the computer program is executable on a processor to implement the four-component borehole strain observation data small window analysis method.

[0028] According to the embodiments provided in the present application, the present application has the following technical effects:

[0029] The application provides a small window analysis method for four-component borehole strain observation data and a related device. BRIEF DESCRIPTION OF DRAWINGS

[0030] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed in the embodiments will be briefly introduced as follows. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without any creative effort on the basis of these drawings.

[0031] Figure 1 A flowchart of a small window analysis method for four-component borehole strain observation data provided in Embodiment 1 of the present application.

[0032] Figure 2 A schematic diagram of four-component borehole strain observation data provided in Embodiment 1 of the present application.

[0033] Figure 3 A schematic diagram of the range of observed stress and first slope provided in Embodiment 1 of the present application.

[0034] Figure 4 A schematic diagram of the change of observed stress in a small time scale provided in Embodiment 1 of the present application.

[0035] Figure 5 A schematic diagram of a linear elastic small unit provided in Embodiment 1 of the present application.

[0036] Figure 6 A structural schematic diagram of a computer device provided in Embodiment 2 of the present application. DETAILED DESCRIPTION

[0037] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work fall within the scope of the present application.

[0038] Embodiment 1

[0039] In an exemplary embodiment, as shown in Figure 1 A small window analysis method of four-component borehole strain observation data is provided, comprising the following steps.

[0040] Step S1, obtaining four-component borehole strain observation data to be analyzed; wherein the four-component borehole strain observation data to be analyzed includes observation stresses at multiple time points, the observation stress is the sum of observation relative tectonic stress and non-tectonic stress caused by each influencing factor, and the influencing factor is a factor that affects the observation of the four-component borehole strain observation data.

[0041] Step S2, segmenting the four-component borehole strain observation data to be analyzed to obtain four-component borehole strain observation data in multiple unit time periods; wherein the observation relative tectonic stress in the four-component borehole strain observation data in the unit time period changes linearly with time, and the non-tectonic stress caused by each influencing factor changes nonlinearly with time.

[0042] Step S3, for each unit time period of four-component borehole strain observation data, performing data analysis on the four-component borehole strain observation data in the unit time period to determine the observation error of the observation relative tectonic stress and the error of the non-tectonic stress caused by each influencing factor, and based on the observation error of the observation relative tectonic stress and the error of the non-tectonic stress caused by each influencing factor, determining the error between the observation stress and the relative tectonic stress.

[0043] Step S4, for each unit time period of four-component borehole strain observation data, based on the error between the observation stress and the relative tectonic stress, determining the value range of the first slope of the linear change of the observation relative tectonic stress with time, selecting a slope in the value range of the first slope as the second slope of the linear change of the relative tectonic stress with time based on the linear elasticity assumption, and determining the relative tectonic stress at each time point in the unit time period based on the second slope.

[0044] By implementing the steps S1 to S4, the embodiment can determine the relative tectonic stress at each moment, which not only has an important reference role for improving the monitoring of the stress tensor and the research method system of the earthquake numerical prediction, but also has important scientific significance and application value for improving the understanding of the regional stress field and the stress evolution characteristics and laws of the fault, and the understanding of the seismic dynamic process and the earthquake preparation.

[0045] The embodiment provides a small window analysis method for four-component borehole strain observation data, which is mainly used for processing the observed stress which may represent the observed relative tectonic stress change in the four-component borehole strain observation data, extracting the change of the relative tectonic stress (referring to a force source which has a key influence on the stress concentration in a small range of the plate, such as the fracture complexity, the deep stratum drag and the like), and in the research related to the earthquake prediction, the key influencing factors which cause the key abnormality are usually sought, which may be temperature, water, other seismic factors and the like, and it is hoped that the key influencing factors which cause the key abnormality are related to the earthquake event, but the existing problem is that there is no applicable method for analyzing the change of the relative tectonic stress in the four-component borehole strain observation data, the small window analysis can be applied to the four-component borehole strain observation data before the earthquake, the change of the relative tectonic stress before the earthquake is extracted, the influence of the larger scale geodynamic force source on the regional relative tectonic stress field is reflected, and the subsequent earthquake prediction research is facilitated.

[0046] The small window analysis in the embodiment refers to the process of obtaining the change of the relative tectonic stress under specific conditions and using specific judgment criteria, and the judgment criteria have originality. The judgment criteria are mainly used for determining the slope of the linear change of the relative tectonic stress with time, and the slope represents the influence of the relative tectonic stress on the next moment under certain background, so as to determine the relative tectonic stress at each moment.

[0047] The small window analysis method for four-component borehole strain observation data provided by the embodiment specifically includes the following steps.

[0048] (1) Obtain the four-component borehole strain observation data to be analyzed.

[0049] The four-component borehole strain observation data will be affected by some non-tectonic stress caused by influencing factors such as temperature, water level, solid tide and borehole stress relaxation. After the non-tectonic stress caused by the influencing factors is removed, the relative tectonic stress can be determined.

[0050] In the embodiment, the four-component borehole strain observation data to be analyzed includes observation stresses at multiple time instants, the observation stress is the sum of the observation relative tectonic stress (i.e., the observed relative tectonic stress, which is equivalent to the relative tectonic stress observation value) and the non-tectonic stress caused by each influencing factor, and the influencing factor is a factor that affects the observation of the four-component borehole strain observation data, and the influencing factor can include temperature, water level, solid tide, and borehole stress relaxation.

[0051] (2) The four-component borehole strain observation data to be analyzed is processed in segments to obtain four-component borehole strain observation data in multiple unit time periods.

[0052] For the four-component borehole strain observation data, the following assumptions (small window assumption) are made in the embodiment:

[0053] 1) The four-component borehole strain observation data observes the change of the relative tectonic stress.

[0054] 2) The change of the relative tectonic stress in a certain time period can be assumed to be linear, and the change of the non-tectonic stress caused by other influencing factors in the time period can be fitted by a nonlinear function.

[0055] Taking the observation stress that can be observed in a period of time as an example, as shown in Figure 2 , the observation result of the observation stress in a period of time, Figure 2 , where σ is the observation stress, t is the time, t0 to t1 is the unit time period to be analyzed, and in the unit time period, the observation relative tectonic stress (affected by more macro factors such as plate movement and fracture system complexity) is assumed to be linear, while the change of the non-tectonic stress caused by the influencing factors such as borehole stress relaxation and temperature is nonlinear. In t0 to t1, the observation stress σ0 at t0 and the observation stress σ1 at t1 are determined, and the relative tectonic stress can be obtained, but in fact, the observation stress at t1 is affected by all influencing factors, including the non-tectonic stress caused by each influencing factor, so if the change of the relative tectonic stress is assumed to be linear, the relative tectonic stress can be extracted by analysis.

[0056] In the embodiment, the four-component borehole strain observation data to be analyzed is processed in segments to obtain four-component borehole strain observation data in multiple unit time periods, wherein the observation relative tectonic stress in the four-component borehole strain observation data in the unit time period changes linearly with time, and the non-tectonic stress caused by each influencing factor changes nonlinearly with time, and the unit time period can be 1 hour or 2 hours.

[0057] (3) For the four-component borehole strain observation data of each unit time period, perform data analysis on the four-component borehole strain observation data of each unit time period, determine the observation error of the relative tectonic stress and the error of the non-tectonic stress caused by each influencing factor, and determine the error between the observed stress and the relative tectonic stress based on the observation error of the relative tectonic stress and the error of the non-tectonic stress caused by each influencing factor.

[0058] In this embodiment, the process of analyzing the four-component borehole strain observation data over a unit time period to determine the observation error of the relative tectonic stress and the error of the non-tectonic stress caused by each influencing factor can be carried out using any existing method, such as wavelet analysis, Fourier transform, etc. This embodiment will not elaborate on this process.

[0059] In this embodiment, the formula for calculating the error between the observed stress and the relative tectonic stress is as follows:

[0060]

[0061] In equation (1), δσ1 is the error between the observed stress and the relative tectonic stress; δσ0 is the observation error of the relative tectonic stress; m is the total number of influencing factors; δσ 1k The error in non-tectonic stress caused by the kth influencing factor.

[0062] (4) For the four-component borehole strain observation data of each unit time period, based on the error between the observed stress and the relative tectonic stress, the range of the first slope of the linear change of the observed relative tectonic stress with time is determined. Using the linear elastic assumption, a slope is selected within the range of the first slope as the second slope of the linear change of the relative tectonic stress with time. Based on the second slope, the relative tectonic stress at each moment in the unit time period is determined.

[0063] Under the influence of a series of factors, such as Figure 3 As shown, Figure 3 In the diagram, the solid black line represents the relative tectonic stress to be calculated, the solid colored line represents the non-tectonic stress caused by other possible influencing factors, and the dashed black line represents the upper and lower limits of the solid black line. The observed stress σ1 at the final time t1 has an upper and lower limit, forming a set of observed stress changes G(σ,t). Assuming the observed stress changes linearly, and the observed stress σ0 is known at the initial time t0, the slope on the y-axis is σ0. Only a slope needs to be determined. Based on the set of observed stress changes G(σ,t), a set of slopes A' can be obtained, i.e. Figure 3A' is a set of slopes, which is the range of the slope of the linear function of the observed stress with time. Then, a proper slope is selected from A' as the second slope of the linear function of the relative tectonic stress with time, which is the best one that excludes all possible influences.

[0064] In this embodiment, the range of the first slope of the linear function of the observed relative tectonic stress with time is determined based on the error between the observed stress and the relative tectonic stress, and specifically includes:

[0065] 1) The upper limit and the lower limit of the observed stress at the end time in a unit time period are determined based on the error between the observed stress and the relative tectonic stress, the upper limit being the sum of the observed stress at the end time and the error between the observed stress and the relative tectonic stress, and the lower limit being the difference between the observed stress at the end time and the error between the observed stress and the relative tectonic stress.

[0066] In formula (1), δσ1 calculated is positive, at this time, the upper limit is σ1+δσ1, and the lower limit is σ1-δσ1.

[0067] 2) The first linear expression of the observed stress with time is established, the slope of the first linear expression on the y-axis being the observed stress at the initial time in a unit time period.

[0068] In the formula, x-axis represents time, and y-axis represents the observed stress, and the first linear expression is as follows:

[0069] G(σ, t) = A't + σ0 (2)

[0070] In formula (2), G(σ, t) is a set of observed stress changes composed of the upper limit and the lower limit of the observed stress; A' is the range of the slope; and σ0 is the observed stress at the initial time in a unit time period.

[0071] 3) The range of the slope of the linear function of the observed stress with time is determined by using the first linear expression with the upper limit and the lower limit as inputs, and the range of the slope of the linear function of the observed stress with time is taken as the range of the first slope of the linear function of the observed relative tectonic stress with time.

[0072] How to find a reasonable slope from A' is to determine the target of the relative tectonic stress. It is very difficult to get the relative tectonic stress at the end time t1, but since there is the assumption that the relative tectonic stress changes linearly with time, δσ changing within δt near the initial time t0 can be determined arbitrarily, as shown in FIG. 3, on the time scale of δt, the rock can be assumed to be linearly elastic, and therefore, the borehole together with the nearby rock mass can be regarded as a linearly elastic small unit, as shown in FIG. 4. Figure 4 Figure 5

[0073] ​​Based on the plane stress assumption and the basic principles of elastic mechanics, the linear elastic small element composed of dx and dy should satisfy:

[0074] σ ij，j +f i = 0 in V (3)

[0075]

[0076] In the above formula, formula (3) is the balance equation in the domain, σ ij,j (i = x, y, j = x, y) represents the stress σ ij in the j direction; f i (i = x, y, z) represents the body force, since the analysis of four-component borehole strain observation data is all based on the plane stress assumption, f i = 0; V represents the linear elastic small element composed of dx and dy to be calculated; formula (4) is the geometric equation, ε ij (i = x, y, j = x, y) represents the strain; u i,j (i = x, y, j = x, y) and u j,i (i = x, y, j = x, y) represent the derivative of displacement; formula (5) is the constitutive equation, σ ij (i = x, y, j = x, y) represents the stress, which describes the relationship between the strain tensor and the strain energy density; w represents the strain energy density; formula (6) is the stress boundary condition, represents the stress of the boundary; S t represents the boundary.

[0077] In a relatively small time interval δt, the linear elastic small element is affected by the relative tectonic stress, and the change of dx direction is δu, and the change of dy direction is δv, and the corresponding strain is: δε x is the x-direction strain, δu is the x-direction displacement, δε y is the y-direction strain, δv is the y-direction displacement, δε xy is the shear strain, and the stress of the boundary of the linear elastic small element is the relative tectonic stress σ x0 at the initial time, σ y0 , σ xy0 , then the following can be obtained:

[0078] δw 外 = ∫(σ x0 + σ xy0 ) δu dx + ∫(σ y0 + σ xy0 ) δv dy (7)

[0079]

[0080] In the above formula, δw外 Work done by external force in δt; σ x0 Normal stress in x direction at initial moment in unit time period; σ xy0 Shear stress at initial moment in unit time period; σ y0 Normal stress in y direction at initial moment in unit time period; δw 内 Work done by internal force in δt; δσ x Normal stress in x direction in δt; δσ y Normal stress in y direction in δt; δσ xy Shear stress in δt.

[0081] Since δw 内 = δw 外 , let σ x0 + σ xy0 = C, σ y0 + σ xy0 = D, C and D are constants, it is easy to obtain: σ x σ y = (2C- σ xy )(2D- σ xy ), σ x is the normal stress in x direction after δt, σ y is the normal stress in y direction after δt, and σ xy is the shear stress after δt.

[0082] Let η = σ x σ y - (2C- σ xy )(2D- σ xy ), a plurality of slopes are selected in the range of the first slope (which can be randomly selected or uniformly selected based on a preset step size), the normal stress in x direction, the normal stress in y direction and the shear stress at the initial moment are used to determine the normal stress in x direction, the normal stress in y direction and the shear stress after δt corresponding to each selected slope based on the selected slope, and η corresponding to each selected slope is further calculated, and min(|η|) is the final second slope selected at A'.

[0083] The index min(|η|) can be used as a judgment criterion as a selection standard of the second slope of the relative tectonic stress in the four-component borehole strain observation data of continuous observation, so that the change of the relative tectonic stress (mainly the plate movement and the given stress in the deep part) can be determined in a relatively short unit time period (assuming that the relative tectonic stress changes linearly).

[0084] In the embodiment, the calculation formula of the second slope of the relative tectonic stress changing linearly with time is:

[0085] η = σx σ y -(2C-σ xy )(2D-σ xy ) (9)

[0086] In formula (9), η is a slope, min|η| is a second slope of the relative tectonic stress linearly changing with time; σ x is a normal stress in the x direction, determined based on a value range of the first slope; σ y is a normal stress in the y direction, determined based on the value range of the first slope; C is a constant, C = σ x0 + σ xy0 , σ x0 is the normal stress in the x direction at an initial moment in a unit time period, σ xy0 is a shear stress at the initial moment in the unit time period; σ xy is a shear stress, determined based on the value range of the first slope; D = σ y0 + σ xy0 , σ y0 is the normal stress in the y direction at the initial moment in the unit time period.

[0087] In this embodiment, based on the second slope, the relative tectonic stress at each moment in a unit time period is determined, specifically including:

[0088] 1) a second linear expression of the relative tectonic stress linearly changing with time is established, a slope of the second linear expression is the second slope, and an abscissa of the second linear expression on the y axis is the relative tectonic stress at an initial moment in a unit time period.

[0089] 2) the relative tectonic stress at each moment in the unit time period is determined by using the second linear expression.

[0090] The second linear expression is:

[0091] L = At + B (10)

[0092] In formula (10), L is the relative tectonic stress, A is the second slope; and B is the relative tectonic stress at the initial moment.

[0093] The embodiment can determine the change of relative tectonic stress in the four-component borehole strain observation data to be analyzed, and verify whether the change of relative tectonic stress before the earthquake is abnormal. The key to the phenomenon such as earthquake is caused by the plate drag in the deep earth, or the relative action of the complex fault system and the plate, which has a relatively long time and may be the root cause of causing underground water anomaly. The relative tectonic stress change has a long period and does not change regularly, but in a relatively short time range (for example, within 1-2 hours) in the smooth observation curve (without other phenomena such as earthquakes), the relative tectonic stress is linearly changed, and other non-tectonic stresses are nonlinearly changed in the relatively short time range. Therefore, the relative tectonic stress can be extracted from the observation phenomenon to determine the change of the relative tectonic stress, which is of great significance to the prediction research of the earthquake in the hours to months before the earthquake.

[0094] The small window analysis method of the embodiment focuses on evaluating the four-component borehole strain observation data which changes dynamically by physical laws. Further, by using the operational research method and comprehensively judging the research on the data influencing factors of the predecessors, an operational and systematic data processing tool can be provided for the whole processing process and big data processing method of the four-component borehole strain observation data, and the connection between the relative tectonic stress and the seismic activity, the tectonic stress change, and the physical nature of the seismic dynamic process. Therefore, by improving and researching the data analysis method of the embodiment, the method has an important reference role for improving the monitoring of the stress tensor, the research method system of the earthquake numerical prediction, and has an important scientific significance and application value for improving the understanding of the regional stress field and the fault stress evolution characteristics and rules, and the understanding of the seismic dynamic process and the earthquake preparation.

[0095] The application also provides an application scenario of the small window analysis method of the four-component borehole strain observation data. Specifically, the small window analysis method of the four-component borehole strain observation data provided by the embodiment can be applied in the earthquake prediction research scenario. The earthquake prediction research scenario includes an analysis link and a prediction link. The analysis link is used for small window analysis of the four-component borehole strain observation data before the earthquake to determine the relative tectonic stress at each moment. The prediction link is used for earthquake prediction research based on the relative tectonic stress at each moment. The small window analysis method of the four-component borehole strain observation data provided by the embodiment belongs to the analysis link.

[0096] Embodiment 2

[0097] In an exemplary embodiment, a computer device can be provided, which can be a server or a terminal, and an internal structure diagram of the computer device can be as shown in Figure 6As shown in the figure. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. Among them, the processor, the memory and the input / output interface are connected through the system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capability. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store data. The input / output interface of the computer device is used to exchange information between the processor and the external device. The communication interface of the computer device is used to communicate with the terminal outside through the network connection. The computer program is executed by the processor to realize a small window analysis method of four-component borehole strain observation data.

[0098] Those skilled in the art can understand that, Figure 6 The structure shown in the figure is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device to which the scheme of the present application is applied. The specific computer device can include more or fewer components than those shown in the figure, or combine certain components, or have a different component arrangement.

[0099] In an exemplary embodiment, a computer device is provided, comprising a memory and a processor, the memory storing a computer program, and the processor executing the computer program to realize the small window analysis method of four-component borehole strain observation data in embodiment 1.

[0100] Embodiment 3

[0101] In an exemplary embodiment, a computer readable storage medium is provided, storing a computer program, which is executed by a processor to realize the small window analysis method of four-component borehole strain observation data in embodiment 1.

[0102] Embodiment 4

[0103] In an exemplary embodiment, a computer program product is provided, comprising a computer program, which is executed by a processor to realize the small window analysis method of four-component borehole strain observation data in embodiment 1.

[0104] It should be noted that the user information (including but not limited to user equipment information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) involved in the present application are all information and data authorized by the user or authorized by all parties, and the collection, use and processing of related data need to comply with relevant regulations.

[0105] The technical features of the above embodiments can be combined in any manner. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described, however, it should be understood that any combination of the technical features is within the scope of the present disclosure as long as there is no contradiction.

[0106] The principles and implementation manners of the present application are described herein by using specific examples, and the above embodiments are only used to help understand the method of the present application and its core idea; meanwhile, for those skilled in the art, the specific implementation manners and application range can be changed according to the idea of the present application. In conclusion, the content of the present description should not be understood as a limitation of the present application.

Claims

1. A small window analysis method of four-component borehole strain observation data, characterized by, The small window analysis method of the four-component borehole strain observation data comprises: Obtaining four-component borehole strain observation data to be analyzed; wherein the four-component borehole strain observation data to be analyzed comprises observation stresses at multiple time points, the observation stress being a sum of an observation relative tectonic stress and a non-tectonic stress caused by each influencing factor, and the influencing factor being a factor affecting the observation of the four-component borehole strain observation data; Segmenting the four-component borehole strain observation data to be analyzed to obtain four-component borehole strain observation data in multiple unit time periods; wherein the observation relative tectonic stress in the four-component borehole strain observation data in the unit time period changes linearly with time, and the non-tectonic stress caused by each influencing factor changes nonlinearly with time; For the four-component borehole strain observation data in each unit time period, performing data analysis on the four-component borehole strain observation data in the unit time period to determine an observation error of the observation relative tectonic stress and an error of the non-tectonic stress caused by each influencing factor, and determining an error between the observation stress and the relative tectonic stress based on the observation error of the observation relative tectonic stress and the error of the non-tectonic stress caused by each influencing factor; For the four-component borehole strain observation data in each unit time period, determining a value range of a first slope of the observation relative tectonic stress changing linearly with time based on the error between the observation stress and the relative tectonic stress, selecting a slope in the value range of the first slope as a second slope of the relative tectonic stress changing linearly with time by using the linear elasticity assumption, and determining the relative tectonic stress at each time point in the unit time period based on the second slope; The value range of the first slope of the observation relative tectonic stress changing linearly with time is determined based on the error between the observation stress and the relative tectonic stress, and specifically comprises: Determining an upper limit and a lower limit of the observation stress at the end time point in the unit time period based on the error between the observation stress and the relative tectonic stress; the upper limit is a sum of the observation stress at the end time point and the error between the observation stress and the relative tectonic stress, and the lower limit is a difference between the observation stress at the end time point and the error between the observation stress and the relative tectonic stress; Establishing a first linear expression of the observation stress changing linearly with time; the slope of the first linear expression on the y-axis is the observation stress at the initial time point in the unit time period; Taking the upper limit and the lower limit as inputs, determining a value range of the slope of the observation stress changing linearly with time by using the first linear expression, and taking the value range of the slope of the observation stress changing linearly with time as the value range of the first slope of the observation relative tectonic stress changing linearly with time; The calculation formula of the second slope of the relative tectonic stress changing linearly with time is: ； wherein is a slope, min is a second slope of the relative tectonic stress varying linearly with time; is a normal stress in the x-direction, determined based on a range of values of the first slope; is a normal stress in the y-direction, determined based on a range of values of the first slope; is a constant, = 0 , is a normal stress in the x-direction at an initial time in a unit time period, is a shear stress at an initial time in a unit time period; is a shear stress, determined based on a range of values of the first slope; , is a normal stress in the y-direction at an initial time in a unit time period.

2. The small window analysis method of four-component borehole strain observation data according to claim 1, characterized in that, The influencing factors include temperature, water level, solid tide and borehole stress relaxation.

3. The small window analysis method of four-component borehole strain observation data according to claim 1, characterized by, The unit time period is 1 hour or 2 hours.

4. The small window analysis method of four-component borehole strain observation data according to claim 1, characterized by, The calculation formula of the error between the observation stress and the relative tectonic stress is: ； wherein, is the error in the observed stress and the relative tectonic stress; is the observation error in the observed relative tectonic stress; is the total number of influencing factors; is the error in the non-tectonic stress caused by the k th influencing factor.

5. The small window analysis method of four-component borehole strain observation data according to claim 1, characterized by, Based on the second slope, the relative tectonic stress at each time point in the unit time period is determined, and specifically comprises: establishing a second linear expression of the relative tectonic crustal stress changing linearly with time; a slope of the second linear expression is a second slope, and an intercept of the second linear expression on a y-axis is the relative tectonic crustal stress at an initial time in a unit time period; determining the relative tectonic crustal stress at each time in the unit time period by using the second linear expression.

6. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method for small window analysis of four-component borehole strain observation data according to any one of claims 1-5.

7. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to implement the method for small window analysis of four-component borehole strain observation data according to any one of claims 1-5.

8. A computer program product comprising a computer program, characterized in that, The computer program is executed by the processor to implement the method for small window analysis of four-component borehole strain observation data according to any one of claims 1-5.

Images