A sand liquefaction influence range evaluation method, device, equipment and medium

Abstract

The present application relates to the technical field of data processing, in particular to a sand liquefaction influence range evaluation method, device, equipment and medium. According to the soil physical parameters obtained at multiple drilling positions in the target area, the actual cyclic stress ratio and the critical cyclic stress ratio of each drilling position are calculated; the base safety factor is calculated according to the ratio of the critical cyclic stress ratio to the actual cyclic stress ratio; the base safety factor is corrected according to the coarse sand mass proportion of each drilling position to obtain the corrected safety factor; the liquefaction index of each drilling position is obtained according to the corrected safety factor range of the corrected safety factor of each drilling position and the liquefaction index corresponding to the corrected safety factor range; the liquefaction index spatial distribution map of the target area is obtained according to the liquefaction index of all drilling positions and the preset geological complexity index. The present application provides a clear spatial range basis for the targeted arrangement of anti-liquefaction measures, and reduces the judgment ambiguity in engineering decision-making.

Description

Technical Field

[0001] This invention relates to the field of data processing technology, specifically to a method, apparatus, equipment, and medium for assessing the impact range of sand liquefaction. Background Technology

[0002] Soil liquefaction is the phenomenon of saturated sand losing strength and stiffness under seismic loading, posing a serious threat to the safety of engineering structures. In site investigation and seismic design, it is necessary to assess the potential area of ​​liquefaction in order to implement targeted anti-liquefaction measures. Currently, widely used liquefaction assessment methods in engineering practice, such as the simplified method based on the standard penetration test (SPT), can provide a judgment on whether liquefaction has occurred at a single borehole location. However, these methods are single-point assessments. For the entire site, the boundary lines of the liquefaction influence area need to be manually delineated based on the assessment results of a limited number of boreholes. Their accuracy and reliability are limited by the engineer's experience and judgment, are highly subjective, and fail to reflect the continuous changing characteristics of liquefaction risk within the site. Summary of the Invention

[0003] The purpose of this invention is to provide a method, apparatus, equipment and medium for assessing the impact range of sand liquefaction, which solves the problems in the prior art.

[0004] This invention is achieved through the following technical solution:

[0005] In a first aspect, embodiments of the present invention provide a method for assessing the impact range of sand liquefaction, including:

[0006] Based on the soil physical parameters obtained from multiple borehole locations within the target area, calculate the actual cyclic stress ratio and critical cyclic stress ratio for each borehole location;

[0007] The basic safety factor is calculated based on the ratio of the critical cyclic stress ratio to the actual cyclic stress ratio.

[0008] The basic safety factor is corrected based on the proportion of coarse sand at each borehole location to obtain the corrected safety factor.

[0009] The liquefaction index for each borehole location is obtained based on the range of corrected safety factors in which the corrected safety factor for each borehole location falls and the liquefaction index corresponding to the range of corrected safety factors.

[0010] Based on the liquefaction index of all borehole locations and the preset geological complexity index, a spatial distribution map of the liquefaction index of the target area is obtained.

[0011] Preferably, the calculation of the actual cyclic stress at each borehole location includes:

[0012] For each borehole location, the product value is obtained based on the ratio of peak ground acceleration to gravitational acceleration in the soil physical parameters corresponding to the borehole location and a preset coefficient.

[0013] The stress ratio is obtained by taking the ratio of the total vertical stress to the effective vertical stress at the borehole depth in the soil physical parameters.

[0014] The depth reduction factor is calculated based on the borehole depth using a preset depth reduction factor formula.

[0015] The actual cyclic stress ratio is obtained by multiplying the product value, the stress ratio, and the depth reduction factor.

[0016] Preferably, calculating the critical cyclic stress ratio at each borehole location includes:

[0017] For each borehole location, a first intermediate value is calculated based on the standard penetration blow count in the soil physical parameters corresponding to the borehole location;

[0018] A second intermediate value is calculated based on the standard penetration number, and the second intermediate value increases linearly with the increase of the standard penetration number.

[0019] A third intermediate value is calculated based on the standard penetration number, and the third intermediate value decreases according to a second preset function as the standard penetration number increases.

[0020] The first intermediate value and the second intermediate value are added together, and then the third intermediate value is subtracted to obtain the basic critical cyclic stress ratio.

[0021] The incremental correction value is calculated based on the fine particle content in the soil physical parameters.

[0022] The critical cyclic stress ratio is obtained based on the basic critical cyclic stress ratio and the incremental correction value.

[0023] Preferably, the step of correcting the basic safety factor based on the proportion of coarse sand at each borehole location to obtain a corrected safety factor includes:

[0024] The total sand mass percentage is calculated based on the mass percentages of fine sand, medium sand, and coarse sand.

[0025] The correction increment is determined based on the ratio of the mass percentage of coarse sand to the mass percentage of total sand.

[0026] Based on the aforementioned correction increment, the particle size distribution correction coefficient is obtained;

[0027] The basic safety factor is corrected based on the particle size distribution correction factor to obtain the corrected safety factor.

[0028] Preferably, obtaining the liquefaction index for each borehole location based on the range of corrected safety factors within which the corrected safety factor for each borehole location falls, and the liquefaction index corresponding to the range of corrected safety factors, includes:

[0029] When the corrected safety factor is less than or equal to the first preset threshold, the liquefaction index takes the first preset value;

[0030] When the safety factor is greater than or equal to the second preset threshold, the liquefaction index takes the second preset value, which is less than the first preset value.

[0031] When the corrected safety factor is between the first preset threshold and the second preset threshold, a proportional value is obtained based on the first difference between the corrected safety factor and the first preset threshold, the second difference between the second preset threshold and the first preset threshold, and the ratio of the first difference and the second difference.

