Method and system for constructing field seismic oscillation calculation model

Abstract

The invention provides a construction method and system of a field seismic oscillation calculation model, and relates to the field of engineering earthquakes. The method comprises the following steps: identifying layering information of each drill hole according to the serial number of the drill hole, and if the layering information of the drill hole does not meet a preset calculation condition, re-dividing the layering of each drill hole according to the maximum layering thickness; according to actually measured wave velocity data in the actually measured wave velocity input file, calculating an equivalent shear wave velocity corresponding to each layer of each drill hole by adopting a numerical equivalent shear wave velocity calculation algorithm; and generating borehole seismic oscillation calculation model data corresponding to each borehole based on the layering information of each borehole and the equivalent shear wave velocity corresponding to the layering, constructing a field seismic oscillation calculation model based on the borehole seismic oscillation calculation model data corresponding to all the boreholes in combination with soil dynamics test information, and generating a field seismic oscillation calculation model file. According to the technical scheme provided by the invention, the construction efficiency and accuracy of the field seismic oscillation calculation model are improved.

Description

Technical Field

[0001] The embodiments of the present application relate to the field of engineering earthquake technology, and in particular, to a method and system for constructing a site earthquake motion calculation model. Background Art

[0002] Domestic and foreign research and earthquake damage examples show that the soil layer covering the site has a strong amplification effect on seismic motion. In the field of engineering earthquakes, the calculation and analysis of site seismic motion is an important part of scientifically determining the seismic fortification parameters of construction projects. The site seismic motion calculation model is mainly constructed based on the type of site covering soil layer, wave velocity structure, soil dynamics properties, etc., among which the wave velocity structure can be described by the equivalent shear wave velocity.

[0003] At present, when constructing a site seismic motion calculation model, technicians generally manually sort out and divide the layers of each borehole according to the histogram of the site's covering soil layer, and then calculate the equivalent shear wave velocity corresponding to each layer. However, this construction method is inefficient, especially when the number of boreholes exceeds 10. The calculation of the equivalent shear wave velocity corresponding to the layers of the boreholes is labor-intensive, time-consuming and prone to errors, which ultimately affects the accuracy of the construction of the site seismic motion calculation model. Summary of the invention

[0004] The embodiments of the present application provide a method and system for constructing a site seismic motion calculation model, so as to solve the problems of low efficiency and poor accuracy in constructing a site seismic motion calculation model in the prior art.

[0005] In a first aspect, an embodiment of the present application provides a method for constructing a site earthquake calculation model, comprising: In response to the user's start-up operation of the target application software, an operation interface is displayed, wherein the operation interface includes parameter items and operation buttons required for constructing the site seismic calculation model; The parameter items are configured, and the parameter items include the path where the input file is located, the maximum layer thickness, the number of boreholes, the hole number prefix, the layer information separator, the soil type name and description separator, the soil type code and name separator, and the measured wave velocity data range. After completing the parameter configuration, in response to the touch operation of the operation button, the soil layer information input file under the path where the input file is located is read, and the soil layer information input file includes: the original layer input file of the borehole, the soil layer type and description input file, and the measured wave velocity input file; According to the borehole number, the stratification information of each borehole is identified from the original stratification input file of the borehole and the soil layer category and description input file. If the stratification information of the borehole does not meet the preset calculation conditions, the stratification of each borehole is re-divided according to the maximum stratification thickness; According to 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 the numerical equivalent shear wave velocity calculation algorithm; Based on the stratification information of each borehole and the equivalent shear wave velocity corresponding to the stratification, the borehole seismic motion calculation model data corresponding to each borehole is generated. Based on the borehole seismic motion calculation model data corresponding to all boreholes and combined with the soil dynamics test information, a site seismic motion calculation model is constructed to generate a site seismic motion calculation model file.

[0006] Optionally, the soil layer information input file further includes: the soil layer information input file further includes: an overburden thickness input file, the overburden thickness input file includes the overburden thickness of each borehole, and the overburden 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, and the thickness of the equally spaced layers is much larger than the fixed thickness used in the numerical equivalent shear wave velocity calculation; The method of calculating 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 includes: Each layer of each drill hole is taken as the target layer; For each of the target layers, according to the measured shear wave velocity of equally spaced layers, a numerical equivalent shear wave velocity calculation algorithm is adopted to discretize the target layer into unit layers of fixed thickness. According to the shear wave velocity corresponding to each unit layer, the travel time of each unit layer is calculated, and the travel time of all unit layers belonging to the same target layer is 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.

[0007] Optionally, after obtaining the equivalent shear wave velocity corresponding to each layer of each borehole, the method further includes: Generate equivalent wave velocity data for each borehole according to the equivalent shear wave velocity corresponding to each layer of each borehole; The equivalent wave velocity data of each borehole is compared with the measured wave velocity data to generate a wave velocity comparison step diagram, which is used to verify the accuracy of the equivalent shear wave velocity corresponding to all layers of the borehole.

[0008] Optionally, the re-dividing the layers of each borehole according to the maximum layer thickness includes: Selecting other factors, wherein the other factors include at least one of the following: soil physical properties, soil dynamic properties and soil type; Each borehole is re-stratified based on the maximum stratum thickness and the other factors described.

