A method and system for constructing a site-specific ground motion calculation model

Abstract

The application provides a construction method and system of a site seismic ground motion calculation model, and relates to the field of engineering seismology. The method comprises the following steps: identifying the stratification information of each borehole according to the number of the borehole; if the stratification information of the borehole does not meet the pre-designed calculation condition, re-dividing the stratification of each borehole according to the maximum stratification thickness; calculating the equivalent shear wave velocity corresponding to each stratification of each borehole by using a numerical equivalent shear wave velocity calculation algorithm according to the measured wave velocity data in a measured wave velocity input file; generating the borehole seismic ground motion calculation model data corresponding to each borehole based on the stratification information and the equivalent shear wave velocity corresponding to the stratification of each borehole; constructing a site seismic ground motion calculation model based on the borehole seismic ground motion calculation model data corresponding to all boreholes and combining soil dynamics test information, and generating a site seismic ground motion calculation model file. The technical scheme provided by the application improves the efficiency and accuracy of the construction of the site seismic ground motion calculation model.

Description

Technical Field

[0001] This application relates to the field of engineering earthquake technology, and in particular to a method and system for constructing a site ground motion calculation model. Background Technology

[0002] Domestic and international research and earthquake damage examples show that the overlying soil layer has a strong amplification effect on seismic motion. In the field of engineering seismology, site ground motion calculation and analysis is an important step in scientifically determining the seismic fortification parameters of construction projects. It mainly involves constructing a site ground motion calculation model based on the type of overlying soil layer, wave velocity structure, and soil dynamic properties. The wave velocity structure can be described by the equivalent shear wave velocity.

[0003] Currently, when constructing a site ground motion calculation model, technicians typically manually analyze and divide each borehole into layers based on a columnar section of the site's overburden soil, and then calculate the equivalent shear wave velocity corresponding to each layer. However, this method is inefficient, especially when the number of boreholes exceeds 10. The calculation of the equivalent shear wave velocity corresponding to each borehole layer is labor-intensive, time-consuming, and prone to errors, ultimately affecting the accuracy of the site ground motion calculation model. Summary of the Invention

[0004] This application provides a method and system for constructing a site ground motion calculation model to solve the problems of low efficiency and poor accuracy in the construction of site ground motion calculation models in the prior art.

[0005] In a first aspect, embodiments of this application provide a method for constructing a site seismic motion calculation model, including:

[0006] In response to the user's launch operation of the target application software, an operation interface is displayed, which includes the parameter items and operation buttons required for building the site ground motion calculation model;

[0007] Configure the parameter items, which include the path of the input file, maximum layer thickness, number of boreholes, borehole number prefix, layer information separator, soil type name and description separator, soil type code and name separator, and measured wave velocity data range. After the parameter configuration is completed, in response to the touch operation of the operation button, read the soil layer information input file under the path of the input file. The soil layer information input file includes: the original borehole layer input file, the soil layer type and description input file, and the measured wave velocity input file.

[0008] Based on the borehole number, the layer information of each borehole is identified from the original layer input file and the soil layer category and description input file. If the layer information of the borehole does not meet the preset calculation conditions, the layer of each borehole is re-divided according to the maximum layer thickness.

[0009] According to the measured wave velocity data in the measured wave velocity input file, the numerical equivalent shear wave velocity calculation algorithm is used to calculate the equivalent shear wave velocity corresponding to each layer of each borehole;

[0010] Based on the layer information and the equivalent shear wave velocity corresponding to each layer of each borehole, the borehole seismic vibration calculation model data corresponding to each borehole is generated, and based on the borehole seismic vibration calculation model data corresponding to all boreholes, a site seismic vibration calculation model is constructed in combination with the soil dynamics test information to generate a site seismic vibration calculation model file.

[0011] Optionally, the soil layer information input file further includes: the soil layer information input file further includes: a cover layer thickness input file, the cover layer thickness input file includes the cover layer thickness of each borehole, and the cover layer thickness is the cumulative value of all target layer thicknesses; the measured wave velocity input file includes the measured wave velocity data of all boreholes, and the measured wave velocity data of each borehole includes the measured shear wave velocity of all equidistant layers, and the thickness of the equidistant layers is much greater than the fixed thickness used in the numerical equivalent shear wave velocity calculation;

[0012] According to the measured wave velocity data in the measured wave velocity input file, the numerical equivalent shear wave velocity calculation algorithm is used to calculate the equivalent shear wave velocity corresponding to each layer of each borehole, including:

[0013] Each layer of each borehole is taken as a target layer respectively;

[0014] For each target layer, according to the measured shear wave velocity of the equidistant layers, the numerical equivalent shear wave velocity calculation algorithm is used to discretize the target layer into unit layers with a fixed thickness, according to the shear wave velocity corresponding to each unit layer, the travel time of each unit layer is calculated, the travel times of all unit layers belonging to the same target layer are accumulated to obtain the target layer travel time, and the target layer thickness is divided by the target layer travel time to obtain the equivalent shear wave velocity corresponding to each layer of each borehole.

[0015] Optionally, after obtaining the equivalent shear wave velocity corresponding to each layer of each borehole, the method further includes:

[0016] According to the equivalent shear wave velocity corresponding to each layer of each borehole, the equivalent wave velocity data of each borehole is generated;

[0017] The equivalent wave velocity data of each borehole is compared with the measured wave velocity data to generate a wave velocity comparison ladder diagram, and the wave velocity comparison ladder diagram is used to verify the accuracy of the equivalent shear wave velocity corresponding to all layers of the borehole.

[0018] Optionally, the redivision of the layers of each borehole according to the maximum layer thickness includes:

[0019] selecting other factors including at least one of the following: soil physical properties, soil dynamic properties, and soil categories;

[0020] repartitioning the layers of each borehole based on the maximum layer thickness and the other factors.

[0021] Optionally, after the borehole seismic ground motion calculation model data corresponding to each borehole is generated, the method further includes:

[0022] According to the number of boreholes and the numbers of the boreholes, a unique borehole folder is established for each borehole, wherein the borehole seismic ground motion calculation model data corresponding to different boreholes are stored in different borehole folders, and the borehole seismic ground motion calculation model data of the same borehole are stored in the same file in the same borehole folder.