[0032] The liquefaction index is obtained by linearly interpolating the ratio value between the first preset value and the second preset value.

[0033] Preferably, the geological complexity index is preset based on the degree of drastic changes in strata, the range of changes in groundwater depth, and the presence or absence of faults or lenses.

[0034] Preferably, obtaining the spatial distribution map of the liquefaction index of the target area based on the liquefaction index of all borehole locations and a preset geological complexity index includes:

[0035] Using the planar coordinates and liquefaction index of all borehole locations as known points, a preset interpolation algorithm is used to perform spatial interpolation to obtain the interpolated liquefaction index of each grid node.

[0036] The final liquefaction index raster value is obtained based on the interpolated liquefaction index of each raster node and the preset geological complexity index.

[0037] Secondly, embodiments of the present invention provide a device for assessing the impact range of sand liquefaction, comprising:

[0038] The first calculation module is used to calculate the actual cyclic stress ratio and critical cyclic stress ratio for each borehole location based on the soil physical parameters obtained from multiple borehole locations within the target area.

[0039] The second calculation module is used to calculate the basic safety factor based on the ratio of the critical cyclic stress ratio to the actual cyclic stress ratio.

[0040] The correction module is used to correct the basic safety factor based on the proportion of coarse sand at each borehole location, so as to obtain the corrected safety factor.

[0041] The liquefaction index module is used to obtain the liquefaction index of each borehole location based on the range of the corrected safety factor where the corrected safety factor of each borehole location is located and the liquefaction index corresponding to the range of the corrected safety factor.

[0042] The distribution map module is used to obtain a spatial distribution map of the liquefaction index of the target area based on the liquefaction index of all borehole locations and a preset geological complexity index.

[0043] Thirdly, embodiments of the present invention provide an electronic device, including: at least one processor, at least one memory, and computer program instructions stored in the memory, which, when executed by the processor, implement the method of the first aspect described above.

[0044] Fourthly, embodiments of the present invention provide a storage medium storing computer program instructions, which, when executed by a processor, implement the method of the first aspect described above.

[0045] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0046] By introducing the proportion of coarse sand mass to correct the basic safety factor, the actual contribution of coarse sand content to the soil's drainage capacity and liquefaction resistance can be quantitatively considered on the basis of the traditional liquefaction safety factor calculation method. This makes the corrected safety factor more accurately reflect the true liquefaction risk state of the soil and avoids overly conservative assessments caused by ignoring the influence of coarse sand.

[0047] By converting the modified safety factor into a liquefaction index, continuous safety factor values ​​are mapped to a standardized probability scale between 0 and 1, making the liquefaction risks between different boreholes comparable. At the same time, it provides input data with a unified dimension for subsequent spatial interpolation, which facilitates risk extrapolation for continuous areas.

[0048] By using spatial interpolation to extend the liquefaction index of discrete borehole points into a continuous grid distribution map of the entire site, the shortcomings of traditional single-point discrimination methods that cannot directly generate a planar map of the liquefaction influence range are overcome. This avoids the experience and arbitrariness of manually drawing boundaries and improves the objectivity and repeatability of the influence range delineation.

[0049] By introducing a geological complexity index to weight and correct the interpolation results, the risk assessment of local areas can be adjusted according to the actual non-uniformity of the site (such as lenses, interlayers, faults, or groundwater level fluctuations). This improves the adaptability and accuracy of the liquefaction index distribution map under different geological conditions, making the assessment results in complex sites more reasonable.

[0050] By combining the preset risk zoning thresholds to output a liquefaction impact range prediction layer, the continuous liquefaction index is transformed into a zoning map that can be directly used in engineering. This provides a clear spatial range basis for the targeted deployment of anti-liquefaction measures and reduces the ambiguity in engineering decision-making. Attached Figure Description

[0051] To more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly described below. It should be understood that the following drawings only show some embodiments of the present invention and should not be considered as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort. In the drawings:

[0052] Figure 1 A flowchart illustrating the method for assessing the impact range of sand liquefaction provided by this invention;

[0053] Figure 2 A schematic diagram of the structure of the sand liquefaction impact range assessment device provided by the present invention;

[0054] Figure 3 This is a schematic diagram of the structure of the electronic device provided by the present invention. Detailed Implementation

[0055] To make the objectives, technical solutions, and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the embodiments and accompanying drawings. The illustrative embodiments and descriptions of the present invention are only used to explain the present invention and are not intended to limit the present invention.

[0056] It should be noted that, in this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising..." does not exclude the presence of additional identical elements in the process, method, article, or apparatus that includes said element.

[0057] It should be noted that all actions involving the acquisition of signals, information, or data in this invention are carried out in compliance with the relevant data protection laws and regulations of the locality and with authorization from the owner of the relevant device.

[0058] Example 1

[0059] Please see Figure 1 This invention provides a method for assessing the impact range of sand liquefaction, including:

[0060] S1. Based on the soil physical parameters obtained from multiple borehole locations within the target area, calculate the actual cyclic stress ratio and critical cyclic stress ratio for each borehole location.

[0061] The actual cyclic stress ratio is a value characterizing the ratio of the cyclic shear stress level generated in the soil under seismic loading to the initial effective stress. It can be calculated using the Seed simplified method formula.

[0062] The critical cyclic stress ratio is a measure of the soil's ability to resist liquefaction, and it can be calculated using the NCEER empirical formula.