[0009] Optionally, after generating the borehole seismic calculation model data corresponding to each borehole, the method further includes: According to the number of boreholes and the borehole numbers, a unique borehole folder is established for each borehole, wherein the borehole seismic calculation model data corresponding to different boreholes are stored in different borehole folders, and the borehole seismic calculation model data of the same borehole are stored in the same file in the same borehole folder.

[0010] Optionally, in the process of configuring the parameter item, the method further includes: It is determined whether the input content of the parameter item meets the corresponding configuration conditions. If not, a prompt box is displayed in a designated area of ​​the operation interface to prompt the user to check and modify the input content to prevent task failure due to misoperation.

[0011] Optionally, there are multiple operation buttons, each operation button corresponds to a link in the process of building a site seismic calculation model, and the operation interface includes a progress status display mark corresponding to each operation button.

[0012] In a second aspect, an embodiment of the present application provides a system for constructing a site earthquake calculation model, including: A display module, for displaying an operation interface in response to a user's startup operation on the target application software, wherein the operation interface includes parameter items and operation buttons required for constructing a site earthquake calculation model; A configuration reading module is used to configure the parameter items, wherein the parameter items include the path where the input file is located, the maximum layer thickness, the number of boreholes, the hole number prefix, the layer information separator, the soil type name and description separator, the soil type code and name separator, and the measured wave velocity data range. After completing the parameter configuration, in response to the touch operation of the operation button, the soil layer information input file under the path where the input file is located is read, wherein the soil layer information input file includes: the original layer input file of the borehole, the soil layer type and description input file, and the measured wave velocity input file; An identification and division module is used to identify the layer information of each borehole from the original layer input file of the borehole 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; A calculation module is used to calculate the equivalent shear wave velocity corresponding to each layer of each borehole according to the measured wave velocity data in the measured wave velocity input file and using a numerical equivalent shear wave velocity calculation algorithm; A construction module is generated, which is used to generate borehole seismic motion calculation model data corresponding to each borehole based on the stratification information of each borehole and the equivalent shear wave velocity corresponding to the stratification, and to construct a site seismic motion calculation model based on the borehole seismic motion calculation model data corresponding to all boreholes in combination with soil dynamics test information, and to generate a site seismic motion calculation model file.

[0013] In a third aspect, an embodiment of the present application provides a computing device, comprising 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 to implement a method for constructing a site seismic calculation model as described in any one of the first aspects.

[0014] In a fourth aspect, an embodiment of the present application provides a computer storage medium storing a computer program, which, when executed by a computer, implements a method for constructing a site seismic calculation model as described in any one of the first aspects.

[0015] In an embodiment of the present application, a method for constructing a site seismic motion calculation model is provided, comprising: in response to a user's startup operation on a target application software, displaying an operation interface, the operation interface including parameter items and operation buttons required for constructing the site seismic motion calculation model; configuring the parameter items, the parameter items including the path where the input file is located, the maximum layer thickness, the number of boreholes, the hole number prefix, the layer information separator, the soil type name and description separator, the soil type code and name separator, and the measured wave velocity data range; after completing the parameter configuration, in response to the touch operation of the operation button, reading the soil layer information input file under the path where the input file is located, the soil layer information input file including: the original layer input file of the borehole, the soil layer category and description input file, and the measured wave velocity input file; according to the borehole number, identify the stratification information of each borehole from the original borehole stratification input file and the soil layer category and description input file; if the stratification information of the borehole does not meet the preset calculation conditions, re-divide the stratification of each borehole according to the maximum stratification thickness; according to the measured wave velocity data in the measured wave velocity input file, adopt the numerical equivalent shear wave velocity calculation algorithm to calculate the equivalent shear wave velocity corresponding to each stratification of each borehole; based on the stratification information of each borehole and the equivalent shear wave velocity corresponding to the stratification, generate the borehole seismic motion calculation model data corresponding to each borehole, based on the borehole seismic motion calculation model data corresponding to all boreholes, combine the soil dynamics test information to construct the site seismic motion calculation model, and generate the site seismic motion calculation model file.

[0016] The embodiment of the present application automatically generates an operation interface by responding to the user's startup operation, and provides a parameter configuration function, reducing the need for manual intervention. The process of reading and parsing the soil layer information input file is automated, greatly shortening the time from data preparation to model generation. Users only need to click the corresponding operation button to complete complex data processing tasks, such as reading files, modeling, etc., simplifying the operation process and improving work efficiency. The soil layer is re-divided based on the maximum layer thickness to ensure that each layer meets the engineering requirements and can accurately reflect the geological characteristics, avoiding simulation errors caused by unreasonable layering. Using a numerical equivalent shear wave velocity calculation algorithm, the target layer is discretized into unit layers of fixed thickness, and the total travel time is calculated by accumulating the travel time and finally obtaining the equivalent shear wave velocity. This method is not only easy to program, 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 calculation model construction method, which can improve the work efficiency and result reliability of site seismic analysis, and provide technical support for the seismic design of buildings.