[0023] Optionally, in the process of configuring the parameter items, the method further includes:

[0024] determining whether the input content of the parameter items meets the corresponding configuration conditions, and if not, displaying a prompt box in a specified area of the operation interface to prompt the user to check and modify the input content, so as to prevent the task from failing due to misoperation.

[0025] Optionally, the number of operation buttons is multiple, each operation button corresponds to one link in the site seismic ground motion calculation model construction process, and the operation interface includes a progress state display mark corresponding to each operation button.

[0026] In a second aspect, an embodiment of the present application provides a site seismic ground motion calculation model construction system, including:

[0027] The display module is configured to display an operation interface in response to a user starting operation on a target application software, wherein the operation interface includes parameter items and operation buttons required for site seismic ground motion calculation model construction.

[0028] The configuration reading module is configured to configure the parameter items, wherein the parameter items include a path of an input file, a maximum layer thickness, a number of boreholes, a hole number prefix, a layer information separator, a soil category name and description separator, a soil category code and name separator, and a measured wave velocity data range; and after the parameter configuration is completed, the configuration reading module is configured to read a soil layer information input file under the path of the input file in response to a touch operation of the operation button, wherein the soil layer information input file includes a borehole original layer input file, a soil layer category and description input file, and a measured wave velocity input file.

[0029] The identification division module is configured to identify the layering information of each borehole from the borehole original layering input file and the soil layer category and description input file according to the number of the borehole, and to redivide the layering of each borehole according to the maximum layering thickness if the layering information of the borehole does not satisfy a pre-designed calculation condition;

[0030] The calculation module is configured to calculate the equivalent shear wave velocity corresponding to each layer of each borehole by using a numerical equivalent shear wave velocity calculation algorithm according to the measured wave velocity data in the measured wave velocity input file.

[0031] The generation and construction module is configured to generate borehole seismic vibration calculation model data corresponding to each borehole based on the layering information and the equivalent shear wave velocity corresponding to each layer of each borehole, to construct a site seismic vibration calculation model based on the borehole seismic vibration calculation model data corresponding to all boreholes and in combination with soil dynamics test information, and to generate a site seismic vibration calculation model file.

[0032] In a third aspect, an embodiment of the present application provides a computing device, including a processing component and a storage component; the storage component stores one or more computer instructions; the one or more computer instructions are used to be called and executed by the processing component, and implement the method for constructing a site seismic vibration calculation model according to any one of the first aspect.

[0033] In a fourth aspect, an embodiment of the present application provides a computer storage medium, which stores a computer program; when the computer program is executed by a computer, the method for constructing a site seismic vibration calculation model according to any one of the first aspect is implemented.

[0034] In the embodiments of the present application, a construction method of a site seismic ground motion calculation model is provided, which comprises: in response to a user starting operation on a target application software, displaying an operation interface, the operation interface comprising parameter items and operation buttons required for construction of the site seismic ground motion calculation model; configuring the parameter items, the parameter items comprising a path of an input file, a maximum layering thickness, a number of drill holes, a hole number prefix, a layering information separator, a soil class name and description separator, a soil class code and name separator, and a measured wave velocity data range; after the parameter configuration is completed, in response to a touch operation of the operation button, reading a soil layer information input file in the path of the input file, the soil layer information input file comprising: a drill hole original layering input file, a soil layer category and description input file, and a measured wave velocity input file; according to the number of the drill holes, identifying the layering information of each drill hole from the drill hole original layering input file and the soil layer category and description input file, and if the layering information of the drill hole does not meet a pre-designed calculation condition, re-dividing the layering of each drill hole according to the maximum layering thickness; according to the measured wave velocity data in the measured wave velocity input file, using a numerical equivalent shear wave velocity calculation algorithm to calculate the equivalent shear wave velocity corresponding to each layering of each drill hole; based on the layering information of each drill hole and the equivalent shear wave velocity corresponding to the layering, generating drill hole seismic ground motion calculation model data corresponding to each drill hole, and based on the drill hole seismic ground motion calculation model data corresponding to all the drill holes, constructing a site seismic ground motion calculation model in combination with soil dynamics test information to generate a site seismic ground motion calculation model file.

[0035] The embodiments of the present application automatically generate an operation interface in response to a user starting operation, and provide a parameter configuration function, thereby reducing the need for manual intervention. The process of reading and analyzing the soil layer information input file realizes automation, greatly shortening the time from data preparation to model generation. The user only needs to click the corresponding operation button to complete a complex data processing task, such as reading a file, modeling, etc., thereby simplifying the operation process and improving the work efficiency. The soil layer is re-divided based on the maximum layering thickness, which ensures that each layering meets the engineering requirements and accurately reflects the geological characteristics, thereby avoiding simulation errors caused by unreasonable layering. The numerical equivalent shear wave velocity calculation algorithm is used to discretize the target layer into unit layers with a fixed thickness, the total travel time is calculated by accumulating the travel time, and the equivalent shear wave velocity is finally obtained. This method not only is easy to program, but also improves the calculation accuracy and reduces human errors. Therefore, the embodiments of the present application provide an efficient, accurate, easy-to-use and high-quality site seismic ground motion calculation model construction method, which can improve the work efficiency and result reliability of site seismic ground motion analysis, and provides technical support for building seismic design.

[0036] These and other aspects of the present application will become more apparent from the following description of the embodiments. BRIEF DESCRIPTION OF DRAWINGS

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

[0038] Figure 1 A flowchart of a construction method of a site seismic ground motion calculation model provided by an embodiment of the present application is shown in the figure.

[0039] Figure 2 A schematic diagram of an operation interface provided by an embodiment of the present application is shown in the figure.

[0040] Figure 3 A data format schematic diagram of a borehole original layering input file provided by an embodiment of the present application is shown in the figure.

[0041] Figure 4 A data format schematic diagram of a soil layer category and description input file provided by an embodiment of the present application is shown in the figure.

[0042] Figure 5 A data format schematic diagram of a measured wave velocity input file provided by an embodiment of the present application is shown in the figure.

[0043] Figure 6 A data format schematic diagram of a cover layer thickness input file provided by an embodiment of the present application is shown in the figure.