[0063] For each borehole location, the obtained soil physical parameters may include: the standard penetration test blow count (SPT) measured and corrected by the standard penetration test; the fine particle content, fine sand mass percentage, medium sand mass percentage, and coarse sand mass percentage obtained by sieve analysis; and the peak ground acceleration, magnitude, borehole depth, and unit weight of the soil layer obtained by seismic hazard analysis.

[0064] By calculating the actual cyclic stress ratio and the critical cyclic stress ratio respectively, the intensity of seismic action and the soil's resistance to liquefaction can be quantitatively compared.

[0065] In some implementations, the calculation of the actual cyclic stress at each borehole location includes:

[0066] For each borehole location, the product value is obtained based on the ratio of peak ground acceleration to gravitational acceleration in the soil physical parameters corresponding to the borehole location and a preset coefficient.

[0067] Peak ground acceleration (PGA) is the numerical value describing the maximum ground motion acceleration in the site ground motion parameters, denoted as . The unit is gravitational acceleration. A multiple of. Gravitational acceleration. These are standard physical constants. The preset coefficient is an empirical constant derived from the Seed simplification method, with a value of 0.65. The dimensionless acceleration ratio is obtained by dividing the peak ground acceleration by the gravitational acceleration, and then multiplied by the preset coefficient 0.65 to obtain the product. This product reflects the fundamental contribution of seismic intensity to the cyclic shear stress of the soil.

[0068] The stress ratio is obtained by taking the ratio of the total vertical stress to the effective vertical stress at the borehole depth in the soil physical parameters.

[0069] Among them, the total vertical stress This is the total pressure generated by the weight of the upper soil column at the borehole depth, including contributions from soil particles and pore water. Effective vertical stress. This is the total vertical stress minus the pore water pressure, representing the actual pressure borne by the soil skeleton. Stress ratio This is used to convert the surface acceleration reduction effect to a specified depth, because the shear stress generated at a deeper depth by the same surface acceleration needs to be amplified by this ratio.

[0070] The depth reduction factor is calculated based on the borehole depth using a preset depth reduction factor formula.

[0071] Among them, the depth reduction factor This is a factor that decreases with increasing depth, used to correct for the distribution error of shear stress along depth under the assumption of a rigid body. The preset depth reduction factor formula is a linear relationship: ,in The value represents the borehole depth in meters. This formula is applicable to depths not exceeding 20 meters. The depth reduction factor calculated using this formula reflects the decrease in shear stress with depth in actual soil during earthquakes.

[0072] The actual cyclic stress ratio is obtained by multiplying the product value, the stress ratio, and the depth reduction factor.

[0073] The actual cyclic stress ratio is obtained by multiplying the product value, stress ratio, and depth reduction factor obtained from the above three steps in sequence. Its complete expression is: This ratio represents the depth under design seismic loading. The ratio of the equivalent cyclic shear stress to the initial effective stress experienced by the soil is a key input parameter for subsequent liquefaction risk assessment.

[0074] Through the above step-by-step calculations, the actual cyclic stress ratio at each borehole location can be accurately calculated based on the site seismic parameters and soil profile data, providing seismic input for determining liquefaction risk.

[0075] In some implementations, calculating the critical cyclic stress ratio for each borehole location includes:

[0076] For each borehole location, a first intermediate value is calculated based on the standard penetration blow count in the soil physical parameters corresponding to the borehole location;

[0077] Among them, standard penetration strikes This is the hammer blow count after overlying pressure correction, used to characterize the compaction of sand. The first intermediate value can be... When the standard penetration number is less than 30, it increases as the number of blows increases.

[0078] A second intermediate value is calculated based on the standard penetration number, and the second intermediate value increases linearly with the increase of the standard penetration number.

[0079] The second intermediate value can be This value increases linearly with the standard penetration blow count, reflecting the positive contribution of the penetration blow count to the resistance to liquefaction intensity.

[0080] A third intermediate value is calculated based on the standard penetration number, and the third intermediate value decreases according to a second preset function as the standard penetration number increases.

[0081] The third intermediate value can be It decreases rapidly as the number of standard penetration blows increases, indicating that the contribution of this term gradually weakens when the number of blows is high.

[0082] The first intermediate value and the second intermediate value are added together, and then the third intermediate value is subtracted to obtain the basic critical cyclic stress ratio.

[0083] The first and second intermediate values ​​are summed, then the third intermediate value is subtracted, and finally a constant (such as 1 / 200) is subtracted to obtain the critical cyclic stress ratio of the foundation. This foundation value represents the soil's resistance to liquefaction without considering the influence of fine particle content.

[0084] The incremental correction value is calculated based on the fine particle content in the soil physical parameters.

[0085] Among them, fine particulate content This represents the percentage of particles with a diameter less than 0.075 mm by mass. The incremental correction value is calculated as follows: when the fine particle content is less than the first threshold, the correction value is proportional to the fine particle content; when the fine particle content exceeds the second threshold, the correction value tends to stabilize. Specifically, a piecewise linear function is used: when the fine particle content is between 5% and 20%, it increases linearly with the increase of the fine particle content; when it is greater than 20%, it remains constant.

[0086] The critical cyclic stress ratio is obtained based on the basic critical cyclic stress ratio and the incremental correction value.

[0087] The critical cyclic stress ratio is obtained by adding the basic critical cyclic stress ratio to the incremental correction value. This value comprehensively reflects the liquefaction resistance of sand under the combined effects of standard penetration test blow count and fine particle content. Through the above step-by-step calculation, the critical cyclic stress ratio at each borehole location can be quickly estimated using standard penetration test data obtained from conventional exploration, serving as the basis for subsequent safety factor calculations.

[0088] S2. Calculate the basic safety factor based on the ratio of the critical cyclic stress ratio to the actual cyclic stress ratio;