[0017] These and other aspects of the present application will become more clearly understood in the description of the following embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, a brief introduction will be given below to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0019] Figure 1 A flowchart of a method for constructing a site earthquake calculation model provided in an embodiment of the present application; Figure 2 A schematic diagram of an operation interface provided in an embodiment of the present application; Figure 3 A schematic diagram of the data format of the original layered input file of the drilling provided in the embodiment of the present application; Figure 4 A schematic diagram of the soil layer classification and data format of the description input file provided in the embodiment of the present application; Figure 5 A schematic diagram of the data format of the measured wave velocity input file provided in the embodiment of the present application; Figure 6 A schematic diagram of the data format of a cover thickness input file provided in an embodiment of the present application; Figure 7 A schematic diagram of the reading process of the soil layer information input file provided in the embodiment of the present application; Figure 8 A schematic diagram of the automated construction process of the drilling seismic motion calculation model provided in an embodiment of the present application; Fig. 9 A schematic diagram of error prompts corresponding to the modeling stage provided in an embodiment of the present application; Fig.10 Schematic diagram of the scenario for solving the equivalent shear wave velocity; Fig.11 A schematic diagram of a process for solving the equivalent shear wave velocity provided in an embodiment of the present application; Fig.12 A schematic diagram of an error prompt corresponding to the equivalent shear wave velocity calculation stage provided in an embodiment of the present application; Fig.13 A schematic diagram of a wave velocity comparison ladder diagram provided in an embodiment of the present application; Fig.14 A schematic diagram of the process flow corresponding to the closing software link provided in the embodiment of the present application; Fig.15 A schematic diagram of the structure of a system for constructing a site earthquake calculation model provided in an embodiment of the present application; Fig.16 A schematic diagram of the structure of a computing device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0020] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application.

[0021] In some of the processes described in the specification and claims of this application and the above-mentioned figures, multiple operations that appear in a specific order are included, but it should be clearly understood that these operations may not be executed in the order in which they appear in this article or executed in parallel. The sequence numbers of the operations, such as 11, 12, etc., are only used to distinguish between different operations, and the sequence numbers themselves do not represent any execution order. In addition, these processes may include more or fewer operations, and these operations may be executed in sequence or in parallel. It should be noted that the descriptions of "first", "second", etc. in this article are used to distinguish different messages, devices, modules, etc., do not represent the order of precedence, and do not limit the "first" and "second" to different types.

[0022] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of this application.

[0023] In order to improve the efficiency and quality of building a site seismic calculation model, the present application embodiment proposes a method for building a site seismic calculation model. For example, Figure 1 A flowchart of a method for constructing a site earthquake calculation model provided in an embodiment of the present application is shown in FIG. Figure 1 As shown, the method includes: S11 . In response to the user's start-up operation on the target application software, an operation interface is displayed, where the operation interface includes parameter items and operation buttons required for constructing a site seismic calculation model.

[0024] It should be understood that the user may be referred to as an operator, the operation interface may be referred to as a drawing operation interface, and the system for constructing a site earthquake calculation model may be referred to as an automatic creation system for a site earthquake calculation model.

[0025] Figure 2 Schematic diagram of the operation interface provided by the embodiment of the present application. The user opens the site seismic calculation model automatic creation system and automatically generates an operation interface. The operation interface displays all the parameters and operation buttons required for the site seismic calculation model automatic creation. Figure 2 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 (meter), the layer information separator, the soil type name and description separator, the soil type code and name separator, the test wave velocity data range (i.e. the measured wave velocity data range), the total number of modeled soil types, the form where the soil layer model is located, the soil layer model data range, etc. The operation buttons include reading the soil layer information input file, modeling, creating a borehole folder, creating a model calculation file, closing the software, etc. The above maximum layer interval is also called the maximum layer thickness.

[0026] Optionally, there are multiple operation buttons, each of which corresponds to a link in the process of building a site earthquake calculation model, and the operation interface includes a progress status display mark corresponding to each operation button. Figure 2As shown, the progress status display mark can be a work stage prompt control, which defaults to a color (such as red, gray, etc.) at the beginning, used to indicate an unprocessed state, and different states are indicated by color changes. For example, when the work stage prompt control is yellow, it is used to indicate that it is in process, and when the work stage prompt control is green, it is used to indicate that the process is completed or the process is finished. It can be seen that the embodiment of the present application has the advantage of clear operation flow and can display the stage in real time. The construction system of the site seismic calculation model divides the modeling into working stages such as "reading soil layer information input file", "modeling", "creating drilling folder", "creating model calculation file", and "closing software" according to the practical experience of soil layer model modeling. When using the construction system of the site seismic calculation model, click these buttons in sequence to complete the modeling. Each stage has a completion prompt, for example, red indicates that the stage is not completed, and green indicates that the stage is completed, so as to prevent the user from forgetting the working stage due to being temporarily busy with other things.

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

[0028] In this embodiment, the soil layer information input file also includes: cover layer thickness input file, etc. In this embodiment, the parameter item is also called parameter. Optionally, commonly used parameters can be displayed in the parameter input box as default values. Users can change them according to their needs. For example, the embodiment of this application can be Figure 1 The parameters are configured in the parameter input box in the operation interface. Therefore, the embodiment of the present application has the following advantages: strong human-computer interaction and more convenient use. The default parameters are directly displayed in the input box of the operation interface, the user can clearly obtain them, and can be adaptively modified according to needs.