[0044] Figure 7 A reading flowchart of a soil layer information input file provided by an embodiment of the present application is shown in the figure.

[0045] Figure 8 An automatic construction flowchart of a borehole seismic ground motion calculation model provided by an embodiment of the present application is shown in the figure.

[0046] Figure 9 A schematic diagram of an error prompt corresponding to a modeling stage provided by an embodiment of the present application is shown in the figure.

[0047] Figure 10 A schematic diagram of a scenario for solving equivalent shear wave velocity is shown in the figure.

[0048] Figure 11 A flowchart of solving equivalent shear wave velocity provided by an embodiment of the present application is shown in the figure.

[0049] Figure 12 A schematic diagram of an error prompt corresponding to an equivalent shear wave velocity calculation stage provided by an embodiment of the present application is shown in the figure.

[0050] Figure 13 A schematic diagram of a wave velocity comparison ladder diagram provided by an embodiment of the present application is shown in the figure.

[0051] Figure 14 This is a flowchart illustrating the software shutdown step provided in an embodiment of this application.

[0052] Figure 15 A schematic diagram of the structure of a system for constructing a site seismic motion calculation model provided in an embodiment of this application;

[0053] Figure 16 This is a schematic diagram of the structure of a computing device provided in an embodiment of this application. Detailed Implementation

[0054] To enable those skilled in the art to better understand the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings.

[0055] In some of the processes described in the specification, claims, and accompanying drawings of this application, multiple operations appearing in a specific order are included. However, it should be clearly understood that these operations may not be executed in the order they appear herein, or may be executed in parallel. The operation numbers, such as 11, 12, etc., are merely used to distinguish different operations and do not themselves represent any execution order. Furthermore, these processes may include more or fewer operations, and these operations may be executed sequentially or in parallel. It should be noted that the descriptions such as "first," "second," etc., in this document are used to distinguish different messages, devices, modules, etc., and do not represent a sequential order, nor do they limit "first" and "second" to different types.

[0056] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0057] To improve the efficiency and quality of constructing site ground motion calculation models, this application proposes a method for constructing such models. For example, Figure 1 A flowchart illustrating a method for constructing a site seismic motion calculation model provided in this application embodiment is shown below. Figure 1 As shown, the method includes:

[0058] S11. In response to the user's startup operation on the target application software, the operation interface is displayed, which includes the parameter items and operation buttons required for building the site ground motion calculation model.

[0059] It should be understood that the user is also called the user. The operation interface is also called the drawing operation interface. The site seismic ground motion calculation model construction system is also called the site seismic ground motion calculation model automatic creation system.

[0060] Figure 2 The schematic diagram of the operation interface provided by the embodiment of the present application. The user opens the site seismic ground motion calculation model automatic creation system to automatically generate the operation interface. The operation interface displays all the parameters and operation buttons required for the site seismic ground motion calculation model automatic creation. As shown in the figure, the input parameters include: the path of the input file, the number of boreholes, the hole number prefix, the maximum layer interval (meters), the layer information separator, the soil class name and description separator, the soil class code and name separator, the test wave velocity data range (i.e. the measured wave velocity data range), the total number of modeling soil classes, the soil layer model table, and the soil layer model data range. The operation buttons include reading the soil layer information input file, modeling, creating a borehole folder, creating a model calculation file, and closing the software. The above maximum layer interval is also called the maximum layer thickness. Figure 2

[0061] Optionally, the number of operation buttons is multiple, each operation button corresponds to a link in the site seismic ground motion calculation model construction process, and the operation interface includes a progress state display mark corresponding to each operation button. As an example, as shown in the figure, the progress state display mark can be a work stage prompt control. The work stage prompt control is initially defaulted to a color (such as red, gray, etc.) for indicating an unprocessed state. Different states are indicated by color changes. As an example, when the work stage prompt control is yellow, it is used to indicate that the processing is in progress. When the work stage prompt control is green, it is used to indicate that the processing is complete or ends. As can be seen, the embodiment of the present application has the advantage of clear operation process and can display the current stage in real time. The site seismic ground motion calculation model construction system divides the modeling into the work stages of “reading the soil layer information input file”, “modeling”, “creating a borehole folder”, “creating a model calculation file”, and “closing the software” according to the soil layer modeling practical experience. When using the site seismic ground motion calculation model construction system, these buttons can be clicked in sequence to complete the modeling. Each stage has a completion prompt, for example, the stage is displayed in red if it is not completed, and in green if it is completed, so as to avoid the user forgetting the current work stage due to being temporarily busy with other things. Figure 2

[0062] ​​S12, configure parameter items, the parameter items include the path where the input file is located, the maximum layered thickness, the number of drill holes, the hole number prefix, the layered information separator, the soil class name and description separator, the soil class code and name separator, and the measured wave velocity data range. After the parameter configuration is completed, the soil layer information input file under the path where the input file is located is read in response to the touch operation of the operation button. The soil layer information input file includes: a drill hole original layered input file, a soil layer category and description input file, and a measured wave velocity input file.

[0063] In the embodiment, the soil layer information input file further includes: a cover layer thickness input file, etc. In the embodiment, the parameter items or parameters can be displayed in the parameter input box in the form of default values. The user can change them according to his own needs. For example, the parameter configuration can be performed through the parameter input box in Figure 1 Therefore, the embodiment has the following advantages: strong man-machine interaction, more convenient to use. The default parameters are directly displayed in the input box of the operation interface, the user can clearly obtain them, and can adaptively modify them according to the needs.