[0089] The basic safety factor is a dimensionless ratio obtained by dividing the critical cyclic stress ratio by the actual cyclic stress ratio. A coefficient greater than 1 indicates that the soil's resistance to liquefaction is greater than that under seismic action, while a coefficient less than 1 indicates that seismic action exceeds its resistance to liquefaction.

[0090] The basic safety factor is obtained by dividing the critical cyclic stress ratio at each borehole location obtained in S1 by the corresponding actual cyclic stress ratio. The smaller this factor, the higher the risk of liquefaction. By calculating the basic safety factor, the liquefaction tendency at each borehole location can be preliminarily determined.

[0091] S3. Based on the proportion of coarse sand at each drilling location, the basic safety factor is corrected to obtain the corrected safety factor;

[0092] The coarse sand mass ratio refers to the percentage of the total mass of the soil sample whose particle size falls within the third preset particle size range.

[0093] The corrected safety factor is the value obtained by multiplying the basic safety factor by the particle size distribution correction factor.

[0094] Obtain the mass percentages of fine sand, medium sand, and coarse sand from sieve analysis. Calculate the particle size distribution correction factor based on the coarse sand mass percentage; this factor increases with increasing coarse sand mass percentage. Multiply the basic safety factor by this correction factor to obtain the corrected safety factor. A higher coarse sand percentage results in a larger corrected safety factor, indicating that the liquefaction risk at the borehole location is considered lower than the basic calculation value after correction. By incorporating the coarse sand mass percentage for correction, the assessment results better reflect the actual physical behavior of the soil.

[0095] In some embodiments, the step of correcting the basic safety factor based on the proportion of coarse sand at each borehole location to obtain a corrected safety factor includes:

[0096] The total sand mass percentage is calculated based on the mass percentages of fine sand, medium sand, and coarse sand.

[0097] Among them, the proportion of fine sand in mass This represents the percentage of the total mass of the soil sample containing particles with a diameter between 0.075 mm and 0.25 mm. Medium sand mass percentage. This represents the percentage by mass of particles with a diameter between 0.25 mm and 0.5 mm. Coarse sand mass percentage. This represents the percentage by mass of particles with a diameter between 0.5 mm and 2.0 mm. Adding these three percentages together gives the total sand mass percentage. This value represents the percentage of the total mass of all sand particles (excluding fine particles).

[0098] The correction increment is determined based on the ratio of the mass percentage of coarse sand to the mass percentage of total sand.

[0099] Among them, the ratio of the mass percentage of coarse sand to the mass percentage of total sand is calculated, i.e. This ratio reflects the relative abundance of coarse sand in the sand grains. The correction increment is this ratio multiplied by a preset second coefficient, for example, 0.5, i.e., increment = When the proportion of coarse sand is high, the correction increment is large; when there is no coarse sand, the correction increment is zero.

[0100] Based on the aforementioned correction increment, the particle size distribution correction coefficient is obtained;

[0101] Among them, the particle size distribution correction coefficient It equals 1 plus the above correction increment, that is The higher the coarse sand content, the larger the correction coefficient, reflecting the physical law that coarse-grained soil has improved drainage capacity and enhanced resistance to liquefaction.

[0102] The basic safety factor is corrected based on the particle size distribution correction factor to obtain the corrected safety factor.

[0103] Among them, the basic safety factor Multiply by particle size distribution correction factor The corrected safety factor is obtained. The basic safety factor is the ratio of the critical cyclic stress ratio to the actual cyclic stress ratio. Multiplying it by a correction factor greater than 1 increases the safety factor, meaning that after considering the beneficial effects of coarse sand, the actual safety margin is improved and the liquefaction risk is reduced. The corrected safety factor will serve as the basis for subsequent calculations of the liquefaction index.

[0104] Through the above correction steps, boreholes with higher coarse sand content are considered safer (i.e., have a higher safety factor) under the same standard penetration test (SPT) blow count and fine grain content conditions, thus reflecting the contribution of coarse sand to resist liquefaction in the liquefaction index and spatial distribution. This correction is based on engineering experience and is applicable to particle size classification indicators available in conventional exploration data.

[0105] S4. Based on the range of corrected safety factors for each borehole location and the liquefaction index corresponding to the range of corrected safety factors, obtain the liquefaction index for each borehole location.

[0106] The liquefaction index is a dimensionless value that maps the modified safety factor to a preset range, used to characterize the probability or degree of liquefaction. The higher the index, the higher the risk of liquefaction.

[0107] The corrected safety factor for each borehole location is compared with a pre-set first and second preset threshold. When the corrected safety factor is less than or equal to the first preset threshold, the liquefaction index is set to the first preset value, indicating a high liquefaction risk. When the corrected safety factor is greater than or equal to the second preset threshold, the liquefaction index is set to the second preset value, indicating a low liquefaction risk. When the corrected safety factor falls between the first and second preset thresholds, the liquefaction index is calculated using a monotonic function interpolation, with the corrected safety factor as the independent variable and the liquefaction index as the dependent variable. By converting continuous corrected safety factors into discrete, graded liquefaction indices, subsequent spatial interpolation and risk zoning are facilitated.

[0108] In some embodiments, obtaining the liquefaction index for each borehole location based on the range of corrected safety factors for each borehole location and the liquefaction index corresponding to the corrected safety factor range includes:

[0109] When the corrected safety factor is less than or equal to the first preset threshold, the liquefaction index takes the first preset value;