[0029] It should be understood that the format of the soil layer information input file is consistent with the format of the histogram information file in the existing survey software. For example, the data format of the original layer input file of the borehole is as follows: Figure 3 As shown in the figure, each row is used to store soil type code, layer bottom depth and other information. Among them, for the first row of each borehole, the corresponding borehole number (or hole number, serial number, etc.) is additionally stored. The data format of the soil layer category and description input file is as follows Figure 4As shown in the figure, each row is used to store information such as soil type code, soil type name, soil type description, etc. Optionally, to enhance organization, the soil layer velocity can be stored in a spreadsheet, and the test velocity data range can be filled in the operation interface for use. The above measured velocity input file is also called the test velocity input file. The data format of the measured velocity input file is as follows: Figure 5 As shown, 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. The measured shear wave velocity may also be called measured wave velocity data, test wave velocity, test shear wave velocity, test wave velocity data, etc. The thickness of the equally spaced layers is 1m. Specifically, the measured wave velocity data (or shear wave velocity value) of a borehole is arranged as a column from top to bottom with increasing depth (generally the depth increment step is 1m), and the measured wave velocity data of all boreholes are arranged from left to right. The overburden thickness input file can be a data (Data, DAT) file. DAT file is a general file format, which is usually used to store various types of data. Its extension ".dat" indicates that the file contains raw data. For example, the overburden thickness information of all boreholes is placed in the data file and arranged from left to right, such as Figure 6 As shown, the thickness of the covering layer in the four boreholes is 13, 19, 12 and 19 respectively.

[0030] For example, Figure 7 This is a schematic diagram of the reading process of the soil layer information input file provided in the embodiment of the present application. Figure 7 As shown, the reading process includes: S71, click the button of "read soil layer information input file". S72, read the original stratification input file of the borehole, the soil layer category and description input file, the measured wave velocity input file, and the cover layer thickness input file. S73, determine whether the hole numbers are arranged in sequence, if so, execute S74, otherwise execute S76. S74, determine whether the buried depth of the layers is arranged in sequence, if so, execute S75, otherwise execute S76. S75, form an array of borehole soil properties, layer bottom depth, measured wave velocity, and position. S76, end "read soil layer information input file".

[0031] Optionally, during the process of configuring parameter items, the method also includes: determining whether the input content of the parameter item meets the corresponding configuration conditions; if not, displaying a prompt box in a designated area of ​​the operation interface to prompt the user to check and modify the input content to prevent task failure due to misoperation.

[0032] S13, according to the borehole number, identify the layer information of each borehole from the original layer input file of the borehole and the soil layer category and description input file, if the layer information of the borehole does not meet the preset calculation conditions, then re-divide the layer of each borehole according to the maximum layer thickness. The borehole numbering starts from 1 and increases by default.

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

[0034] S15. Generate borehole seismic motion calculation model data corresponding to each borehole based on the stratification information of each borehole and the equivalent shear wave velocity corresponding to the stratification. Construct a site seismic motion calculation model based on the borehole seismic motion calculation model data corresponding to all boreholes and in combination with soil dynamics test information, and generate a site seismic motion calculation model file.

[0035] From the description of S13 to S15, it can be known that the main idea of ​​automatically constructing the borehole seismic motion calculation model is to re-layer the soil layer of the borehole according to the input soil layer properties, layer bottom depth, test wave velocity and other information according to the set maximum layer thickness, establish a layer suitable for calculation, generate a fine layered borehole model matrix, and then calculate the equivalent shear wave velocity corresponding to each layer of the borehole according to the test wave velocity, generate the borehole seismic motion calculation model data, and output the borehole seismic motion calculation model data file. Exemplarily, the automatic construction process of the borehole seismic motion calculation model is as follows: Figure 8 As shown, it includes: S81, generating an initial drilling model matrix according to input information. S82, re-stratifying according to factors such as maximum layer thickness. S83, generating a fine stratified drilling model matrix. S84, judging whether the layer thickness is uniform, if so, executing S86, otherwise executing S85. S85, adjusting the stratification to regenerate a fine stratified drilling model matrix. S86, calculating the equivalent shear wave velocity corresponding to each layer. S87, generating a stratified drilling model matrix with wave velocity. S88, performing bedrock layer processing. S89, generating a final stratified drilling model matrix. S90, outputting a drilling seismic calculation model data file.

[0036] Afterwards, the user can assign soil type numbers and other data to each type of soil in a spreadsheet based on the processed dynamic triaxial test data (which can be referred to as dynamic triaxial data for short), and form a complete seismic motion calculation model spreadsheet file. Therefore, the embodiment of the present application has the advantage of being in line with the actual working situation. 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 more flexible. If calculated directly by computer, it may not be completely consistent with human expectations, and the time spent on modifying the model may exceed the time required for direct manual processing. Therefore, the selection, classification and merging of dynamic triaxial data are handled by professional and technical personnel. Perform operations in the spreadsheet to assign soil type information to the soil layer.

