Near-fault seismic oscillation baseline correction method and device matched with target permanent displacement and computer storage medium
By introducing baseline correction function and displacement slope function fitting technology in the near fault earthquake baseline correction method, the problems of deviation and insufficient recording of baseline correction results in the existing technology are solved, and the precise matching goal of near fault earthquake records is achieved. Permanent displacement correction.
Patent Information
- Application Number
- CN202411349221.X
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-09-26
- Publication Date
- 2025-02-14
AI Technical Summary
The existing near-fault earthquake baseline correction methods have limitations in removing noise and retaining permanent displacement information, resulting in a deviation from the GPS co-seismic displacement recorded after correction, and insufficient earthquake recording that meets different permanent displacement levels.
A near-fault earthquake baseline correction method that matches the target's permanent displacement is adopted. Through initial baseline correction and displacement slope function fitting, the baseline correction function is calculated to correct the baseline drift function, and the original acceleration record is re-corrected to ensure the permanent displacement recorded after correction. Precisely match preset target displacement.
Accurate baseline correction for near-fault earthquake recording is achieved, so that the permanent displacement recorded after correction matches the GPS co-seismic displacement or preset target displacement level, significantly improving the accuracy and applicability of the correction results.
Smart Images
Figure CN119439274A_ABST
Abstract
Description
Technical Field
[0001] The present application belongs to the technical field of disaster prevention and reduction and engineering earthquakes, and in particular, relates to a near-fault ground motion baseline correction method, device and computer storage medium for matching target permanent displacement. Background Art
[0002] Near-fault seismic records generally have obvious concentration and directionality due to the short distance between faults, and show the characteristics of large velocity pulses and permanent ground displacement. Near-fault seismic records generally show a linear trend deviating from the equilibrium position at the end of the velocity time history, and more serious nonlinear divergence at the end of the displacement time history. Therefore, in the baseline correction of near-fault seismic motions, it is necessary not only to remove the noise in the seismic records, but also to retain the permanent displacement information to the maximum extent, which significantly increases the difficulty of baseline correction.
[0003] The existing near-fault baseline correction method based on the Iwan two-stage correction principle (such as patent CN116224441A) can restore the permanent displacement form well, but due to the limitations of the correction parameter selection, the accuracy and stability of the restoration of the permanent displacement are affected. Therefore, the permanent displacement recorded after correction and the GPS co-seismic displacement (the closer the GPS measuring point is to the strong earthquake station, the higher the authenticity) still have a certain degree of deviation, which needs further correction. On the other hand, in the study of cross-fault bridges, seismic motions with different permanent displacement levels are required as input, and the number of existing records that meet the conditions is seriously insufficient. Therefore, how to obtain reasonable input that meets the requirements through appropriate correction of actual near-fault records is a key issue that needs to be carried out urgently.
[0004] Therefore, it is of great significance to develop a baseline correction method that can accurately match the preset target displacement by reasonably improving the principle of near-fault ground motion baseline correction. Summary of the invention
[0005] The purpose of the present application is to overcome the problems of the prior art and disclose a near-fault seismic baseline correction method, device and computer storage medium for matching target permanent displacement. The present application can accurately match the preset target permanent displacement on the basis of effectively restoring the near-fault seismic permanent displacement, so that the permanent displacement recorded near the fault after correction accurately matches the GPS co-seismic displacement or different target permanent displacement levels.
[0006] On the one hand, the purpose of this application is achieved through the following technical solutions:
[0007] A near-fault ground motion baseline correction method for matching target permanent displacement, the near-fault ground motion baseline correction method comprising:
[0008] S1: Perform initial baseline correction and fit the corrected displacement time course;
[0009] S2: Calculate the baseline correction function according to the target displacement to obtain the corrected baseline drift function;
[0010] S3: Recalibrate the original acceleration record.
[0011] In one or more possible embodiments, step S1 includes:
[0012] S11: Perform initial baseline correction based on Iwan's two-stage baseline correction principle to obtain a m 、a f , T1, T2, D p and the corrected displacement time course, thus obtaining the initial two-stage baseline drift function uniquely; where a m 、a f Denote the average displacement of the permanent displacement segment and the strong earthquake segment, T1 and T2 denote the starting time of the strong earthquake segment and the permanent displacement segment, respectively. p Indicates permanent displacement;
[0013] S12: Fit the displacement time history after initial correction by the displacement ramp function to obtain t0 and t P , where t0 represents the starting offset time of permanent displacement, t P Indicates the moment when the permanent displacement is first reached.
[0014] In one or more possible embodiments, in step S11, the baseline correction is performed on the original data based on the Iwan two-stage baseline correction principle to obtain the average drift a of the strong earthquake segment and the permanent offset segment. m and a f and its corresponding drift starting time T1 and T2, and then the corrected acceleration record is obtained by removing the initial baseline drift in the original acceleration record, and then the corrected displacement time history is obtained by performing two integrations.
[0015] In one or more possible embodiments, step S2 includes:
[0016] S21: Introduce the pulse function as the baseline correction function, and set the baseline correction function starting time t1 = t0 and ending time t2 = t P , and then by the target displacement D target With permanent displacement D p The pulse amplitude A is calculated by the difference between them, that is, the three key parameters t1, t2 and A of the baseline correction function have been determined, thereby obtaining a unique baseline correction function;
[0017] S22: Correcting the initial two-segment baseline drift function using the obtained baseline correction function to obtain a corrected baseline drift function.
[0018] In one or more possible embodiments, in step S21, the pulse amplitude A is calculated as:
[0019] In one or more possible embodiments, step S3 includes: removing the corrected baseline drift from the original acceleration record to obtain a corrected acceleration time history, and then integrating to obtain a corrected velocity and displacement time history, and the correction is completed.
[0020] On the other hand, the present application also discloses:
[0021] A near-fault ground motion baseline correction device, the near-fault ground motion baseline correction device comprising a processor and a memory;
[0022] The memory is used to store computer programs or instructions.
[0023] The processor is used to execute the computer program or instructions in the memory, and when the computer program or instructions are executed, it is used to implement the above-mentioned near-fault ground motion baseline correction method.
[0024] On the other hand, the present application also discloses:
[0025] A computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the computer program implements the aforementioned near-fault ground motion baseline correction method.
[0026] The aforementioned main scheme of the present application and its further options can be freely combined to form multiple schemes, all of which are schemes that can be adopted and claimed for protection in the present application. After understanding the scheme of the present application, those skilled in the art can understand that there are multiple combinations based on the prior art and common knowledge, all of which are technical schemes to be protected by the present application, and they are not exhaustively listed here.
[0027] Beneficial effects of this application:
[0028] Based on the effective restoration of the near-fault seismic permanent displacement, the present application can accurately match the preset target permanent displacement, so that the permanent displacement recorded near the fault after correction accurately matches the GPS co-seismic displacement or different target permanent displacement levels. BRIEF DESCRIPTION OF THE DRAWINGS
[0029] Figure 1 It is a schematic diagram of the process flow of this application method;
[0030] Figure 2 It is a schematic diagram of Iwan's two-stage initial baseline drift function;
[0031] Figure 3 is a schematic diagram of the ramp function;
[0032] Figure 4 is a schematic diagram of the baseline correction function;
[0033] Figure 5 is a schematic diagram of the baseline drift function after correction;
[0034] Figure 6 It is the correction result after matching the GPS co-seismic displacement in the EW direction of TCU052 station;
[0035] Figure 7 It is the correction result after matching the GPS co-seismic displacement in the EW direction of TCU078 station;
[0036] Figure 8 It is the correction result after matching the GPS co-seismic displacement in the NS direction of TCU137 station;
[0037] Fig. 9 It is the correction result after matching the GPS co-seismic displacement in the EW direction of TCU052 station;
[0038] Fig.10 It is a comparison of different permanent displacement levels recorded in the EW direction at station TCU052. DETAILED DESCRIPTION
[0039] The following describes the embodiments of the present application through specific examples, and those skilled in the art can easily understand other advantages and effects of the present application from the contents disclosed in this specification. The present application can also be implemented or applied through other different specific embodiments, and the details in this specification can also be modified or changed in various ways based on different viewpoints and applications without departing from the spirit of the present application. It should be noted that the following embodiments and features in the embodiments can be combined with each other without conflict.
[0040] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, further definition and explanation thereof is not required in subsequent drawings.
[0041] In the description of the present application, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc. indicate positions or positional relationships based on the positions or positional relationships shown in the accompanying drawings, or the positions or positional relationships in which the product of the application is usually placed when in use. They are only for the convenience of describing the present application and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific position, be constructed and operated in a specific position, and therefore cannot be understood as a limitation on the present application. In addition, the terms "first", "second", "third", etc. are only used to distinguish the description, and cannot be understood as indicating or implying relative importance.
[0042] In addition, the terms "horizontal", "vertical", "overhanging" and the like do not mean that the components are required to be absolutely horizontal or overhanging, but can be slightly tilted. For example, "horizontal" only means that its direction is more horizontal than "vertical", and does not mean that the structure must be completely horizontal, but can be slightly tilted.
[0043] In the description of this application, it should also be noted that, unless otherwise clearly specified and limited, the terms "set", "install", "connect", and "connect" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be the internal communication of two elements. For ordinary technicians in this field, the specific meanings of the above terms in this application can be understood according to specific circumstances.
[0044] In addition, the present application would like to point out that, in the present application, unless the specific structure, connection relationship, positional relationship, power source relationship, etc. are specifically written out, the structure, connection relationship, positional relationship, power source relationship, etc. involved in the present application are all known by those skilled in the art on the basis of the prior art without creative work.
[0045] refer to Figure 1 As shown in the figure, a near-fault seismic baseline correction method for matching target permanent displacement is shown, and the near-fault seismic baseline correction method includes the following steps.
[0046] Step S1: Perform initial baseline correction and fit the corrected displacement time course.
[0047] In one or more possible embodiments, step S1 includes:
[0048] Step S11: Perform initial baseline correction based on Iwan's two-stage baseline correction principle to obtain a m 、a f , T1, T2, D p and the corrected displacement time course, thus uniquely obtaining the initial two-stage baseline drift function, such as Figure 2 As shown. m 、a f Denote the average displacement of the permanent displacement segment and the strong earthquake segment, T1 and T2 denote the starting time of the strong earthquake segment and the permanent displacement segment, respectively. p Indicates permanent displacement.
[0049] Preferably, in step S11, the baseline correction is performed on the original data based on the Iwan two-stage baseline correction principle to obtain the average drift a of the strong earthquake segment and the permanent offset segment. m and a fand its corresponding drift starting time T1 and T2, and then the corrected acceleration record is obtained by removing the initial baseline drift in the original acceleration record, and then two integrations are performed to obtain the corrected displacement time history, which is called the Iwan initial correction result.
[0050] Figure 2 Middle f represents the average offset of the permanent offset segment, which can be obtained by linearly fitting the permanent offset segment of the velocity time history, as shown in formula (1), m It represents the average displacement of the strong earthquake segment and can be obtained by formula (2).
[0051] V f (t) = V0 + a f t, t>T2 (1)
[0052]
[0053] Where T1 and T2 represent the starting time of the strong earthquake segment and the permanent displacement segment, respectively, and V0 is equal to the velocity value corresponding to the time T2.
[0054] Step S12: Fit the displacement time history after initial correction by using the displacement ramp function to obtain t0 and t P , where t0 represents the starting offset time of permanent displacement, t P Indicates the moment when the permanent displacement is first reached.
[0055] Specifically, the introduction Figure 3 The ramp function shown in FIG. 1 is used to describe the overall change trend of the permanent displacement, that is, the Iwan initial correction result obtained in step S12 is optimally fitted by the ramp function to obtain the key parameter values t0 and t P t0~t P The segment describes the rising segment shape of the displacement time history.
[0056] Step S2: Calculate the baseline correction function according to the target displacement to obtain a corrected baseline drift function.
[0057] Preferably, step S2 comprises:
[0058] Step S21: Introduce the pulse function as the baseline correction function, and set the baseline correction function to start at t1 = t0 and end at t2 = t P , and then by the target displacement D target With permanent displacement D p The pulse amplitude A is calculated by the difference between them, that is, the three key parameters t1, t2 and A of the baseline correction function have been determined, thereby obtaining a unique baseline correction function.
[0059] Specifically, the pulse function is introduced as the baseline correction function, such as Figure 4 In order to ensure the continuity and rationality of the displacement time history after the final correction, it is assumed Figure 4 The effective interval t1 to t2 of the baseline correction function shown is exactly the same as Figure 3 The rising section t0~t P coincidence, that is, let t1 = t0, t2 = t P , that is to say Figure 2 The correction of the baseline drift function shown only occurs in the rising phase of the displacement time course and does not affect the displacement phase before t0 and t P The overall shape of the post-permanent displacement section. Then:
[0060]
[0061] Where t1 represents the pulse start time, t2 represents the pulse end time, and A represents the pulse amplitude.
[0062] Calculate the target displacement D target With permanent displacement D p The difference ΔD between them is shown in equation (6); then, let ΔD be Figure 4 The displacement peaks of the baseline correction function are equal, as shown in equation (7), and finally the pulse amplitude A is obtained, as shown in equation (8). So far, the three key parameters t1, t2 and A of the baseline correction function have been determined, thereby obtaining a unique baseline correction function.
[0063] ΔD=D target -D p (6)
[0064]
[0065] Step S22: Correct the initial two-stage baseline drift function by using the obtained baseline correction function to obtain a corrected baseline drift function. Figure 5 As shown, Figure 5 T1, T2 and t1, t2 in are calculated by step S11 and step S21 respectively, and there is no definite size relationship between the two. Figure 5 Just an example.
[0066] Step S3: Recalibrate the original acceleration record. Specifically, step S3 includes: removing the corrected baseline drift from the original acceleration record to obtain the corrected acceleration time history, and then integrating to obtain the corrected velocity and displacement time history, and the calibration is completed.
[0067] Based on the same technical concept, the present application also discloses a near-fault seismic baseline correction device, which includes a processor and a memory, and optionally, also includes: a communication interface and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus.
[0068] A computer program is stored in the memory, and when the program is executed by the processor, the processor is enabled to complete the steps of the above-mentioned near-fault ground motion baseline correction method.
[0069] The communication bus mentioned in the above-mentioned near-fault ground motion baseline correction device can be a Peripheral Component Inter connect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. The communication bus can be divided into an address bus, a data bus, a control bus, etc. The communication interface is used for communication between the above-mentioned beam-column node design device and other devices.
[0070] The memory may include a random access memory (RAM) or a non-volatile memory (NVM), such as at least one disk memory. Optionally, the memory may also be at least one storage device located away from the aforementioned processor.
[0071] The above-mentioned processor can be a general-purpose processor, including a central processing unit, a network processor (Network Processor, NP), etc.; it can also be a digital signal processing processor (Digital Signal Processing, DSP), an application-specific integrated circuit, a field programmable gate array or other programmable logic device, a discrete gate or transistor logic device, a discrete hardware component, etc.
[0072] Based on the same technical concept, the present application also discloses: a computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps of the near-fault ground motion baseline correction method as described above are implemented.
[0073] The above-mentioned computer-readable storage medium can be any available medium or data storage device that can be accessed by the processor in the near-fault ground motion baseline correction device, including but not limited to magnetic storage such as floppy disks, hard disks, magnetic tapes, magneto-optical disks (MO), etc., optical storage such as CDs, DVDs, BDs, HVDs, etc., and semiconductor storage such as ROMs, EPROMs, EEPROMs, non-volatile memories (NANDFLASH), solid-state drives (SSDs), etc.
[0074] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.
[0075] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.
[0076] These computer program instructions may also be stored in a computer readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.
[0077] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process in the computer or other programmable device. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.
[0078] Example 1
[0079] Taking a typical near-fault record in an earthquake as an example, the GPS co-seismic displacement is used as the target displacement. Figure 6 , Figure 7 , Figure 8 , Fig. 9 The correction results of the matching target permanent displacement obtained by the method of this application are given, and compared with the original records and the results obtained by the Iwan baseline correction method. It can be seen that the correction results obtained by the patented method can accurately match the target GPS co-seismic displacement without changing the main morphological characteristics of the displacement time history, which significantly increases the effectiveness and applicability of the correction results.
[0080] Example 2
[0081] Taking the EW direction records of TCU052 and TCU102 stations in the earthquake as an example, the baseline correction method of this application is used to pre-set different target displacement levels (100cm, 150cm, 200cm, 250cm, 300cm, 350cm), and the corresponding correction results can be obtained, such as Fig.10 As shown in the figure, it can be seen that the permanent displacement of the near-fault ground motion record after correction can accurately match the preset target displacement level and can better retain the main morphological characteristics of the displacement time history.
[0082] The above description is only a preferred embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present application should be included in the protection scope of the present application.
Claims
1. A near-fault ground motion baseline correction method for matching target permanent displacement, characterized in that: The near-fault ground motion baseline correction method comprises: S1: Perform initial baseline correction and fit the corrected displacement time course; S2: Calculate the baseline correction function according to the target displacement to obtain the corrected baseline drift function; S3: Recalibrate the original acceleration record.
2. The near-fault ground motion baseline correction method according to claim 1, characterized in that: Step S1 includes: S11: Perform initial baseline correction based on Iwan's two-stage baseline correction principle to obtain a m 、a f , T1, T2, D p and the corrected displacement time course, thus obtaining the initial two-stage baseline drift function uniquely; where a m 、a f Denote the average displacement of the permanent displacement segment and the strong earthquake segment, T1 and T2 denote the starting time of the strong earthquake segment and the permanent displacement segment, respectively. p Indicates permanent displacement; S12: Fit the displacement time history after initial correction by the displacement ramp function to obtain t0 and t P , where t0 represents the starting offset time of permanent displacement, t P Indicates the moment when the permanent displacement is first reached.
3. The near-fault ground motion baseline correction method according to claim 2, characterized in that: In step S11, the baseline correction of the original data is performed based on the Iwan two-stage baseline correction principle to obtain the average drift a of the strong earthquake segment and the permanent offset segment. m and a f and its corresponding drift starting time T1 and T2, and then the corrected acceleration record is obtained by removing the initial baseline drift in the original acceleration record, and then the corrected displacement time history is obtained by performing two integrations.
4. The near-fault ground motion baseline correction method according to claim 2, characterized in that: Step S2 includes: S21: Introduce the pulse function as the baseline correction function, and set the baseline correction function starting time t1 = t0 and ending time t2 = t P , and then by the target displacement D target With permanent displacement D p The pulse amplitude A is calculated by the difference between them, that is, the three key parameters t1, t2 and A of the baseline correction function have been determined, thereby obtaining a unique baseline correction function; S22: Correcting the initial two-segment baseline drift function using the obtained baseline correction function to obtain a corrected baseline drift function.
5. The near-fault ground motion baseline correction method according to claim 4, characterized in that: In step S21, the pulse amplitude A is calculated as:
6. The near-fault ground motion baseline correction method according to claim 4, characterized in that: Step S3 includes: removing the corrected baseline drift from the original acceleration record to obtain the corrected acceleration time history, and then integrating to obtain the corrected velocity and displacement time history, and the correction is completed.
7. A near-fault ground motion baseline correction device, characterized in that: The near-fault ground motion baseline correction device comprises a processor and a memory; The memory is used to store computer programs or instructions. The processor is used to execute the computer program or instructions in the memory, and when the computer program or instructions are executed, it is used to implement the near-fault ground motion baseline correction method according to any one of claims 1 to 6.
8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the near-fault ground motion baseline correction method according to any one of claims 1 to 6 is implemented.