[0110] Among them, the corrected safety factor This is the safety factor after adjusting for the proportion of coarse sand; a higher value indicates greater safety. The first preset threshold is the critical value for judging high liquefaction risk and can be less than 1, such as 0.8. When the adjusted safety factor is lower than this threshold, it indicates that the soil's resistance to liquefaction is severely insufficient, and the risk of liquefaction is extremely high. In this case, the liquefaction index is taken as the first preset value, such as 0.95, representing a high probability of liquefaction.

[0111] When the safety factor is greater than or equal to the second preset threshold, the liquefaction index takes the second preset value, which is less than the first preset value.

[0112] The second preset threshold is a critical value for judging low liquefaction risk, and can be greater than 1, such as 1.2. When the corrected safety factor is higher than this threshold, it indicates that the soil has sufficient resistance to liquefaction and the liquefaction risk is extremely low. At this time, the liquefaction index is taken as the second preset value, such as 0.05, representing a low liquefaction probability.

[0113] When the corrected safety factor is between the first preset threshold and the second preset threshold, a proportional value is obtained based on the first difference between the corrected safety factor and the first preset threshold, the second difference between the second preset threshold and the first preset threshold, and the ratio of the first difference and the second difference.

[0114] Specifically, when the corrected safety factor falls within the middle range, a linear interpolation method is used to continuously map the safety factor to the liquefaction index. Specifically, the difference between the corrected safety factor and a first preset threshold is calculated as the first difference value. The difference between a second preset threshold and the first preset threshold is calculated as the second difference value. The first difference value is divided by the second difference value to obtain a proportional value. This proportional value is between 0 and 1; when the corrected safety factor is close to the first preset threshold, the proportional value is close to 0, and when it is close to the second preset threshold, the proportional value is close to 1.

[0115] The liquefaction index is obtained by linearly interpolating the ratio value between the first preset value and the second preset value.

[0116] The liquefaction index is obtained by multiplying the proportional value by the difference between the second and first preset values, and then adding the first preset value. That is, liquefaction index = first preset value + proportional value × (second preset value - first preset value). Since the second preset value is less than the first preset value, this interpolation result decreases as the safety factor increases. This linear interpolation allows for a smooth transition between high and low risk levels.

[0117] The piecewise function described above converts the continuous safety factor into a liquefaction index between a first preset value and a second preset value. This index quantifies the liquefaction probability or risk level at each borehole location, facilitating subsequent spatial interpolation and risk zoning.

[0118] S5. Based on the liquefaction index of all borehole locations and the preset geological complexity index, obtain the spatial distribution map of the liquefaction index of the target area.

[0119] The geological complexity index is a quantitative parameter reflecting the degree of drastic changes in site strata, the range of changes in groundwater level depth, and the presence of faults or lenses. The geological complexity index is pre-set based on the degree of drastic changes in strata, the range of changes in groundwater level depth, and the presence of faults or lenses. For example, a first preset index value is used for simple sites, and a second preset index value is used for complex sites, where:

[0120] The criteria for classifying a site as simple are as follows: the soil layers within the site are uniformly and continuously distributed in the depth direction, without obvious pinch-outs or lenses; the main soil layer is a single layer or a sand layer with stable bedding, with minimal variation in layer thickness; the variation in groundwater level depth within the site is less than a preset depth variation threshold, and there are no local groundwater level anomalies (such as confined water or perched water zones); no active faults pass through the site, or the faults are located at a distance greater than a preset safety distance from the site boundary. When all the above conditions are met, the site is classified as a simple site, and the geological complexity index is set to the first preset index value. The criteria for classifying a site as complex are as follows: the site contains at least one of the following characteristics: uneven soil layer distribution, with the presence of lenses, interlayers, pinch-outs, or drastic phase transitions; the variation in groundwater level depth exceeds a preset depth variation threshold, or there are local groundwater level anomalies; an active fault passes through the site, or the faults are located at a distance less than a preset safety distance from the site boundary. When any one of these conditions is met, the site is classified as a complex site, and the geological complexity index is set to the second preset index value. In specific assessments, stratigraphic comparisons are conducted using stratigraphic profiles revealed by boreholes. The top and bottom depth variations of the same stratigraphic level in each borehole are statistically analyzed, and the standard deviation or coefficient of variation of the stratigraphic thickness is calculated. When the coefficient of variation exceeds a preset threshold, the stratigraphic level is deemed complex. Groundwater level changes are evaluated by the ratio of the range of stable water level measurements from each borehole to the average burial depth; if this ratio exceeds a preset threshold, the groundwater level is deemed complex. Fault activity is determined based on regional geological data and field investigations.

[0121] The liquefaction index spatial distribution map is a raster map obtained by extending the liquefaction index of discrete borehole points to the entire continuous region through spatial interpolation.

[0122] Spatial interpolation extends the evaluation results of discrete borehole points to the entire site, overcoming the limitation of traditional single-point discrimination in generating continuous distribution maps. Multiplying by a geological complexity index allows for weighted correction of local areas under complex geological conditions, making the distribution map more consistent with actual geological characteristics.

[0123] In some embodiments, obtaining the spatial distribution map of the liquefaction index of the target area based on the liquefaction index of all borehole locations and a preset geological complexity index includes:

[0124] Using the planar coordinates and liquefaction index of all borehole locations as known points, a preset interpolation algorithm is used to perform spatial interpolation to obtain the interpolated liquefaction index of each grid node.

[0125] The planar coordinates of the borehole locations are the two-dimensional coordinates of each borehole's projection onto the Earth's surface, obtained through measurement, and are denoted as... The preset interpolation algorithm can be a spatial statistical method used to estimate the value at an unknown location based on the values ​​at discrete known points. This includes inverse distance weighting and Kriging. The principle of spatial interpolation is that the closer a location is to a known point, the closer its estimated value is to that known point; the farther away it is, the smaller the influence. The entire target area is divided into regular grids, each grid cell having a central node. For each grid node, the interpolated liquefaction index is calculated by weighted averaging based on the liquefaction indices of known borehole points around it. This process yields a continuous probability surface covering the entire region, where each grid location corresponds to an interpolated liquefaction risk estimate.