[0064] It should be understood that the format of the soil layer information input file is consistent with the format of the columnar chart information file in the existing survey software. For example, the data format of the drill hole original layered input file is as shown in Figure 3 Each row is used to store the soil class code, the layer bottom depth, etc. For the first row of each drill hole, the corresponding drill hole number (or hole number, number, etc.) is additionally stored. The data format of the soil layer category and description input file is as shown in Figure 4 Each row is used to store the soil class code, the soil class name, the soil class description, etc. Optionally, to enhance the orderliness, the soil layer wave velocity can be stored in an electronic form, and the test wave velocity data range is filled in the operation interface to be used. The above measured wave velocity input file or test wave velocity input file. The data format of the measured wave velocity input file is as shown in Figure 5As shown, the measured wave velocity input file includes the measured wave velocity data of all the boreholes, and the measured wave velocity data of each borehole includes the measured shear wave velocity of all the equally-spaced layers, the measured shear wave velocity or the measured wave velocity data, the test wave velocity, the test shear wave velocity, the test wave velocity data, etc., and the thickness of the equally-spaced layers is 1 m. Specifically, the measured wave velocity data (or the shear wave velocity value) of one borehole is arranged as a column, and arranged from top to bottom as the depth increases (generally, the depth increases by 1 m), and the measured wave velocity data of all the boreholes is arranged from left to right. The overburden thickness input file can be a data (Data, DAT) file, and the DAT file is a common file format, which is usually used to store various types of data. The extension ".dat" indicates that the file contains raw data. For example, the overburden thickness information of all the boreholes is arranged in the data file from left to right, as shown in the following table. Figure 6 As shown, the overburden thicknesses of the four boreholes are 13, 19, 12 and 19 respectively.

[0065] For example, Figure 7 The reading process of the soil layer information input file provided by the embodiment of the present application is shown in the following table. Figure 7 As shown, the reading process includes: S71, clicking the button of "reading the soil layer information input file". S72, reading the borehole original layer input file, the soil layer category and description input file, the measured wave velocity input file and the overburden thickness input file. S73, judging whether the borehole numbers are arranged in order, if yes, executing S74, otherwise executing S76. S74, judging whether the depths of the layers are arranged in order, if yes, executing S75, otherwise executing S76. S75, forming the borehole soil property array, the layer bottom depth array, the measured wave velocity array and the position array. S76, ending the "reading of the soil layer information input file".

[0066] Optionally, in the process of configuring the parameter item, the method further includes: judging whether the input content of the parameter item meets the corresponding configuration condition, if not, displaying a prompt box in the specified area of the operation interface to prompt the user to check and modify the input content, so as to prevent the task from failing due to misoperation.

[0067] S13, according to the number of the borehole, identifying the layer information of each borehole from the borehole original layer input file and the soil layer category and description input file, and if the layer information of the borehole does not meet the pre-designed calculation condition, re-dividing the layers of each borehole according to the maximum layer thickness. The number of the borehole is increased by default from 1.

[0068] S14, according to the measured wave velocity data in the measured wave velocity input file, calculating the equivalent shear wave velocity corresponding to each layer of each borehole by using the numerical equivalent shear wave velocity calculation algorithm.

[0069] S15. Based on the layer information of each borehole and the equivalent shear wave velocity corresponding to each layer, generate borehole ground motion calculation model data for each borehole. Based on the borehole ground motion calculation model data for all boreholes, and combined with soil dynamics test information, construct a site ground motion calculation model and generate a site ground motion calculation model file.

[0070] As described in S13-S15, the main idea of ​​automatically constructing a borehole seismic motion calculation model is to re-layer the soil layers in the borehole based on the input soil properties, bottom layer depth, and test wave velocity, according to the set maximum layer thickness, to establish layers suitable for calculation, generate a fine-layered borehole model matrix, then calculate the equivalent shear wave velocity corresponding to each layer of the borehole based on the test wave velocity, generate borehole seismic motion calculation model data, and output the borehole seismic motion calculation model data file. For example, the automated construction process of the borehole seismic motion calculation model is as follows: Figure 8 As shown, the process includes: S81, generating an initial borehole model matrix based on input information. S82, re-stratifying the borehole based on factors such as maximum layer thickness. S83, generating a refined layered borehole model matrix. S84, determining whether the layer thickness is uniform; if so, proceed to S86; otherwise, proceed to S85. S85, adjusting the layers to regenerate the refined layered borehole model matrix. S86, calculating the equivalent shear wave velocity corresponding to each layer. S87, generating a layered borehole model matrix with wave velocities. S88, performing bedrock layer processing. S89, generating the final layered borehole model matrix. S90, outputting the borehole seismic calculation model data file.

[0071] Afterwards, users can assign soil type numbers and other data to each soil type in a spreadsheet based on the processed dynamic triaxial test data (hereinafter referred to as dynamic triaxial data), forming a complete seismic motion calculation model spreadsheet file. Therefore, the embodiments of this application have the advantage of being closely aligned with actual working conditions. Considering that the selection, classification, and merging of dynamic triaxial data required for site seismic motion calculation modeling requires the experience of technical personnel, it is relatively flexible. If calculated directly by computer, it is very likely that it will not be completely consistent with human expectations, and the time spent modifying the model may exceed the time required for direct manual processing. Therefore, the selection, classification, and merging of dynamic triaxial data is handled by professional technical personnel. The operation is performed in the spreadsheet, assigning soil type information to the soil layers.

[0072] To prevent program errors caused by incorrect input parameters (which could refer to the parameters mentioned above), configure a try...catch statement to promptly remind the user to check and modify the input parameters, such as... Figure 9 As shown.

[0073] By executing steps S11-S15, the columnar chart information input file of the overburden soil layer can be obtained, and the columnar chart information of the overburden soil layer of the site can be read from the columnar chart information input file in one key. The format of the columnar chart information input file is consistent with the format of the existing survey software. After the overburden soil layer survey of the site is completed, the columnar chart information file thereof can be used as the input file of the construction system of the site seismic vibration calculation model, without the need for further editing, thereby reducing the repetitive work. All the information can be read by clicking the corresponding button with the mouse. Similarly, the soil layer information input file of the overburden soil layer is obtained, and the stratification information of each borehole is automatically identified from the soil layer information input file. According to the soil layer information input file, the stratification information of each borehole is identified according to the hole number, which lays a foundation for re-dividing the stratification of each borehole and finding the corresponding soil layer wave velocity and calculating the equivalent shear wave velocity.