[0037] To avoid program errors caused by input parameter errors (which can refer to the above parameters), configure a try...catch statement to promptly remind users to check and modify the input parameters, such as Fig. 9 shown.

[0038] By executing steps S11 to S15, it can be known that the present embodiment can obtain the histogram information input file of the covering soil layer, and read the histogram information of the covering soil layer of the site from the histogram information input file with one click. The format of the histogram information input file is consistent with the format of the existing survey software. After completing the survey of the covering soil layer of the site, the present application embodiment can use its histogram information file as the input file of the construction system of the site seismic calculation model, without the need for re-editing, thereby reducing duplication of work. All information can be read in by simply clicking the corresponding button with the mouse. Similarly, the soil layer information input file of the covering soil layer is obtained, and the stratification information of the 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 the foundation for re-dividing the stratification of each borehole, finding the corresponding soil layer wave velocity, and calculating the equivalent shear wave velocity.

[0039] The embodiment of the present application automatically generates an operation interface in response to the user's startup operation and provides a parameter configuration function, thereby reducing the need for manual intervention. The process of reading and parsing the soil layer information input file is automated, greatly shortening the time from data preparation to model generation. Users only need to click the corresponding operation button to complete complex data processing tasks, such as reading files, modeling, etc., which simplifies the operation process and improves work efficiency. The soil layer is re-divided based on the maximum layer thickness to ensure that each layer meets the engineering requirements and can accurately reflect the geological characteristics, 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 of fixed thickness (for example, 1E-6m), and the total travel time is calculated by accumulating the travel time and finally obtaining the equivalent shear wave velocity. This method is not only easy to program, 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 method for constructing a site seismic motion calculation model, which can improve the work efficiency and result reliability of site seismic motion analysis, and provide technical support for the seismic design of buildings.

[0040] In some possible embodiments, S13, re-dividing the layers of each borehole according to the maximum layer thickness, includes: Step 131: Select other factors, where the other factors include at least one of the following: soil physical properties, soil dynamic properties, and soil type.

[0041] Step 132: re-dividing the layers of each borehole based on the maximum layer thickness and other factors.

[0042] By executing steps 131 and 132, the embodiment of the present application can realize automatic division of soil layers. Figure 1 Generally, the soil is divided into layers according to the soil type at different locations. This traditional layering method generally does not meet the layering requirements of the site seismic calculation model. Therefore, the embodiment of the present application can re-layer according to the soil properties and layer thickness requirements. The construction system of the site seismic calculation model can automatically re-divide the soil layers for each input borehole.

[0043] Fig.10 The schematic diagram of the scene for solving the equivalent shear wave velocity is used to solve the common equivalent shear wave velocity calculation problem. The corresponding traditional equivalent shear wave velocity calculation formula is shown in equations (1) to (3).

[0044] (1) (2) (3) in, , and It is used to express the equivalent shear wave velocity of layer 1, layer 2 and layer 3 of a borehole obtained by the traditional equivalent shear wave velocity calculation formula. The layer thickness of layer 1 is , in m, the travel time of layer 1 is , the layer thickness of layer 2 is , in m, the travel time of layer 2 is , the layer thickness of layer 3 is , in m, the travel time of layer 3 is .in, , , , and is the test wave velocity in equally spaced layer 1, equally spaced layer 2, equally spaced layer 3, equally spaced layer 4 and equally spaced layer 5, is the distance difference between the lower line of layer 1 and the lower line of equally spaced layer 1 of the equivalent shear wave velocity model to be determined, is the distance difference between the lower line of layer 2 and the lower line of equally spaced layer 3 of the equivalent shear wave velocity model to be determined, is the distance difference between the lower line of layer 3 and the lower line of equally spaced layer 4 of the equivalent shear wave velocity model to be determined.

[0045] From equations (1) to (3), the calculation of the equivalent shear wave velocity requires the layer thickness and travel time to be calculated according to different stratification conditions, and then divided. The complexity is mainly reflected in the calculation process of the travel time. For different stratum divisions, the formulas to be solved are very different, so it is more complicated to solve this problem with conventional programming ideas. To simplify this problem, the embodiment of the present application proposes a numerical equivalent shear wave velocity calculation algorithm. The main idea of ​​the algorithm is to discretize all the soil layers of the borehole into a series of thin layers (or unit layers), and calculate the shear wave travel time of each thin layer in advance according to the test wave velocity. In this way, the calculation of the travel time is simplified to the sum of the travel times of all thin layers corresponding to the stratum to be calculated, and the equivalent shear wave velocity can be obtained by directly dividing by the layer thickness. The more thin layers there are, the more accurate the value obtained. Considering the improvement of computer performance, sufficient calculation accuracy can be obtained when each thin layer reaches 1E-6m. Therefore, the soil layer information input file also includes: an overburden thickness input file, the overburden thickness input file includes the overburden thickness of each borehole, the overburden thickness is the cumulative value of the thickness of all target layers, and is used to determine the upper and lower limits of the thickness of all target layers; the measured wave velocity input file includes the measured shear wave velocity of the equally spaced layers, the thickness of the equally spaced layers is much (for example, 1m) greater than the fixed thickness (for example, 1E-6m), and much greater 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, and there is an obvious 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, calculate the equivalent shear wave velocity corresponding to each layer of each borehole, including: Step 141: Take each layer of each borehole as a target layer.

