Disclosure of Invention
Technical problem
The invention provides a method for calculating the proportionality coefficient of a soil level reaction coefficient and an identification system for solving the problem that retrieval is inconvenient after a ground service work order is filed in airport operation; the surface wave is a geophysical prospecting method, has the characteristics of no damage, strong field adaptability and the like, the wave speed of characteristic parameters of the surface wave is closely related to the physical and mechanical properties of a soil body, and the surface wave is obtained under the conditions that the in-situ stress state and the natural water content are unchanged and the primary structure is not disturbed.
Technical scheme
The first purpose of the invention is to provide a method for calculating a proportionality coefficient of a soil level reaction force coefficient, which comprises the following steps:
s1, data acquisition: carrying out field surface wave exploration and collecting transient surface wave data;
s2, data processing: obtaining the wavelet velocity data within the calculated depth range by processing the transient surface wave data;
s3, dividing soil layers;
s4, determining soil layer influence factors a, b, c, d and e in a calculated depth range;
s5, determining a depth influence coefficient lambda;
s6, determining an area experience correction coefficient psi;
s7, according to the formula m = ψ λ (av)b +cv2 -dv + e) calculating the proportionality coefficient of the horizontal reaction coefficient of the soilWherein:
m is a proportional coefficient of a horizontal reaction force coefficient of the soil;
psi is a regional experience correction coefficient;
and v is the soil layer surface wave speed in the calculated depth range, and the unit is m/s.
Preferably, S1 is specifically: the method comprises the steps that a seismometer is used for carrying out surface wave exploration, after the geophones are arranged, the geophones are fixed and static, a shot point moves from one side of the geophone arrangement to a position at a certain horizontal distance away from the geophone arrangement, the moving track of the shot point is parallel to the geophone arrangement, the moving distance of the shot point every time is equal to the geophone distance, transient surface wave data are collected once every time the shot point moves once until the shot point moves to the position at the other side of the geophone arrangement at the same horizontal distance away from the geophone arrangement.
Preferably, S2 is specifically: and taking each detector as a common central point, extracting all transient surface wave data taking the detector as the common central point for superposition, and obtaining surface wave velocity distribution in a calculated depth range through dispersion curve extraction and velocity inversion.
Preferably, S3 is specifically: and dividing soil layers according to the change of the surface wave velocity along with the depth, and determining the type of each layer of soil and the surface wave velocity v in the calculated depth range.
Preferably, in S4: the soil layer influence factors a, b, c, d and e in the calculated depth range are determined by the following method: when the soil layer is clay, a =0.003, b =1.752, c =0.002, d =0.295, e =14.6; when the soil layer is sandy soil, a =0.0023, b =1.236, c =0.063, d =0.303, e =10.1; when the soil layer is sandy soil, a =0.0015, b =0.835, c =0.125, d =0.516, e =16.7.
Preferably, in S5: by the formula λ =1.08 (σ)k /σpk ) Calculating a depth influence coefficient lambda; wherein σk Calculating the standard value of the additional stress in the soil at the point under the action of load, wherein the unit is kPa; sigmapk Calculating the vertical total stress generated by the self weight of the soil at a point, wherein the unit is kPa; wherein, when lambda<At 0.1, λ =0.1 is taken.
Preferably, in S6: the region empirical correction coefficient psi is determined according to horizontal displacement data of the region foundation pit retaining member and experience, and psi =1.0 is taken in inexperienced regions.
A second object of the present invention is to provide a system for calculating a proportionality coefficient of a soil level reaction force coefficient, comprising:
a data acquisition module: carrying out field surface wave exploration and collecting transient surface wave data;
a data processing module: obtaining the wavelet velocity data within the calculated depth range by processing the transient surface wave data;
the soil layer dividing module: dividing soil layers;
a parameter determination module: determining soil layer influence factors a, b, c, d and e in a calculated depth range; determining a depth influence coefficient lambda; determining an area experience correction coefficient psi;
a calculation module: receiving data of the data acquisition module, the data processing module, the soil layer division module and the parameter determination module according to the formula m = ψ λ (av)b +cv2 Dv + e) calculating a proportionality coefficient of the horizontal reaction force coefficient of the soil, wherein:
m is a proportional coefficient of a horizontal reaction force coefficient of the soil;
psi is a regional experience correction coefficient;
and v is the soil layer surface wave speed in the calculated depth range, and the unit is m/s.
A third object of the present invention is to provide a computer program for implementing the above method of calculating a proportionality coefficient of a soil level reaction force coefficient.
A fourth object of the present invention is to provide an information data processing terminal that realizes the above-described method of calculating a proportionality coefficient of a soil level reaction force coefficient.
A fifth object of the present patent is to provide a computer-readable storage medium comprising instructions which, when run on a computer, cause the computer to perform the above-described method of calculating a proportionality coefficient of soil level reaction force coefficients.
The invention has the advantages and positive effects that:
the method for calculating the proportionality coefficient of the horizontal reaction coefficient of the soil based on the surface wave exploration technology directly uses the transient surface wave test result to calculate the proportionality coefficient of the horizontal reaction coefficient of the soil, the calculation process is simple and quick, the parameters required by calculation are all from transient surface wave exploration, drilling, horizontal loading, a cross plate shearing test and the like are not needed, the inaccuracy of the test result caused by soil disturbance in drilling sampling and indoor geotechnical tests can be avoided, the defects of long horizontal loading test period, high cost and limited application range of the cross plate shearing test are overcome, and the method has obvious advantages particularly in areas with dense buildings, large human interference factors and difficult implementation of the drilling, horizontal loading and cross plate shearing tests in cities and towns. The method has strong field adaptability, can be used for different types of soil bodies such as clay, sandy soil, gravelly soil and the like, is nondestructive investigation, has small damage to the field, simple used equipment, low cost and high efficiency, still has higher calculation precision under the condition of not using empirical correction coefficients, and has larger popularization and use values and wide application prospects.
Detailed Description
In order to further understand the contents, features and effects of the present invention, the following embodiments are illustrated and described in detail with reference to the accompanying drawings.
Referring to fig. 1 of the drawings, a drawing,
a method for calculating a proportionality coefficient of a soil level reaction coefficient comprises the following steps:
step 1: carrying out field surface wave exploration and collecting transient surface wave data;
in the step 1, a seismometer is used for carrying out surface wave exploration, after the geophones are arranged, the geophones are fixed, shot points start to move from one side of the geophone arrangement at a certain horizontal distance away from the geophone arrangement, the moving track of the shot points is parallel to the geophone arrangement, the moving distance of the shot points each time is equal to the geophone distance (namely, the channel distance), and transient surface wave data are collected once each time the shot points move, until the shot points move to the other side of the geophone arrangement at the same horizontal distance away from the geophone arrangement.
Step 2: and (5) carrying out transient surface wave data processing to obtain the surface wave velocity distribution in the calculated depth range.
In step 2, each detector is taken as a common central point, all transient surface wave data taking the detector as the common central point are extracted and superposed, and surface wave velocity distribution in a calculated depth range is obtained through frequency dispersion curve extraction and velocity inversion.
And step 3: dividing soil layers:
in step 3, dividing soil layers according to the change of the surface wave velocity along with the depth, and determining the type and the surface wave velocity v of each layer of soil in the calculated depth range.
And 4, step 4: determining influence factors a, b, c, d and e of each soil layer in a calculation depth range;
in step 4, the soil layer influence factors a, b, c, d and e in the calculated depth range are determined by the following method: when the soil layer is clay, a =0.003, b =1.752, c =0.002, d =0.295, e =14.6; when the soil layer is sandy soil, a =0.0023, b =1.236, c =0.063, d =0.303, e =10.1; when the soil layer is sandy soil, a =0.0015, b =0.835, c =0.125, d =0.516, e =16.7.
And 5: determining a depth influence coefficient lambda;
in step 5, by the formula λ =1.08 (σ)k /σpk ) Calculating a depth influence coefficient lambda; wherein σk Calculating the standard value of the additional stress in the soil at the point under the action of load, wherein the unit is kPa; sigmapk Calculating the vertical total stress generated by the self weight of the soil at a point, wherein the unit is kPa; wherein, when lambda<At 0.1, λ =0.1 is taken.
Step 6: determining a regional empirical correction factor psi
In step 6, the area experience correction coefficient psi is determined according to the horizontal displacement data of the area foundation pit retaining member and experience, and psi =1.0 is taken in inexperienced areas.
And 7: according to the formula m = ψ λ (aV)b +cV2 dV + e) calculating a proportionality coefficient of the horizontal reaction force coefficient of the soil, wherein:
m is a proportional coefficient of a horizontal reaction force coefficient of the soil;
psi is a regional experience correction coefficient;
λ is a depth influence coefficient;
a. b, c, d and e are soil layer influence factors in the calculated depth range;
v is the soil layer surface wave speed in the calculated depth range, and the unit is m/s; .
A system for calculating a proportionality coefficient of a soil level reaction force coefficient, comprising:
a data acquisition module: carrying out field surface wave exploration and collecting transient surface wave data; the method specifically comprises the following steps: the method comprises the following steps of utilizing a seismograph to conduct field surface wave exploration and collecting transient surface wave data, wherein the specific process is as follows: after the detectors are arranged, the detectors are fixed, the shot points start to move from one side of the arrangement of the detectors to a position with a certain horizontal distance away from the arrangement of the detectors, the moving track of the shot points is parallel to the arrangement of the detectors, the moving distance of the shot points at each time is equal to the distance between the detectors (namely, the track distance), and transient surface wave data are collected once when the shot points move once until the shot points move to the position with the same horizontal distance away from the arrangement of the detectors at the other side of the arrangement of the detectors;
a data processing module: obtaining the wavelet velocity data within the calculated depth range by processing the transient surface wave data; the method specifically comprises the following steps: processing the transient surface wave data acquired by the data acquisition module, taking each detector as a common central point, extracting all transient surface wave data taking the detector as the common central point for superposition, and obtaining the distribution of surface wave velocity along with depth through dispersion curve extraction and velocity inversion;
the soil layer dividing module: dividing soil layers; and dividing soil layers according to the change of the surface wave velocity along with the depth, and determining the type and the surface wave velocity v of each layer of soil in the calculated depth range.
A parameter determination module: determining soil layer influence factors a, b, c, d and e in a calculated depth range; determining a depth influence coefficient lambda; determining an area experience correction coefficient psi; preferably: the method for calculating the soil layer influence factors a, b, c, d and e in the depth range is determined by the following method: when the soil layer is clay, a =0.003, b =1.752, c =0.002, d =0.295, e =14.6; when the soil layer is sandy soil, a =0.0023, b =1.236, c =0.063, d =0.303, e =10.1; when the soil layer is sandy soil, a =0.0015, b =0.835, c =0.125,d =0.516,e =16.7. By the formula λ =1.08 (σ)k /σpk ) Calculating a depth influence coefficient lambda; wherein σk Calculating the standard value of the additional stress in the soil at the point under the action of load, wherein the unit is kPa; sigmapk Calculating the vertical total stress generated by the self weight of the soil at a point, wherein the unit is kPa; wherein, when lambda<At 0.1, λ =0.1 is taken. Determining an empirical region correction coefficient psi according to horizontal displacement data of the region foundation pit retaining member and experience, and taking psi =1.0 in an inexperienced region
A calculation module: receiving data of the data acquisition module, the data processing module, the soil layer division module and the parameter determination module according to the formula m = ψ λ (av)b +cv2 Dv + e) calculating a proportionality coefficient of the horizontal reaction coefficient of the soil, wherein:
m is a proportional coefficient of a horizontal reaction force coefficient of the soil;
psi is a regional experience correction coefficient;
and v is the surface wave speed of the soil layer in the calculated depth range, and the unit is m/s.
The method for calculating the proportional coefficient of the soil horizontal reaction coefficient based on the surface wave exploration technology has the advantages of simple and quick calculation process, capability of calculating required parameters from transient surface wave exploration without drilling, horizontal load, cross plate shearing test and the like, strong field adaptability, capability of being used for different types of soil bodies such as clay, sandy soil, gravelly soil and the like, nondestructive exploration, small damage to the field, simple used equipment, low cost and high efficiency, still has higher calculation precision under the condition of not using experience to correct the coefficient, obvious advantages in areas with difficult drilling and in-situ test implementation in cities with dense buildings and the like, and great popularization and use values and wide application prospects.
A computer program for implementing the above-described method of calculating a proportionality coefficient of a soil level reaction force coefficient.
An information data processing terminal for realizing the method for calculating the proportionality coefficient of the soil level reaction force coefficient.
A computer-readable storage medium comprising instructions which, when executed on a computer, cause the computer to perform the above-described method of calculating a proportionality coefficient of soil level reaction force coefficients.
In the above embodiments, the implementation may be wholly or partially realized by software, hardware, firmware, or any combination thereof. When used in whole or in part, can be implemented in a computer program product that includes one or more computer instructions. When loaded or executed on a computer, cause the flow or functions according to embodiments of the invention to occur, in whole or in part. The computer may be a general purpose computer, a special purpose computer, a network of computers, or other programmable device. The computer instructions may be stored in a computer readable storage medium or transmitted from one computer readable storage medium to another, for example, the computer instructions may be transmitted from one website site, computer, server, or data center to another website site, computer, server, or data center via wire (e.g., coaxial cable, fiber optic, digital Subscriber Line (DSL), or wireless (e.g., infrared, wireless, microwave, etc.)). The computer readable storage medium can be any available medium that can be accessed by a computer or a data storage device, such as a server, a data center, etc., that includes one or more available media. The usable medium may be a magnetic medium (e.g., floppy Disk, hard Disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., solid State Disk (SSD)), among others.
The above description is only for the preferred embodiment of the present invention, and is not intended to limit the present invention in any way, and all simple modifications, equivalent changes and modifications made to the above embodiment according to the technical spirit of the present invention are within the scope of the technical solution of the present invention.