[0074] The embodiment of the present application automatically generates an operation interface in response to a user starting operation, and provides a parameter configuration function, thereby reducing the need for manual intervention. The process of reading and analyzing the soil layer information input file realizes automation, and greatly shortens the time from data preparation to model generation. The user can complete complex data processing tasks such as reading files and modeling by clicking the corresponding operation button, thereby simplifying the operation process and improving the work efficiency. The soil layer is re-divided based on the maximum stratification thickness, which ensures that each stratification meets the engineering requirements and accurately reflects the geological characteristics, and avoids simulation errors caused by unreasonable stratification. The numerical equivalent shear wave velocity calculation algorithm is adopted to discretize the target layer into unit layers with a fixed thickness (for example, 1E-6m), the total travel time is calculated by accumulating the travel time, and the equivalent shear wave velocity is finally obtained. This method not only facilitates programming implementation, but also improves the calculation accuracy and reduces human errors. Therefore, the embodiment of the present application provides an efficient, accurate, easy-to-use and high-quality site seismic vibration calculation model construction method, which can improve the work efficiency and result reliability of site seismic vibration analysis, and provides technical support for building seismic design.

[0075] In some possible embodiments, S13, re-dividing the stratification of each borehole based on the maximum stratification thickness, comprises:

[0076] Step 131, selecting other factors, the other factors comprising at least one of: soil physical properties, soil dynamic properties, and soil categories.

[0077] Step 132, re-dividing the stratification of each borehole based on the maximum stratification thickness and the other factors.

[0078] By executing steps 131-132, the embodiment of the present application can realize automatic division of the soil layer. As understood, the traditional columnar chart information input file of the overburden soil layer is manually edited to obtain the stratification information of each borehole, which is time-consuming and laborious. The embodiment of the present application can realize automatic division of the soil layer, thereby improving the work efficiency. Figure 1Generally, the soil layers are divided according to the soil types at different positions, and the traditional layering method generally does not meet the layering requirements of the site seismic motion calculation model. Therefore, the embodiments of the application can re-layer according to the soil properties and layering thickness requirements. The construction system of the site seismic motion calculation model can automatically re-divide the soil layers for each input drill hole.

[0079] Figure 10 For the scene diagram of solving the equivalent shear wave velocity, the common equivalent shear wave velocity calculation problem is solved. The corresponding traditional equivalent shear wave velocity calculation formula is shown in formula (1)~(3).

[0080] (1)

[0081] (2)

[0082] (3)

[0083] wherein, , and are used to represent the equivalent shear wave velocities of the layer 1, the layer 2 and the layer 3 of a certain drill hole obtained by the traditional equivalent shear wave velocity calculation formula, the layering thickness of the layer 1 is , the unit is m, the travel time of the layer 1 is , the layering thickness of the layer 2 is , the unit is m, the travel time of the layer 2 is , the layering thickness of the layer 3 is , the unit is m, and the travel time of the layer 3 is . Wherein, , , , and are the test wave velocities in the equal-interval layer 1, the equal-interval layer 2, the equal-interval layer 3, the equal-interval layer 4 and the equal-interval layer 5, is the distance difference between the lower line of the layer 1 of the to-be-solved equivalent shear wave velocity model and the lower line of the equal-interval layer 1, is the distance difference between the lower line of the layer 2 of the to-be-solved equivalent shear wave velocity model and the lower line of the equal-interval layer 3, is the distance difference between the lower line of the layer 3 of the to-be-solved equivalent shear wave velocity model and the lower line of the equal-interval layer 4.

[0084] From formula (1)~(3), to calculate the equivalent shear wave velocity, the layer thickness and travel time need to be calculated respectively according to different stratification, and then divided, and the complexity mainly reflects on the calculation process of travel time. For different layer division, the difference of the formula to be solved is great, so it is relatively complex to solve this problem by using the conventional programming idea. In order to simplify this problem, the embodiment of the application provides a numerical equivalent shear wave velocity calculation algorithm. The main idea of the algorithm is to discretize all soil layers of the borehole into a series of thin layers (or unit layers), and according to the test wave velocity, the shear wave travel time of each thin layer is calculated in advance, so that the calculation of travel time is simplified to the sum of the travel time of all thin layers corresponding to the layer to be solved, and the equivalent shear wave velocity can be obtained by directly dividing the layer thickness. The more thin layers, the more accurate the numerical value obtained. Considering the improvement of computer performance, each thin layer reaches 1E-6m to obtain sufficient calculation accuracy. Therefore, the soil layer information input file also includes: a cover layer thickness input file, the cover layer thickness input file includes the cover layer thickness of each borehole, and the cover layer thickness is the cumulative value of all target layer thicknesses, which is used to determine the upper and lower limits of all target layer thicknesses; the measured wave velocity input file includes the measured shear wave velocity of equal interval stratification, and the thickness of equal interval stratification is much larger (for example, 1m) than the fixed thickness (for example, 1E-6m), and much larger means that the order of magnitude of the thickness of equal interval stratification is different from the order of magnitude of the fixed thickness, and there is a significant difference between the two orders of magnitude; step S14, according to the measured wave velocity data in the measured wave velocity input file, using the numerical equivalent shear wave velocity calculation algorithm, calculating the equivalent shear wave velocity corresponding to each layer of each borehole, including:

[0085] Step 141, each layer of each borehole is taken as a target layer respectively.

[0086] Step 142, for each target layer, according to the measured shear wave velocity of equal interval stratification, using the numerical equivalent shear wave velocity calculation algorithm, discretizing the target layer into unit layers with fixed thickness, calculating the travel time of each unit layer according to the shear wave velocity corresponding to each unit layer, accumulating the travel time of all unit layers belonging to the same target layer to obtain the target layer travel time, and dividing the target layer thickness by the target layer travel time to obtain the equivalent shear wave velocity corresponding to each layer of each borehole.

[0087] Figure 11 The flowchart for solving the equivalent shear wave velocity provided by the embodiment of the application is shown in FIG. 1. Figure 11As shown, the method comprises the following steps: S111, reading the buried depth array of each layer to be calculated. S112, reading the test wave velocity array. S113, discretizing the test wave velocity into thin layers. S114, calculating the travel time array of the discrete thin layers. S115, determining whether i is less than or equal to the number of layers to be calculated. If yes, executing S116, otherwise executing S121. S116, calculating the thickness of the i th layer. S117, finding the travel time array index corresponding to the upper and lower buried depths of the i th layer. S118, calculating the sum of the elements in the travel time array index range to obtain the travel time of the i th layer. S119, dividing the travel time by the thickness of the i th layer to obtain the equivalent wave velocity of the i th layer. S120, i = i + 1, and continue to execute S115.