[0126] The final liquefaction index raster value is obtained based on the interpolated liquefaction index of each raster node and the preset geological complexity index.

[0127] Among them, the geological complexity index It is a pre-defined parameter that varies with spatial location, reflecting the degree of heterogeneity of the strata at that location. In simple sites, Take the first preset index value, such as 1.0, which means the interpolation result can be used directly. In complex sites, A second preset index value, such as 1.5, is chosen to indicate that these complex factors increase the uncertainty or local differences in liquefaction risk, thus requiring amplification and weighting of the interpolation results. Final liquefaction index raster value. The result is obtained by multiplying the interpolated liquefaction index by the geological complexity index: After multiplication, if the result exceeds 1.0, it is rounded to 1.0. This allows for a more conservative assessment of risk in complex geological areas, where relatively high risk estimates are amplified to reflect uncertainties; in simpler areas, the interpolation results are used directly. The resulting raster map is the spatial distribution map of the liquefaction index for the target area, and the value of each raster cell can be directly used for subsequent risk zoning and impact delineation. By extending the liquefaction assessment results of discrete borehole points to a continuous distribution across the entire site, and introducing a geological complexity index for weighted correction of local areas, the spatial distribution map better reflects the actual geological characteristics, providing an intuitive and quantitative basis for engineering decisions.

[0128] The following is a specific example to further illustrate this embodiment.

[0129] Drill holes in a grid pattern within the target area, with a spacing of 50m between adjacent holes. Each hole should be drilled to a depth at least 5m below the expected liquefaction depth. Standard penetration test (SPT) and sieve test should be performed in each hole.

[0130] Standard penetration test: Record the measured number of hammer blows. (hit / 30cm), corrected to .

[0131] Sieve analysis: Determine the content of fine particles with a diameter less than 0.075 mm. (%); and the mass percentage of particles in the three particle size ranges of 0.075mm~0.25mm, 0.25mm~0.5mm, and 0.5mm~2.0mm are respectively recorded as the mass percentage of fine sand. Medium sand quality ratio Coarse sand mass ratio ,and .

[0132] In this embodiment, a total of 5 boreholes were drilled, and their coordinates and test data are as follows:

[0133] Drill hole ZK1: coordinates (0m, 0m); ; ; , , .

[0134] Drill hole ZK2: coordinates (100m, 0m); ; ; , , .

[0135] Drill hole ZK3: coordinates (100m, 100m); ; ; , , .

[0136] Drill hole ZK4: coordinates (0m, 100m); ; ; , , .

[0137] Drill hole ZK5: coordinates (50m, 50m); ; ; , , .

[0138] The groundwater level at the site is 2.0m deep. The unit weight of the soil above water is taken as 18kN / m³, and the unit weight underwater is taken as 9.8kN / m³. The depth of each borehole is calculated as 7m.

[0139] Calculate the basic safety factor for each borehole location.

[0140] (1) Calculate the actual cyclic stress ratio

[0141] Using the Seed simplified formula:

[0142] ;

[0143] in:

[0144] The peak ground acceleration is 0.20g for this site.

[0145] The acceleration due to gravity is taken as 9.81 m / s².

[0146] For depth Total vertical stress at the location (kPa);

[0147] For depth Effective vertical stress (kPa) at the location;

[0148] For depth reduction factor, take ( m).

[0149] depth m, calculate the total vertical stress:

[0150] ;

[0151] Effective vertical stress:

[0152] ;

[0153] Depth reduction factor:

[0154] ;

[0155] have to:

[0156] ;

[0157] All boreholes are at the same depth, therefore Both are 0.290.

[0158] (2) Calculate the critical cyclic stress ratio

[0159] The NCEER empirical formula is used. When When the critical cyclic stress ratio of the foundation is reached, it is calculated using the following formula:

[0160] ;

[0161] like Then, incremental correction:

[0162] ;

[0163] Critical cyclic stress ratio:

[0164] ;

[0165] Calculate each borehole separately:

[0166] Drilling ZK1: , .

[0167] ;

[0168] , Therefore .

[0169] Drilling ZK2: , .

[0170] ;

[0171] , .

[0172] Drilling ZK3: , .

[0173] ;

[0174] , .

[0175] Drilling ZK4: , .

[0176] ;

[0177] , .

[0178] Drilling ZK5: , .

[0179] ;

[0180] , .

[0181] (3) Calculate the basic safety factor

[0182] Basic safety factor .

[0183] ZK1: ;

[0184] ZK2: ;

[0185] ZK3: ;

[0186] ZK4: ;

[0187] ZK5: ;

[0188] Correcting the safety factor using particle classification

[0189] Introducing particle size distribution correction factor Considering the effect of coarse sand mass ratio on improving liquefaction resistance:

[0190] ;

[0191] Because of each drill hole The calculation is as follows:

[0192] ZK1: , ;

[0193] ZK2: , ;

[0194] ZK3: , ;

[0195] ZK4: , ;

[0196] ZK5: , .

[0197] Corrected safety factor :

[0198] ZK1: ;

[0199] ZK2: ;

[0200] ZK3: ;

[0201] ZK4: ;

[0202] ZK5: .

[0203] Step S103: Calculate the liquefaction index at each borehole location.

[0204] Define liquefaction index Values ​​between 0 and 1, with larger values ​​indicating a higher risk of liquefaction. A piecewise linear function is used.