[0046] Step 142: For each target layer, based on the measured shear wave velocity 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. The travel time of each unit layer is calculated based on the shear wave velocity corresponding to each unit layer. The travel time of all unit layers belonging to the same target layer is accumulated to obtain the travel time of the target layer. The target layer thickness is divided by the travel time of the target layer to obtain the equivalent shear wave velocity corresponding to each layer of each borehole.

[0047] Fig.11 The flowchart of solving the equivalent shear wave velocity provided in the embodiment of the present application is as follows. Fig.11As shown, the following steps are included: S111. Read the burial depth array of the layer to be determined. S112. Read the test wave velocity array. S113. Discretize the test wave velocity layer into thin layers. S114. Calculate the discrete thin layer travel time array. S115. Determine whether i≤the number of layers to be determined. If so, execute S116, otherwise execute S121. S116. Calculate the thickness of the i-th layer. S117. Find the thin layer travel time array index corresponding to the upper and lower burial depths of the i-th layer. S118. Calculate the sum of the elements within the travel time array index range to obtain the travel time of the i-th layer. S119. Divide 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, continue to execute S115.

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

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

[0050] In some possible embodiments, after obtaining the equivalent shear wave velocity corresponding to each layer of each borehole, the method further includes: generating equivalent wave velocity data of each borehole according to the equivalent shear wave velocity corresponding to each layer of each borehole; comparing the equivalent wave velocity data of each borehole with the measured wave velocity data 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. 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.

[0051] In order to verify the rationality of the equivalent shear wave velocity calculation, the construction system of the site seismic calculation model provides a function for comparing the measured shear wave velocity with the equivalent shear wave velocity. Fig.13 Schematic diagram of the wave velocity comparison ladder diagram provided in the embodiment of the present application. Fig.13 As shown, the thick line represents the equivalent shear wave velocity, and the thin line represents the measured shear wave velocity.

[0052] By executing step 143, the embodiment of the present application can perform intuitive comparative display based on the automatic calculation of the equivalent shear wave velocity. The construction system of the site seismic calculation model 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 based on the measured shear wave velocity, and can put the calculated equivalent shear wave velocity and the measured shear wave velocity on a schematic diagram for comparative display.

[0053] In some possible embodiments, after generating the borehole seismic calculation model data corresponding to each borehole, the method further includes: Step 16: Create a unique borehole folder for each borehole according to the number of boreholes and the borehole numbers, wherein the borehole seismic calculation model data corresponding to different boreholes are stored in different borehole folders, and the borehole seismic calculation model data of the same borehole are stored in the same file in the same borehole folder.

[0054] In other words, the construction system of the site seismic motion calculation model generates a borehole folder according to the number and serial number of the boreholes, and stores the borehole seismic motion calculation model data. To facilitate the borehole seismic motion calculation, the construction system of the site seismic motion calculation model can generate a borehole seismic motion calculation model file according to the seismic motion calculation model spreadsheet file, and store it in each borehole folder. The file format is consistent with the format required by the existing earthquake response analysis program. Each borehole seismic motion calculation file includes soil type, thickness, wave velocity, density, etc.

[0055] Optionally, the embodiment of the present application may provide step 17 after step 16: Step 17: Create an index file in each drilling folder to record the names and locations of all files in the drilling folder through the index file.

[0056] By executing steps 16 to 17, the embodiment of the present application generates a borehole folder with one click. According to the borehole information, a folder is created for each borehole to store the site seismic motion calculation model file of each borehole to avoid confusion caused by placing them all in one folder. In addition, the embodiment of the present application can generate a site soil layer seismic motion calculation input file with one click. According to the stratification and other information of each borehole, the stratification thickness, soil type, shear wave velocity, density and other information of each borehole are written into the site seismic motion calculation model file of the corresponding borehole folder.

[0057] The embodiment of the present application provides a file path setting function, and the user can set the software working path according to the path of his own folder. After setting the path, reading input files and storing files are all carried out according to this path. In addition, the embodiment of the present application can realize real-time error reminders. Considering that incorrect manual input parameters may cause task failure, a parameter input error reminder function is set. If an error occurs, there will be a pop-up window to remind you.

[0058] It should be understood that after the user clicks the "Close Software" button, the system automatically determines whether there is a graphics window (referred to as a figure window) that has not been closed. If so, the figure window is closed to avoid occupying computer resources and the screen. Then the software data is cleared and the software is closed. For example, Fig.14 This is a flow chart of the corresponding closing software link provided in the embodiment of the present application. Fig.14As shown, the process corresponding to the closing software link includes: S141, click the "Close Software" button. S142, determine whether there is a window that is not closed. If so, execute S143, otherwise execute S144. S143, close the window. S144, clear the software data.

[0059] Fig.15 A schematic diagram of a system for constructing a site earthquake calculation model provided in an embodiment of the present application is shown in FIG. Fig.15 As shown, the system includes: The display module 151 is used to display an operation interface in response to a user's start-up operation on the target application software. The operation interface includes parameter items and operation buttons required for constructing the site seismic calculation model.