[0088] Figure 12 The schematic diagram of the error prompt corresponding to the equivalent shear wave velocity calculation stage provided by the embodiment of the present application. To avoid program errors caused by user input parameter errors, a try...catch statement is configured to remind the user to check and modify the input parameters in time.

[0089] By executing steps 141-142, the numerical equivalent shear wave velocity calculation algorithm provided by the embodiment of the present application is easy to program and implement, and is not prone to errors.

[0090] In some possible embodiments, after obtaining the equivalent shear wave velocity corresponding to each layer of each borehole, the method further comprises: generating equivalent wave velocity data of each borehole according to the equivalent shear wave velocity corresponding to each layer of each borehole; and comparing the equivalent wave velocity data of each borehole with the measured wave velocity data to generate a wave velocity comparison ladder diagram, which is used to verify the accuracy of the equivalent shear wave velocity corresponding to all layers of the borehole. It should be understood that the equivalent wave velocity data of a borehole includes the equivalent shear wave velocity corresponding to all layers of the borehole.

[0091] To verify the reasonableness of the equivalent shear wave velocity calculation, the site seismic motion calculation model construction system provides a comparison function of the measured shear wave velocity and the equivalent shear wave velocity. Figure 13 The schematic diagram of the wave velocity comparison ladder diagram provided by the embodiment of the present application. As shown in Figure 13 The thick line represents the equivalent shear wave velocity, and the thin line represents the measured shear wave velocity.

[0092] By executing step 143, the embodiment of the present application can visually compare and display on the basis of automatically calculating the equivalent shear wave velocity. The site seismic motion calculation model construction system provides a method for automatically calculating the equivalent shear wave velocity, which can automatically and quickly calculate the equivalent shear wave velocity corresponding to each layer of each borehole according to the measured shear wave velocity, and can compare and display the calculated equivalent shear wave velocity and the measured shear wave velocity on a schematic diagram.

[0093] In some possible embodiments, after generating borehole seismic motion calculation model data for each borehole, the method further includes:

[0094] Step 16: Based on the number of boreholes and their numbers, create a unique borehole folder for each borehole. The borehole seismic calculation model data for different boreholes are stored in different borehole folders, and the borehole seismic calculation model data for the same borehole are stored in the same file within the same borehole folder.

[0095] In other words, the system for constructing the site's seismic motion calculation model generates borehole folders based on the number and number of boreholes to store the borehole seismic motion calculation model data. To facilitate borehole seismic motion calculations, the system can generate borehole seismic motion calculation model files from the seismic motion calculation model spreadsheet file and store them in each borehole folder. The file format is consistent with the format required by existing seismic response analysis programs. Each borehole seismic motion calculation file includes soil type, thickness, wave velocity, density, etc.

[0096] Optionally, embodiments of this application may provide step 17 after step 16:

[0097] Step 17: Create an index file in each borehole folder to record the name and location of all files in the borehole folder.

[0098] By executing steps 16 and 17, this embodiment of the application generates borehole folders with a single click. Based on the borehole information, a folder is created for each borehole to store the site ground motion calculation model file for each borehole, avoiding confusion caused by storing everything in one folder. Furthermore, this embodiment of the application can generate site soil layer ground motion calculation input files with a single click. Based on the layering information of each borehole, the layer thickness, soil type, shear wave velocity, density, and other information for each borehole are written into the site ground motion calculation model file of the corresponding borehole folder.

[0099] This application provides a file path setting function, allowing users to set the software's working path according to their own folder paths. After setting the path, reading input files and storing files will both follow this path. Furthermore, this application provides real-time error alerts. Considering that incorrect manual input parameters may lead to task failure, a parameter input error alert function is included. If an error occurs, a pop-up alert will appear.

[0100] It should be understood that after the user clicks the "Close Software" button, the system automatically checks if any graphical windows (or simply window windows) are closed. If so, it closes the window windows to avoid consuming computer resources and screen space. Then, it clears the software data and closes the software. For example, Figure 14 This is a flowchart illustrating the software shutdown step provided in an embodiment of this application.Figure 14 As shown, the process of closing the corresponding software link includes: S141, clicking the button of closing the software. S142, judging whether there are windows not closed. If yes, executing S143, otherwise executing S144. S143, closing the windows. S144, clearing the software data.

[0101] Figure 15 A structural schematic diagram of a construction system of a site seismic ground motion calculation model provided by an embodiment of the present application is shown in the figure, and the system includes: Figure 15

[0102] The display module 151 is configured to display an operation interface in response to a user starting operation of target application software, and the operation interface includes parameter items and operation buttons required for construction of the site seismic ground motion calculation model.

[0103] The configuration reading module 152 is configured to configure the parameter items, and the parameter items include a path of an input file, a maximum layering thickness, a number of drill holes, a hole number prefix, a layering information separator, a soil class name and description separator, a soil class code and name separator, and a measured wave velocity data range. After the parameter configuration is completed, the configuration reading module 152 is configured to read a soil layer information input file in the path of the input file in response to a touch operation of the operation button, and the soil layer information input file includes a drill hole original layering input file, a soil layer category and description input file, and a measured wave velocity input file.

[0104] The identification and division module 153 is configured to identify layering information of each drill hole from the drill hole original layering input file and the soil layer category and description input file according to the number of the drill hole, and if the layering information of the drill hole does not meet a pre-designed calculation condition, the identification and division module 153 is configured to re-divide the layering of each drill hole according to the maximum layering thickness.

[0105] The calculation module 154 is configured to calculate an equivalent shear wave velocity corresponding to each layering of each drill hole by using a numerical equivalent shear wave velocity calculation algorithm according to measured wave velocity data in the measured wave velocity input file.

[0106] The generation construction module 155 is configured to generate drill hole seismic ground motion calculation model data corresponding to each drill hole based on the layering information and the equivalent shear wave velocity corresponding to each layering of each drill hole, to construct a site seismic ground motion calculation model based on the drill hole seismic ground motion calculation model data corresponding to all drill holes and in combination with soil dynamics test information, and to generate a site seismic ground motion calculation model file.