[0205] when hour, (High risk);

[0206] when hour, (Low risk);

[0207] when Linear interpolation:

[0208] ;

[0209] calculate:

[0210] ZK1: ;

[0211] ZK2: ;

[0212] ZK3: ;

[0213] ZK4: ;

[0214] ZK5: .

[0215] Step S104: Spatial interpolation of geological complexity data

[0216] The site has homogeneous strata and a stable groundwater level, classifying it as a simple site with a geological complexity index of [missing information]. For complex sites, take 1.5.

[0217] Using the planar coordinates of all boreholes and liquefaction index As known points, spatial interpolation was performed using ordinary kriging, with a spherical model chosen for the variogram. The interpolation grid resolution was 50m × 50m. The initial liquefaction index for each grid node was obtained through interpolation, and then multiplied by... The final raster value is obtained. Because... The final value is equal to the interpolation result.

[0218] The following estimation uses the inverse distance weighted method, with a power exponent of 2, to obtain the interpolation liquefaction index for some raster nodes:

[0219] Node (25,25): Close to ZK1 and ZK5, interpolation value is approximately 0.92;

[0220] Node (75,25): Close to ZK2 and ZK5, interpolation value is approximately 0.93;

[0221] Node (25,75): Close to ZK4, interpolation approximately 0.52;

[0222] Node (75,75): Close to ZK3, interpolation value is approximately 0.85.

[0223] Output liquefaction influence range prediction layer

[0224] Based on the liquefaction index raster chart, set risk zoning thresholds:

[0225] Liquefaction Zone (High Risk):

[0226] Potential liquefaction zone (medium risk):

[0227] Non-liquefiable area (low risk):

[0228] Based on the above interpolation results, a liquefaction influence range map is generated:

[0229] Liquefaction Zone: Covers the area around ZK1, ZK2, ZK5 and the (75,25) and (75,75) node areas, forming a large area.

[0230] Possible liquefaction zone: Only appears in a small area near node (25,75), with an area of ​​approximately 2500 m².

[0231] Non-liquefiable zone: A circular area with a radius of approximately 20m centered on ZK4.

[0232] The final output is a GIS-formatted map of the liquefaction impact area, overlaid with borehole locations and topographic base maps, and accompanied by a table of raster probability values.

[0233] Example 2

[0234] Please see Figure 2 This invention provides a device for assessing the impact range of sand liquefaction, comprising:

[0235] The first calculation module is used to calculate the actual cyclic stress ratio and critical cyclic stress ratio for each borehole location based on the soil physical parameters obtained from multiple borehole locations within the target area.

[0236] The second calculation module is used to calculate the basic safety factor based on the ratio of the critical cyclic stress ratio to the actual cyclic stress ratio.

[0237] The correction module is used to correct the basic safety factor based on the proportion of coarse sand at each borehole location, so as to obtain the corrected safety factor.

[0238] The liquefaction index module is used to obtain the liquefaction index of each borehole location based on the range of the corrected safety factor where the corrected safety factor of each borehole location is located and the liquefaction index corresponding to the range of the corrected safety factor.

[0239] The distribution map module is used to obtain a spatial distribution map of the liquefaction index of the target area based on the liquefaction index of all borehole locations and a preset geological complexity index.

[0240] It should be noted that each module and unit in the sand liquefaction impact range assessment device in this embodiment corresponds one-to-one with each step in the sand liquefaction impact range assessment method in the aforementioned embodiment. Therefore, the specific implementation of this embodiment can refer to the implementation of the aforementioned sand liquefaction impact range assessment method, and will not be repeated here.

[0241] Example 3

[0242] Please see Figure 3 This embodiment provides an electronic device, including at least one processor 301 and a memory 302. Optionally, the device further includes a communication component 303. The processor 301, memory 302, and communication component 303 are connected via a bus 304.

[0243] In a specific implementation, at least one processor 301 executes computer execution instructions stored in memory 302, causing at least one processor 301 to perform the above-described method.

[0244] The specific implementation process of processor 301 can be found in the above method embodiments, and its implementation principle and technical effect are similar. It will not be repeated here.

[0245] In the above embodiments, it should be understood that the processor can be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), etc. The general-purpose processor can be a microprocessor or any conventional processor. The steps of the method disclosed in this invention can be directly implemented by a hardware processor, or implemented by a combination of hardware and software modules within the processor.

[0246] The memory may include random access memory (RAM) and may also include non-volatile memory (NVM), such as at least one disk storage device.

[0247] The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. Buses can be categorized as address buses, data buses, control buses, etc. For ease of illustration, the buses shown in the accompanying drawings are not limited to a single bus or a single type of bus.

[0248] This application also provides a computer program product, including a computer program that, when executed by a processor, implements the above-described method.

[0249] This application also provides a computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, implement the above-described method.

[0250] The aforementioned readable storage medium can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as 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. The readable storage medium can be any available medium accessible to a general-purpose or special-purpose computer.

[0251] An exemplary readable storage medium is coupled to a processor, enabling the processor to read information from and write information to the readable storage medium. Of course, the readable storage medium can also be a component of the processor. The processor and the readable storage medium can reside in an Application Specific Integrated Circuit (ASIC). Alternatively, the processor and the readable storage medium can exist as discrete components in the device.

[0252] The division of units is merely a logical functional division; in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be indirect coupling or communication connection through some interfaces, devices, or units, and may be electrical, mechanical, or other forms.

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

[0254] In addition, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0255] If a function is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0256] Those skilled in the art will understand that all or part of the steps of the above-described method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When executed, the program performs the steps of the above-described method embodiments; and the aforementioned storage medium includes various media capable of storing program code, such as ROM, RAM, magnetic disks, or optical disks.