[0060] The configuration reading module 152 is used to configure parameter items, which include the path where the input file is located, the maximum layer thickness, the number of drill holes, the hole number prefix, the layer information separator, the soil type name and description separator, the soil type code and name separator, and the measured wave velocity data range. After completing the parameter configuration, in response to the touch operation of the operation button, the soil layer information input file under the path where the input file is located is read. The soil layer information input file includes: the original layer input file of the drill hole, the soil layer category and description input file, and the measured wave velocity input file.

[0061] The identification and division module 153 is used to identify the stratification information of each borehole from the original stratification input file of the borehole and the soil layer category and description input file according to the borehole number. If the stratification information of the borehole does not meet the preset calculation conditions, the stratification of each borehole is re-divided according to the maximum stratification thickness.

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

[0063] Generate a construction module 155, which is used to generate borehole seismic motion calculation model data corresponding to each borehole based on the stratification information of each borehole and the equivalent shear wave velocity corresponding to the stratification, build a site seismic motion calculation model based on the borehole seismic motion calculation model data corresponding to all boreholes and combine the soil dynamics test information to generate a site seismic motion calculation model file.

[0064] Fig.15 The construction system of the site earthquake calculation model can execute Figure 1The implementation principle and technical effect of the method for constructing the site seismic calculation model described in the illustrated embodiment will not be described in detail. The specific manner in which each module and unit performs operations in the construction system of the site seismic calculation model in the above embodiment has been described in detail in the embodiment of the method, and will not be described in detail here.

[0065] In one possible design, Fig.15 The construction system of the site earthquake calculation model of the embodiment shown can be implemented as a computing device, such as Fig.16 As shown, the computing device may include a storage component 151 and a processing component 162 .

[0066] The storage component 161 stores one or more computer instructions, wherein the one or more computer instructions are invoked and executed by the processing component 162 .

[0067] The processing component 162 is used to: in response to the user's start-up operation of the target application software, display an operation interface, the operation interface includes parameter items and operation buttons required for constructing the site seismic calculation model; configure parameter items, the parameter items include the path where the input file is located, the maximum layer thickness, the number of boreholes, the hole number prefix, the layer information separator, the soil type name and description separator, the soil type code and name separator, and the measured wave velocity data range, after completing the parameter configuration, in response to the touch operation of the operation button, read the soil layer information input file under the path where the input file is located, the soil layer information input file includes: the original layer input file of the borehole, the soil layer category and description input file, and the measured wave velocity input file; according to the borehole The stratification information of each borehole is identified from the original stratification input file of the borehole and the soil layer category and description input file. If the stratification information of the borehole does not meet the preset calculation conditions, the stratification of each borehole is re-divided according to the maximum stratification thickness; the equivalent shear wave velocity corresponding to each stratum of each borehole is calculated using the numerical equivalent shear wave velocity calculation algorithm according to the measured wave velocity data in the measured wave velocity input file; based on the stratification information of each borehole and the equivalent shear wave velocity corresponding to the stratification, the borehole seismic motion calculation model data corresponding to each borehole is generated; based on the borehole seismic motion calculation model data corresponding to all boreholes, the site seismic motion calculation model is constructed in combination with the soil dynamics test information, and the site seismic motion calculation model file is generated.

[0068] The processing component 162 may 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 may also be implemented by one or more application-specific integrated circuits (ASIC), digital signal processors (DSP), digital signal processing devices (DSPD), programmable logic devices (PLD), field programmable gate arrays (FPGA), controllers, microcontrollers, microprocessors or other electronic components to perform the above method.

[0069] The storage component 161 is configured to store various types of data to support operations at the terminal. The storage component can be implemented by any type of volatile or non-volatile storage device or a combination thereof, 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 memory, flash memory, magnetic disk or optical disk.

[0070] Of course, the computing device may also include other components, such as input / output interfaces, display components, communication components, etc.

[0071] The input / output interface provides an interface between the processing component and the peripheral interface module, which may be an output device, an input device, etc.

[0072] The communication component is configured to facilitate, among other things, wired or wireless communications between the computing device and other devices.

[0073] Among them, the computing device can be a physical device or an elastic computing host provided by a cloud computing platform, etc. In this case, the computing device can refer to a cloud server, and the above-mentioned processing components, storage components, etc. can be basic server resources rented or purchased from the cloud computing platform.

[0074] The present application also provides a computer storage medium storing a computer program, wherein the computer program can achieve the above-mentioned Figure 1 The method for constructing a site seismic calculation model of the illustrated embodiment.

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

[0076] The system embodiment described above is merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art may understand and implement it without creative work.