[0107] Figure 15 The construction system of the site seismic ground motion calculation model can perform Figure 1 ​The implementation principle and technical effects of the construction method of the site seismic ground motion calculation model in the illustrated embodiment will not be described again. The specific manner in which each module and unit in the construction system of the site seismic ground motion calculation model in the above-described embodiment performs operations has been described in detail in the embodiments related to the method, and will not be described in detail here.

[0108] In one possible design, Figure 15 The construction system of the site seismic ground motion calculation model in the illustrated embodiment can be implemented as a computing device, such as Figure 16 As illustrated, the computing device can include a storage component 161 and a processing component 162.

[0109] The storage component 161 stores one or more computer instructions, wherein the one or more computer instructions are called for execution by the processing component 162.

[0110] The processing component 162 is configured to: in response to a user starting operation on a target application software, display an operation interface, the operation interface including parameter items and operation buttons required for construction of a site seismic ground motion calculation model; configure the parameter items, the parameter items including a path of an input file, a maximum layering thickness, a number of drill holes, a hole number prefix, a layering information separator, a soil class name and description separator, a soil class code and name separator, and a range of measured wave velocity data, after completing the parameter configuration, in response to a touch operation of the operation button, read a soil layer information input file in the path of the input file, the soil layer information input file including: a drill hole original layering input file, a soil layer category and description input file, and a measured wave velocity input file; according to the number of the drill holes, identify the layering information of each drill hole from the drill hole original layering input file and the soil layer category and description input file, if the layering information of the drill hole does not meet a pre-designed calculation condition, redivide the layering of each drill hole according to the maximum layering thickness; according to the measured wave velocity data in the measured wave velocity input file, calculate the equivalent shear wave velocity corresponding to each layering of each drill hole by using a numerical equivalent shear wave velocity calculation algorithm; based on the layering information of each drill hole and the equivalent shear wave velocity corresponding to the layering, generate drill hole seismic ground motion calculation model data corresponding to each drill hole, and based on the drill hole seismic ground motion calculation model data corresponding to all the drill holes, construct a site seismic ground motion calculation model in combination with soil dynamics test information to generate a site seismic ground motion calculation model file.

[0111] The processing component 162 can include one or more processors to execute computer instructions to complete all or part of the steps in the above method. Of course, the processing component can also be one or more Application-Specific Integrated Circuits (ASICs), Digital Signal Processors (DSPs), Digital Signal Process Device (DSPD), Programmable Logic Devices (PLDs), Field Programmable Gate Arrays (FPGAs), controllers, microcontrollers, microprocessors or other electronic components, for executing the above method.

[0112] The storage component 161 is configured to store various types of data to support the operation of the terminal. The storage component can be implemented by any type of volatile or non-volatile storage devices or their combinations, such as Random Access Memory (RAM), Static Random-Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic storage, flash memory, magnetic disk or optical disk.

[0113] Of course, the computing device can also include other components, such as an input / output interface, a display component, a communication component, etc.

[0114] The input / output interface provides an interface between the processing component and peripheral interface modules, which can be output devices, input devices, etc.

[0115] The communication component is configured to facilitate wired or wireless communication between the computing device and other devices, etc.

[0116] The computing device can be a physical device or an elastic computing host provided by a cloud computing platform, and the computing device can refer to a cloud server, and the processing component, the storage component, etc. can be basic server resources rented or purchased from the cloud computing platform.

[0117] The embodiment of the application further provides a computer storage medium, which stores a computer program. Figure 1 The application further provides a construction method of the site ground motion calculation model of the embodiment.

[0118] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the system and the units described above can refer to the corresponding processes in the foregoing method embodiments, which will not be described here.

[0119] The system embodiments described above are merely illustrative, wherein the units described as separate components can or can not be physically separated, and the components shown as units can or can not be physical units, that is, can be located in one place or can be distributed on multiple network units. Part or all of the modules can be selected to achieve the purpose of the embodiment scheme according to actual needs. Those skilled in the art can understand and implement without creative labor.

[0120] Through the foregoing description of the embodiments, those skilled in the art can clearly understand that each embodiment can be realized by means of software and a necessary general hardware platform, and of course, can also be realized by hardware. Based on such understanding, the foregoing technical solutions can be embodied in the form of a software product, which can be stored in a computer readable storage medium, such as a ROM / RAM, a magnetic disk, an optical disk, etc., and includes a plurality of instructions to make a computer device (which can be a personal computer, a server, or a network device, etc.) execute the method described in each embodiment or some parts of the embodiment.

[0121] Finally, it should be noted that: the foregoing embodiments are only used to illustrate the technical solutions of the application, rather than limit them; although the application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: it can still modify the technical solutions recorded in the foregoing embodiments, or make equivalent replacement to some technical features; and these modifications 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 application.

Claims

1. A method for constructing a site seismic motion calculation model, characterized in that, include: In response to the user's launch operation of the target application software, an operation interface is displayed, which includes the parameter items and operation buttons required for building the site ground motion calculation model; Configure the parameter items, which include the path of the input file, maximum layer thickness, number of boreholes, borehole number prefix, layer information separator, soil type name and description separator, soil type code and name separator, and measured wave velocity data range. After the parameter configuration is completed, in response to the touch operation of the operation button, read the soil layer information input file under the path of the input file. The soil layer information input file includes: the original borehole layer input file, the soil layer type and description input file, and the measured wave velocity input file. Based on the borehole number, the layer information of each borehole is identified from the original layer input file and the soil layer category and description input file. If the layer information of the borehole does not meet the preset calculation conditions, the layer of each borehole is re-divided according to the maximum layer thickness. Based on the measured wave velocity data in the measured wave velocity input file, the equivalent shear wave velocity corresponding to each layer of each borehole is calculated using a numerical equivalent shear wave velocity calculation algorithm. Based on the layer information of each borehole and the equivalent shear wave velocity corresponding to each layer, borehole ground motion calculation model data corresponding to each borehole is generated. Based on the borehole ground motion calculation model data corresponding to all boreholes, combined with soil dynamics test information, a site ground motion calculation model is constructed, and a site ground motion calculation model file is generated. The soil layer information input file also includes: a cover layer thickness input file, which includes the cover layer thickness of each borehole, and the cover layer thickness is the cumulative value of all target layer thicknesses; the measured wave velocity input file includes the measured wave velocity data of all boreholes, and the measured wave velocity data of each borehole includes the measured shear wave velocity of all equally spaced layers, where the thickness of the equally spaced layers is much larger than the fixed thickness used in the numerical equivalent shear wave velocity calculation; "much larger" means that the order of magnitude of the thickness of the equally spaced layers is different from the order of magnitude of the fixed thickness, where the order of magnitude of the thickness of the equally spaced layers is on the order of meters, and the order of magnitude of the fixed thickness is on the order of micrometers; The step involves calculating the equivalent shear wave velocity for each layer of each borehole using the measured wave velocity data from the measured wave velocity input file and employing a numerical equivalent shear wave velocity calculation algorithm. This includes: Each layer of each borehole is designated as the target layer; For each target layer, based on the measured shear wave velocities of equally spaced layers, a numerical equivalent shear wave velocity calculation algorithm is used to discretize the target layer into unit layers of fixed thickness. Based on the shear wave velocity corresponding to each unit layer, the travel time of each unit layer is calculated. The travel times of all unit layers belonging to the same target layer are accumulated to obtain the target layer travel time. The target layer thickness is divided by the target layer travel time to obtain the equivalent shear wave velocity corresponding to each layer of each borehole.