[0257] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for assessing the impact range of sand liquefaction, characterized in that, include: Based on the soil physical parameters obtained from multiple borehole locations within the target area, calculate the actual cyclic stress ratio and critical cyclic stress ratio for each borehole location; The basic safety factor is calculated based on the ratio of the critical cyclic stress ratio to the actual cyclic stress ratio. The basic safety factor is corrected based on the proportion of coarse sand at each borehole location to obtain the corrected safety factor. The liquefaction index for each borehole location is obtained based on the range of corrected safety factors in which the corrected safety factor for each borehole location falls and the liquefaction index corresponding to the range of corrected safety factors. Based on the liquefaction index of all borehole locations and the preset geological complexity index, a spatial distribution map of the liquefaction index of the target area is obtained.

2. The method for assessing the impact range of sand liquefaction according to claim 1, characterized in that, Calculate the actual cyclic stress at each borehole location, including: For each borehole location, the product value is obtained based on the ratio of peak ground acceleration to gravitational acceleration in the soil physical parameters corresponding to the borehole location and a preset coefficient. The stress ratio is obtained by taking the ratio of the total vertical stress to the effective vertical stress at the borehole depth in the soil physical parameters. The depth reduction factor is calculated based on the borehole depth using a preset depth reduction factor formula. The actual cyclic stress ratio is obtained by multiplying the product value, the stress ratio, and the depth reduction factor.

3. The method for assessing the impact range of sand liquefaction according to claim 1, characterized in that, Calculate the critical cyclic stress ratio for each borehole location, including: For each borehole location, a first intermediate value is calculated based on the standard penetration blow count in the soil physical parameters corresponding to the borehole location; A second intermediate value is calculated based on the standard penetration number, and the second intermediate value increases linearly with the increase of the standard penetration number. A third intermediate value is calculated based on the standard penetration number, and the third intermediate value decreases according to a second preset function as the standard penetration number increases. The first intermediate value and the second intermediate value are added together, and then the third intermediate value is subtracted to obtain the basic critical cyclic stress ratio. The incremental correction value is calculated based on the fine particle content in the soil physical parameters. The critical cyclic stress ratio is obtained based on the basic critical cyclic stress ratio and the incremental correction value.

4. The method for assessing the impact range of sand liquefaction according to claim 1, characterized in that, The basic safety factor is corrected based on the proportion of coarse sand at each borehole location to obtain a corrected safety factor, including: The total sand mass percentage is calculated based on the mass percentages of fine sand, medium sand, and coarse sand. The correction increment is determined based on the ratio of the mass percentage of coarse sand to the mass percentage of total sand. Based on the aforementioned correction increment, the particle size distribution correction coefficient is obtained; The basic safety factor is corrected based on the particle size distribution correction factor to obtain the corrected safety factor.

5. The method for assessing the impact range of sand liquefaction according to claim 1, characterized in that, The step of obtaining the liquefaction index for each borehole location based on the range of corrected safety factors for each borehole location and the liquefaction index corresponding to the range of corrected safety factors includes: When the corrected safety factor is less than or equal to the first preset threshold, the liquefaction index takes the first preset value; When the safety factor is greater than or equal to the second preset threshold, the liquefaction index takes the second preset value, which is less than the first preset value. When the corrected safety factor is between the first preset threshold and the second preset threshold, a proportional value is obtained based on the first difference between the corrected safety factor and the first preset threshold, the second difference between the second preset threshold and the first preset threshold, and the ratio of the first difference and the second difference. The liquefaction index is obtained by linearly interpolating the ratio value between the first preset value and the second preset value.

6. The method for assessing the impact range of sand liquefaction according to claim 1, characterized in that, The geological complexity index is preset based on the degree of drastic changes in strata, the range of changes in groundwater depth, and the presence or absence of faults or lenses.

7. The method for assessing the impact range of sand liquefaction according to claim 1, characterized in that, The step of obtaining a spatial distribution map of the liquefaction index of the target area based on the liquefaction index of all borehole locations and a preset geological complexity index includes: Using the planar coordinates and liquefaction index of all borehole locations as known points, a preset interpolation algorithm is used to perform spatial interpolation to obtain the interpolated liquefaction index of each grid node. The final liquefaction index raster value is obtained based on the interpolated liquefaction index of each raster node and the preset geological complexity index.

8. A device for assessing the impact range of sand liquefaction, characterized in that, include: The first calculation module is used to calculate the actual cyclic stress ratio and critical cyclic stress ratio for each borehole location based on the soil physical parameters obtained from multiple borehole locations within the target area. The second calculation module is used to calculate the basic safety factor based on the ratio of the critical cyclic stress ratio to the actual cyclic stress ratio. The correction module is used to correct the basic safety factor based on the proportion of coarse sand at each borehole location, so as to obtain the corrected safety factor. The liquefaction index module is used to obtain the liquefaction index of each borehole location based on the range of the corrected safety factor where the corrected safety factor of each borehole location is located and the liquefaction index corresponding to the range of the corrected safety factor. The distribution map module is used to obtain a spatial distribution map of the liquefaction index of the target area based on the liquefaction index of all borehole locations and a preset geological complexity index.

9. An electronic device, characterized in that, include: At least one processor, at least one memory, and computer program instructions stored in the memory, which, when executed by the processor, implement the method as described in any one of claims 1-7.

10. A computer-readable storage medium having computer program instructions stored thereon, characterized in that, The method as described in any one of claims 1-7 is implemented when the computer program instructions are executed by the processor.

Images