[0077] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0078] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit it. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for constructing a site earthquake calculation model, characterized in that: include: In response to the user's start-up operation of the target application software, an operation interface is displayed, wherein the operation interface includes parameter items and operation buttons required for constructing the site seismic calculation model; The parameter items are configured, and the parameter items include the path where the input file is located, the maximum layer thickness, the number of boreholes, the hole number prefix, the layer information separator, the soil type name and description separator, the soil type code and name separator, and the measured wave velocity data range. After completing the parameter configuration, in response to the touch operation of the operation button, the soil layer information input file under the path where the input file is located is read, and the soil layer information input file includes: the original layer input file of the borehole, the soil layer type and description input file, and the measured wave velocity input file; According to the borehole number, the stratification information of each borehole is identified from the original stratification input file of the borehole and the soil layer category and description input file. If the stratification information of the borehole does not meet the preset calculation conditions, the stratification of each borehole is re-divided according to the maximum stratification thickness; According to 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 the numerical equivalent shear wave velocity calculation algorithm; Based on the stratification information of each borehole and the equivalent shear wave velocity corresponding to the stratification, the borehole seismic motion calculation model data corresponding to each borehole is generated. Based on the borehole seismic motion calculation model data corresponding to all boreholes and combined with the soil dynamics test information, a site seismic motion calculation model is constructed to generate a site seismic motion calculation model file.

2. The method according to claim 1, characterized in that The soil layer information input file also includes: an overburden thickness input file, the overburden thickness input file includes the overburden thickness of each borehole, and the overburden thickness is the cumulative value of the thickness of all target layers; 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, and the thickness of the equally spaced layers is much larger than the fixed thickness used in the numerical equivalent shear wave velocity calculation; the 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, the order of magnitude of the thickness of the equally spaced layers is meter level, and the order of magnitude of the fixed thickness is micrometer level; The method of calculating 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 includes: Each layer of each drill hole is taken as the target layer; For each of the target layers, according to the measured shear wave velocity of equally spaced layers, a numerical equivalent shear wave velocity calculation algorithm is adopted to discretize the target layer into unit layers of fixed thickness. According to the shear wave velocity corresponding to each unit layer, the travel time of each unit layer is calculated, and the travel time of all unit layers belonging to the same target layer is 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.

3. The method according to claim 2, characterized in that After obtaining the equivalent shear wave velocity corresponding to each layer of each borehole, the method further includes: Generate equivalent wave velocity data for each borehole according to the equivalent shear wave velocity corresponding to each layer of each borehole; The equivalent wave velocity data of each borehole is compared with the measured wave velocity data to generate a wave velocity comparison step diagram, which is used to verify the accuracy of the equivalent shear wave velocity corresponding to all layers of the borehole.

4. The method according to claim 1, characterized in that The re-dividing the layers of each borehole according to the maximum layer thickness comprises: Selecting other factors, wherein the other factors include at least one of the following: soil physical properties, soil dynamic properties and soil type; Each borehole is re-stratified based on the maximum stratum thickness and the other factors described.

5. The method according to claim 1, characterized in that After generating the borehole seismic calculation model data corresponding to each borehole, the method further includes: According to the number of boreholes and the borehole numbers, a unique borehole folder is established for each borehole, wherein the borehole seismic calculation model data corresponding to different boreholes are stored in different borehole folders, and the borehole seismic calculation model data of the same borehole are stored in the same file in the same borehole folder.

6. The method according to claim 1, characterized in that In the process of configuring the parameter items, the method further includes: It is determined whether the input content of the parameter item meets the corresponding configuration conditions. If not, a prompt box is displayed in a designated area of ​​the operation interface to prompt the user to check and modify the input content to prevent task failure due to misoperation.

7. The method according to claim 1, characterized in that There are multiple operation buttons, each of which corresponds to a link in the process of building a site seismic calculation model, and the operation interface includes a progress status display mark corresponding to each operation button.

8. A system for constructing a site earthquake calculation model, characterized in that: include: A display module, for displaying an operation interface in response to a user's startup operation on the target application software, wherein the operation interface includes parameter items and operation buttons required for constructing a site earthquake calculation model; A configuration reading module is used to configure the parameter items, wherein the parameter items include the path where the input file is located, the maximum layer thickness, the number of boreholes, the hole number prefix, the layer information separator, the soil type name and description separator, the soil type code and name separator, and the measured wave velocity data range. After completing the parameter configuration, in response to the touch operation of the operation button, the soil layer information input file under the path where the input file is located is read, wherein the soil layer information input file includes: the original layer input file of the borehole, the soil layer type and description input file, and the measured wave velocity input file; An identification and division module is used to identify the layer information of each borehole from the original layer input file of the borehole 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; A calculation module is used to calculate the equivalent shear wave velocity corresponding to each layer of each borehole according to the measured wave velocity data in the measured wave velocity input file and using a numerical equivalent shear wave velocity calculation algorithm; A construction module is generated, which is used to generate borehole seismic motion calculation model data corresponding to each borehole based on the stratification information of each borehole and the equivalent shear wave velocity corresponding to the stratification, and to construct a site seismic motion calculation model based on the borehole seismic motion calculation model data corresponding to all boreholes in combination with soil dynamics test information, and to generate a site seismic motion calculation model file.

9. A computing device, characterized in that It comprises 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 to implement a method for constructing a site seismic calculation model as described in any one of claims 1 to 7.

10. A computer storage medium, characterized in that: A computer program is stored, and when the computer program is executed by a computer, a method for constructing a site seismic calculation model as described in any one of claims 1 to 7 is implemented.

Images