2. The method according to claim 1, characterized in that, After obtaining the equivalent shear wave velocity corresponding to each layer of each borehole, the method further includes: Based on the equivalent shear wave velocity corresponding to each layer of each borehole, the equivalent wave velocity data of each borehole is generated. The equivalent wave velocity data and the measured wave velocity data of each borehole are compared to generate a wave velocity comparison ladder diagram. The wave velocity comparison ladder diagram is used to verify the accuracy of the equivalent shear wave velocity corresponding to all layers of the borehole.

3. The method according to claim 1, characterized in that, The process of re-dividing each borehole into layers based on the maximum layer thickness includes: Select other factors, which include at least one of the following: soil physical properties, soil dynamic properties, and soil type; The layers of each borehole are redefined based on the maximum layer thickness and the other factors mentioned above.

4. The method according to claim 1, characterized in that, After generating borehole seismic motion calculation model data for each borehole, the method further includes: Based on the number of boreholes and their numbers, a unique borehole folder is created for each borehole. The borehole seismic motion calculation model data for different boreholes are stored in different borehole folders, and the borehole seismic motion calculation model data for the same borehole are stored in the same file within the same borehole folder.

5. The method according to claim 1, characterized in that, In configuring the parameter items, the method further includes: The system determines whether the input content of the parameter item meets the corresponding configuration conditions. If it does not meet the conditions, a prompt box is displayed in a designated area of ​​the operation interface to prompt the user to check and modify the input content, thus preventing task failure due to misoperation.

6. The method according to claim 1, characterized in that, The number of operation buttons is multiple, and each operation button corresponds to a step in the process of building a site ground motion calculation model. The operation interface includes a progress status display indicator for each operation button.

7. A system for constructing a site seismic motion calculation model, characterized in that, include: The display module is used to display the operation interface in response to the user's launch operation of the target application software. The operation interface includes the parameter items and operation buttons required for building the site ground motion calculation model. The configuration reading module is used to configure the parameter items, which include the path of the input file, maximum layer thickness, number of boreholes, borehole number prefix, layer information separator, soil type name and description separator, soil type code and name separator, and the range of measured wave velocity data. After the parameter configuration is completed, in response to the touch operation of the operation button, the module reads the soil layer information input file under the path of the input file. The soil layer information input file includes: the original borehole layer input file, the soil layer type and description input file, and the measured wave velocity input file. The identification and division module is used to identify the layer information of each borehole from the original layer input file and the soil layer category and description input file according to the borehole number. If the layer information of the borehole does not meet the preset calculation conditions, the layer of each borehole is re-divided according to the maximum layer thickness. The calculation module is used to calculate the equivalent shear wave velocity corresponding to each layer of each borehole based on the measured wave velocity data in the measured wave velocity input file and using a numerical equivalent shear wave velocity calculation algorithm. The generation module is used to generate borehole ground motion calculation model data for each borehole based on the layer information and the equivalent shear wave velocity corresponding to each layer. Based on the borehole ground motion calculation model data for all boreholes, and combined with soil dynamics test information, the site ground motion calculation model is constructed, and the site ground motion calculation model file is generated. The soil layer information input file also includes: a cover layer thickness input file, which includes the cover layer thickness of each borehole, and the cover layer thickness is the cumulative value of all target layer thicknesses; the measured wave velocity input file includes the measured wave velocity data of all boreholes, and the measured wave velocity data of each borehole includes the measured shear wave velocity of all equally spaced layers, where the thickness of the equally spaced layers is much larger than the fixed thickness used in the numerical equivalent shear wave velocity calculation; "much larger" means that the order of magnitude of the thickness of the equally spaced layers is different from the order of magnitude of the fixed thickness, where the order of magnitude of the thickness of the equally spaced layers is on the order of meters, and the order of magnitude of the fixed thickness is on the order of micrometers; The step involves calculating the equivalent shear wave velocity for each layer of each borehole using the measured wave velocity data from the measured wave velocity input file and employing a numerical equivalent shear wave velocity calculation algorithm. This includes: Each layer of each borehole is designated as the target layer; For each target layer, based on the measured shear wave velocities of equally spaced layers, a numerical equivalent shear wave velocity calculation algorithm is used to discretize the target layer into unit layers of fixed thickness. Based on the shear wave velocity corresponding to each unit layer, the travel time of each unit layer is calculated. The travel times of all unit layers belonging to the same target layer are accumulated to obtain the target layer travel time. The target layer thickness is divided by the target layer travel time to obtain the equivalent shear wave velocity corresponding to each layer of each borehole.

8. A computing device, characterized in that, It includes a processing component and a storage component; the storage component stores one or more computer instructions; the one or more computer instructions are invoked and executed by the processing component to implement a method for constructing a site seismic motion calculation model as described in any one of claims 1 to 6.

9. A computer storage medium, characterized in that, The system contains a computer program that, when executed by a computer, implements a method for constructing a site ground motion calculation model as described in any one of claims 1 to 6.

Images