Mud geological building pile foundation positioning method

Abstract

The application provides a silt geology building pile foundation positioning method, and belongs to the technical field of pile foundation construction. The silt geology building pile foundation positioning method comprises the following steps: determining a main target detection area and a subordinate target area of silt geology; controlling a seismic source to send seismic waves from the subordinate target area to the main area; obtaining a first first-arrival time of the seismic waves from any geophone; determining a first target velocity of seismic wave propagation at each position in the main target detection area and a velocity inversion graph by using a seismic wave tomography method according to the first first-arrival time; determining the range of a silt geology high stress area in the velocity inversion graph according to the first target velocity and the corresponding relationship between velocity and stress; comparing the values of the vibration sensors in each detection hole in the high stress area range to determine the best piling position. The application solves the problems that the existing silt geology building pile foundation positioning method cannot quickly determine the range of the high stress area and cannot determine the best piling position.

Description

Technical Field

[0001] This invention belongs to the field of pile foundation construction technology, and more specifically, relates to a method for locating building pile foundations in silty geological conditions. Background Technology

[0002] With the significant increase in people's demands for living environments and the role of water conservancy and transportation in promoting industrial development, a large number of buildings have emerged in coastal, riverside, and lakeside areas. Due to the special geographical location, many buildings are located in silt or silty soil areas. However, silt is a soft soil with a natural water content greater than the fluidity limit and a porosity greater than 1.5. It is a modern sediment formed under conditions of slow water flow in bays, lakes, or river bends with the participation of microorganisms. It is rich in organic matter, often grayish-black in color, has low mechanical strength, and high compressibility. The presence of silt in the foundation can easily cause building subsidence. Therefore, when constructing buildings in silty soil areas, it is generally necessary to establish pile foundations in the silty soil to support the building. However, due to the special properties of silty soil, some silty soils have high stress and are suitable for pile driving, while others have low stress and are not suitable for pile drivers. The traditional method of pile foundation positioning is to measure the silty soil at each location to determine where it is suitable for pile driving. However, since the number of piles required for a building is very large, if the traditional pile foundation positioning method is used for each location, it is not possible to quickly determine the range of high stress areas and to quickly find the positioning points of the pile foundations, resulting in low work efficiency.

[0003] Currently, existing methods for locating pile foundations in silty geological structures often suffer from the inability to quickly determine the range of high-stress areas and the inability to determine the optimal pile driving location. Summary of the Invention

[0004] In view of this, the present invention provides a method for locating building pile foundations in silty geology, which aims to solve the problems of existing methods for locating building pile foundations in silty geology, which cannot quickly determine the range of high stress areas and cannot determine the optimal pile driving position.

[0005] This invention is implemented as follows:

[0006] This invention provides a method for locating building pile foundations in silty geological conditions, comprising the following steps:

[0007] S10: Determine the target area for pre-positioning the pile foundation and the area for placing the seismic source in the silty geology, wherein a first number of detection holes are set in the target area for pre-positioning the pile foundation, and a first number of geophones and a second number of vibration sensors are set in the first number of detection holes.

[0008] S20: Control the source to send seismic waves from the source placement area to the pre-positioning target area of ​​the pile foundation;

[0009] S30: The detector collects the time and wavelength of the first seismic wave received in each probe hole, i.e. the first arrival time and wavelength of the seismic wave, and uploads the collected time data to the host computer.

[0010] S40: The host computer determines the wave velocity and wavelength inversion map based on the first arrival time collected by the detector and the first arrival travel time tomography method of seismic waves.

[0011] S50: Based on the wave velocity, wavelength, and the relationship between wave velocity and stress, as well as the relationship between wavelength and stress, determine the extent of high-stress zones in silty geology in the wave velocity and wavelength inversion diagram.

[0012] S60: Compare the vibration sensor readings in each probe hole within the high-stress zone to determine the optimal piling location;

[0013] Among them, the pre-positioning target area of ​​the pile foundation is the silty geological area where pile foundation positioning detection is required; the source placement area is the location where the source is placed; the wave velocity and wavelength inversion map is used to indicate the first target seismic wave velocity and the first target seismic wave wavelength at each location in the pre-positioning target area of ​​the pile foundation.

[0014] Based on the above technical solution, the method for locating building pile foundations in silty geology according to the present invention can be further improved as follows:

[0015] In step S10, the pre-positioning target area of ​​the pile foundation is divided into a third number of grid cells.

[0016] Grid cells can be squares, rectangles, triangles, etc.

[0017] The location of the detection holes is selected using a polar coordinate system, and the number of detection holes is N, wherein:

[0018] The coordinates of the probe hole are P. i (r i θ i ), ;

[0019] The value of i is: 1, 2, 3, ..., N

[0020] Where i is the sequence number of the probe hole, r i Let θ be the polar diameter of the i-th probe aperture. i Let be the polar angle of the i-th probe.

[0021] In step S20, the earthquake source placement area includes four areas: the first earthquake source placement area, the second earthquake source placement area, the third earthquake source placement area, and the fourth earthquake source placement area.

[0022] Furthermore, in step S20, the transmission of seismic waves from the seismic source placement area to the pre-positioning target area of ​​the pile foundation specifically includes the following steps:

[0023] Step 1: The seismic source sends seismic waves from the first seismic source placement area to the pre-positioned target area of ​​the pile foundation;

[0024] Step 2: After time T1, the seismic source sends seismic waves from the second seismic source placement area to the pre-positioned target area of ​​the pile foundation;

[0025] Step 3: After T2, the seismic source sends seismic waves from the third seismic source placement area to the pre-positioned target area of ​​the pile foundation;

[0026] Step 4: After time T3, the seismic source sends seismic waves from the fourth seismic source placement area to the pre-positioned target area of ​​the pile foundation;

[0027] Wherein, T1 is the time interval from when the seismic source sends seismic waves from the first seismic source placement area to the pre-positioning target area of ​​the pile foundation until the last vibration sensor receives the data of vibration; T2 is the time interval from when the seismic source sends seismic waves from the second seismic source placement area to the pre-positioning target area of ​​the pile foundation until the last vibration sensor receives the data of vibration; T3 is the time interval from when the seismic source sends seismic waves from the third seismic source placement area to the pre-positioning target area of ​​the pile foundation until the last vibration sensor receives the data of vibration.

[0028] In step S40, the wave velocity-wavelength inversion map is determined based on the seismic wave first arrival travel time tomography method, specifically including the following steps:

[0029] S401: Establish a seismic wave analysis model, wherein each grid cell in the pre-positioning target area of ​​the pile foundation has a corresponding seismic wave velocity and seismic wave wavelength;

[0030] S402: Perform multiple rounds of iterative inversion calculations on the seismic wave analysis model based on the seismic wave velocity, the seismic wave first arrival time, the seismic wave wavelength, and the seismic wave first arrival wavelength, so as to update the initial seismic wave velocity in the seismic wave analysis model to the first target seismic wave velocity and the initial seismic wave wavelength in the seismic wave analysis model to the first target seismic wave wavelength.

[0031] S403: Based on the wave velocity and wavelength of the first target seismic wave, obtain the wave velocity and wavelength inversion map of the pre-positioning target area of ​​the pile foundation.

[0032] In step S50, based on the wave velocity of the first target seismic wave and the correspondence between velocity and stress, the range of the high-stress zone in the silty geology is determined in the wave velocity-wavelength inversion map. This specifically includes the following steps:

[0033] Step 1: Obtain the reference velocity and reference wavelength, wherein the reference velocity is the velocity of the seismic wave when it is destroyed near the silty geology, and the reference wavelength is the wavelength of the seismic wave when it is destroyed near the silty geology;

[0034] Step 2: Based on the difference between the reference velocity and the second target seismic wave velocity of any grid cell and the difference between the reference wavelength and the second target seismic wave wavelength of any grid cell, determine the first wave velocity difference value and the first wavelength difference value of any grid cell, wherein the second target seismic wave velocity indicates the wave velocity of the seismic wave in any grid cell when the silty geology is not mined, and the second target seismic wave wavelength indicates the wavelength of the seismic wave in any grid cell when the silty geology is not mined;

[0035] Step 3: Determine the second velocity difference value of any grid cell based on the difference between the first target seismic wave velocity and the second target seismic wave velocity; determine the second wavelength difference value of any grid cell based on the difference between the first target seismic wave wavelength and the second target seismic wave wavelength.

[0036] Step 4: Determine the wave velocity anomaly coefficient and wavelength anomaly coefficient of any grid cell based on the ratio of the second wave velocity difference to the first wave velocity difference and the ratio of the second wavelength difference to the first wavelength difference.

[0037] Step 5: In response to the wave velocity anomaly coefficient of any grid cell being greater than the second set threshold, and in response to the length anomaly coefficient of any grid cell being greater than the second set threshold, determine that any grid cell is a high-stress grid cell;

[0038] Step Six: Based on each of the high-stress grid cells, determine the extent of the high-stress zone in the silt geology in the wave velocity-wavelength inversion diagram.

[0039] Furthermore, in step S402, based on the initial wave velocity, the first arrival time of the seismic wave, the initial wavelength, and the first arrival wavelength of the seismic wave, multiple rounds of iterative inversion calculations are performed on the seismic wave analysis model to update the initial wave velocity in the seismic wave analysis model to the first target seismic wave velocity and the initial wavelength to the first target seismic wave wavelength. Specifically, this includes the following steps:

[0040] Step 1: For any round of iterative inversion calculation, perform ray tracing on each ray path according to the initial wave velocity and initial wavelength used in this round, and determine the ray path length of each ray path in this round, each grid cell passed through by each ray path in this round, and the reference duration and reference wavelength required for the seismic wave to propagate along each ray path in this round.

[0041] Step 2: Based on the first arrival time of the seismic wave and the reference duration of this round, as well as the first arrival wavelength of the seismic wave and the reference wavelength of this round, determine the travel time residual and the loss function of this round;

[0042] If the value of the loss function in this round is greater than the first preset threshold, the following steps are executed:

[0043] Step 3: Based on the travel time residual and the ray path length mentioned in this round, determine the slowness update amount of any grid cell in this round, wherein the slowness update amount mentioned in this round is used to indicate the initial slowness update amount in this round; the initial slowness in this round is the reciprocal of the initial wave velocity in this round.

[0044] Step 4: Update the initial wave velocity to the reciprocal of the sum of the slowness update amount and the initial slowness, and use the updated initial wave velocity as the initial wave velocity for the next round;

[0045] Step 5: If the value of the loss function in this round is not greater than the first set threshold, stop iterative inversion of the seismic wave analysis model, and use the updated initial wave velocity as the target seismic wave velocity.

[0046] The seismic wave analysis model can be either a constant wave velocity wavelength model or a variable wave velocity wavelength model.

[0047] The constant wave velocity wavelength model is:

[0048] P1=C1(a,b)+λ1(a,b)

[0049] C1(a,b)=k c1 λ1(a,b)=k λ1

[0050] The pre-positioning target area of ​​the pile foundation is divided into A rows and B columns of grid cells, where a = 1, 2, ..., A, b = 1, 2, ..., B, P1 is the constant wave velocity wavelength model function, C1(a,b) is the constant wave velocity function corresponding to the grid cell in the a-th row and b-th column of the pre-positioning target area of ​​the pile foundation, λ1(a,b) is the constant wavelength function corresponding to the grid cell in the a-th row and b-th column of the pre-positioning target area of ​​the pile foundation, and k c1 Let k be the constant wave velocity corresponding to the grid cell in the a-th row and b-th column of the pre-positioning target area of ​​the pile foundation. λ1 The fixed wavelength k corresponds to the grid cell in the a-th row and b-th column of the pre-positioning target area of ​​the pile foundation. c1 With k λ1 All are constants.

[0051] The variable wave speed wavelength model is as follows:

[0052] P2=C2(a,b)+λ2(a,b)

[0053] C2(a,b)=c0+d(b,b), λ2(a,b)=λ0+e(a,b)

[0054] In this design, the pre-positioning target area of ​​the pile foundation is divided into grid cells of row A and column B, where a is the horizontal coordinate of the grid cell, b is the vertical coordinate of the grid cell, a = 1, 2, ..., A, b = 1, 2, ..., B, P2 is the variable wave velocity wavelength model function, C2(a,b) is the variable wave velocity function corresponding to the grid cell in row a and column b of the pre-positioning target area of ​​the pile foundation, λ2(a,b) is the variable wavelength function corresponding to the grid cell in row a and column b of the pre-positioning target area of ​​the pile foundation, d(a,b) is the variable wave velocity corresponding to the grid cell in row a and column b of the pre-positioning target area of ​​the pile foundation, e(a,b) is the variable wave velocity corresponding to the grid cell in row a and column b of the pre-positioning target area of ​​the pile foundation, c0 is the initial value of the variable wave velocity function, λ0 is the initial value of the variable wavelength function, and c0 and λ0 are both constants.

[0055] In step S60, the vibration sensor values ​​in each detection hole within the high-stress zone are compared to determine the optimal piling location. This specifically includes the following steps:

[0056] Step 1: The vibration sensor collects vibration data from each detection hole within the high-stress zone and uploads the data to the host computer;

[0057] Step 2: The host computer compares and analyzes the vibration values ​​in each detection hole to determine the optimal piling location;

[0058] The optimal piling position is when the vibration value inside the detection hole is less than the third set threshold.

[0059] The seismic source is an electromagnetic shock source, which consists of an energy storage main control box and an electromagnetic shock hammer.

[0060] Compared with existing technologies, the beneficial effects of the method for locating building pile foundations in silty geology provided by this invention are as follows: It determines the main target detection area and subordinate target areas in the silty geology; it controls the seismic source to send seismic waves from the subordinate target area to the main area; it obtains the first arrival time of the seismic waves from any detector; based on the first arrival time, it uses seismic wave tomography to determine the first target velocity and velocity inversion map of the seismic wave propagation at each location within the main target detection area; based on the first target velocity and the correspondence between velocity and stress, it determines the range of the high-stress zone in the silty geology in the velocity inversion map, improving work efficiency; and it compares the values ​​of vibration sensors in each detection hole within the high-stress zone to determine the optimal pile driving position. Attached Figure Description

[0061] To more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the description of the embodiments of the present invention will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0062] Figure 1 A flowchart for the method of locating pile foundations in silty geological structures;

[0063] Figure 2 A flowchart of a tomographic imaging method for locating pile foundations in silty geological conditions; Detailed Implementation

[0064] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0065] Therefore, the following detailed description of the embodiments of the invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely to illustrate selected embodiments of the invention. All other embodiments obtained by those skilled in the art based on the embodiments of the invention without inventive effort are within the scope of protection of the invention.

[0066] It should be noted that similar labels and letters in the following figures indicate similar items. Therefore, once an item is defined in one figure, it does not need to be further defined and explained in subsequent figures.

[0067] In the description of this invention, it should be understood that the terms "center," "longitudinal," "lateral," "length," "width," "thickness," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," "clockwise," and "counterclockwise," etc., indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are only for the convenience of describing this invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on this invention.

[0068] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of this invention, "a plurality of" means two or more, unless otherwise explicitly specified.

[0069] like Figure 1 The image shows a first embodiment of a method for locating building pile foundations in silty geology provided by the present invention. This embodiment includes the following steps:

[0070] S10: Determine the target area for pre-positioning of pile foundations and the area for placement of seismic sources in silty geology, wherein a first number of detection holes are set in the target area for pre-positioning of pile foundations, and a first number of geophones and a second number of vibration sensors are set in the first number of detection holes.

[0071] S20: Control the source to send seismic waves from the source placement area to the pre-positioning target area of ​​the pile foundation;

[0072] S30: The detector collects the time and wavelength of the first seismic wave received in each probe hole, i.e. the first arrival time and wavelength of the seismic wave, and uploads the collected time data to the host computer.

[0073] S40: The host computer determines the wave velocity and wavelength inversion map based on the first arrival time collected by the detector and the first arrival travel time tomography method of seismic waves.

[0074] S50: Based on the wave velocity, wavelength, and the relationship between wave velocity and stress, as well as the relationship between wavelength and stress, determine the extent of high-stress zones in silty geology in the wave velocity and wavelength inversion diagram.

[0075] S60: Compare the vibration sensor readings in each probe hole within the high-stress zone to determine the optimal piling location;

[0076] Among them, the pre-positioning target area of ​​the pile foundation is the silty geological area where pile foundation positioning detection is required; the source placement area is the location where the source is placed; the wave velocity and wavelength inversion map is used to indicate the first target seismic wave velocity and the first target seismic wave wavelength at each location in the pre-positioning target area of ​​the pile foundation.

[0077] In the above technical solution, in step S10, the target area for pre-positioning the pile foundation is divided into a third number of grid cells.

[0078] In the above technical solution, the location of the detection holes is selected using a polar coordinate system, and the number of detection holes is N, where:

[0079] The coordinates of the probe hole are P i (r i θ i ), ;

[0080] The value of i can be: 1, 2, 3, ..., N

[0081] Where i is the sequence number of the probe hole, r i Let θ be the polar diameter of the i-th probe aperture. i Let be the polar angle of the i-th probe.

[0082] In the above technical solution, step S20 includes four earthquake source placement areas, namely the first earthquake source placement area, the second earthquake source placement area, the third earthquake source placement area, and the fourth earthquake source placement area.

[0083] Furthermore, in the above technical solution, step S20, which involves the source sending seismic waves from the source placement area to the pre-positioning target area of ​​the pile foundation, specifically includes the following steps:

[0084] Step 1: The seismic source sends seismic waves from the first seismic source placement area to the pre-positioned target area of ​​the pile foundation;

[0085] Step 2: After time T1, the seismic source sends seismic waves from the second seismic source placement area to the pre-positioned target area of ​​the pile foundation;

[0086] Step 3: After T2, the seismic source sends seismic waves from the third seismic source placement area to the pre-positioned target area of ​​the pile foundation;

[0087] Step 4: After time T3, the seismic source sends seismic waves from the fourth seismic source placement area to the pre-positioned target area of ​​the pile foundation;

[0088] Wherein, T1 is the time interval from when the seismic source sends seismic waves from the first seismic source placement area to the pre-positioning target area of ​​the pile foundation until the last vibration sensor receives the data of vibration; T2 is the time interval from when the seismic source sends seismic waves from the second seismic source placement area to the pre-positioning target area of ​​the pile foundation until the last vibration sensor receives the data of vibration; T3 is the time interval from when the seismic source sends seismic waves from the third seismic source placement area to the pre-positioning target area of ​​the pile foundation until the last vibration sensor receives the data of vibration.

[0089] Among them, such as Figure 2 As shown, in the above technical solution, step S40, determining the wave velocity-wavelength inversion map based on the seismic wave first arrival travel time tomography method, specifically includes the following steps:

[0090] S401: Establish a seismic wave analysis model, in which each grid cell in the pre-positioning target area of ​​the pile foundation has a corresponding seismic wave velocity and seismic wave wavelength;

[0091] S402: Perform multiple rounds of iterative inversion calculations on the seismic wave analysis model based on the seismic wave velocity, the first arrival time of the seismic wave, and the wavelength of the seismic wave and the first arrival wavelength of the seismic wave, so as to update the initial seismic wave velocity in the seismic wave analysis model to the first target seismic wave velocity, and update the initial seismic wave wavelength in the seismic wave analysis model to the first target seismic wave wavelength.

[0092] S403: Based on the wave velocity and wavelength of the first target seismic wave, obtain the wave velocity and wavelength inversion map of the pre-positioning target area of ​​the pile foundation.

[0093] In the above technical solution, step S50, based on the wave velocity of the first target seismic wave and the correspondence between velocity and stress, determines the range of the high-stress zone in the silty geology in the wave velocity-wavelength inversion map, specifically including the following steps:

[0094] Step 1: Obtain the reference velocity and reference wavelength, where the reference velocity is the velocity of the seismic wave when it is destroyed near the silty geology, and the reference wavelength is the wavelength of the seismic wave when it is destroyed near the silty geology.

[0095] Step 2: Based on the difference between the reference velocity and the second target seismic wave velocity of any grid cell and the difference between the reference wavelength and the second target seismic wave wavelength of any grid cell, determine the first wave velocity difference and the first wavelength difference of any grid cell. The second target seismic wave velocity indicates the wave velocity of the seismic wave in any grid cell when the silty geology is not mined, and the second target seismic wave wavelength indicates the wavelength of the seismic wave in any grid cell when the silty geology is not mined.

[0096] Step 3: Based on the difference between the wave velocity of the first target seismic wave and the wave velocity of the second target seismic wave, determine the second wave velocity difference value of any grid cell; based on the difference between the wavelength of the first target seismic wave and the wavelength of the second target seismic wave, determine the second wavelength difference value of any grid cell.

[0097] Step 4: Determine the wave velocity anomaly coefficient and wavelength anomaly coefficient of any grid cell based on the ratio of the second wave velocity difference to the first wave velocity difference and the ratio of the second wavelength difference to the first wavelength difference.

[0098] Step 5: If the wave velocity anomaly coefficient of any grid cell is greater than the second set threshold, and if the length anomaly coefficient of any grid cell is greater than the second set threshold, then any grid cell is determined to be a high-stress grid cell.

[0099] Step 6: Determine the extent of the high-stress zone in the silt geology based on each high-stress grid cell in the wave velocity-wavelength inversion diagram.

[0100] Furthermore, in the above technical solution, in step S402, based on the initial wave velocity, the first arrival time of the seismic wave, the initial wavelength, and the first arrival wavelength of the seismic wave, multiple rounds of iterative inversion calculations are performed on the seismic wave analysis model to update the initial wave velocity in the seismic wave analysis model to the first target seismic wave velocity and the initial wavelength to the first target seismic wave wavelength. Specifically, this includes the following steps:

[0101] Step 1: For any round of iterative inversion calculation, perform ray tracing on each ray path according to the initial wave velocity and initial wavelength used in this round, and determine the ray path length of each ray path in this round, each grid cell passed through by each ray path in this round, and the reference duration and reference wavelength required for the seismic waves to propagate along each ray path in this round.

[0102] Step 2: Based on the first arrival time of the seismic wave and the reference duration of this cycle, as well as the first arrival wavelength of the seismic wave and the reference wavelength of this cycle, determine the travel time residual and the loss function of this cycle;

[0103] If the value of the loss function in this round is greater than the first set threshold, perform the following steps:

[0104] Step 3: Based on the current travel time residual and ray path length, determine the slowness update amount for any grid cell in the current round. The current slowness update amount is used to indicate the update amount of the initial slowness in the current round; the initial slowness in the current round is the reciprocal of the initial wave velocity in the current round.

[0105] Step 4: Update the initial wave velocity to the reciprocal of the sum of the slowness update amount and the initial slowness, and use the updated initial wave velocity as the initial wave velocity for the next round.

[0106] Step 5: If the value of the loss function in this round is not greater than the first set threshold, stop iterative inversion of the seismic wave analysis model, and use the updated initial wave velocity as the target seismic wave velocity.

[0107] The slow update function is:

[0108]

[0109] In the formula, ΔS(a,b) is the slow update function, and ΔT j Let L be the travel time residual corresponding to the j-th ray path. j Let Q be the length of the first ray path of the j-th ray path, Q be the number of rays passing through the grid cell at position (a,b), and j be the ray number.

[0110] The loss function is:

[0111]

[0112] In the formula, F is the loss function, and T 0jLet ΔT be the initial arrival time of the seismic wave from any seismic source area to any geophone. j Let T be the travel time residual corresponding to the j-th ray path. j T is the first reference time for the seismic wave to travel from any source location area to any detector. j The range of values ​​for is T 0j ~T j , where j = 1, 2, ..., Z, and Z refers to the number of ray paths.

[0113] The first set threshold can be preset, for example, the first set threshold is 1.0e-2 (or 0.01), 1.0e-3 (or 0.001), etc., and this disclosure does not limit it; the second set threshold can be preset, for example, the second set threshold is 0.7, 0.75, etc., and this disclosure does not limit it; the third set threshold can be preset, for example, the third set threshold is 1.0e-2, 1.0e-3, etc.

[0114] In the above technical solution, step S60, comparing the vibration sensor values ​​in each detection hole within the high-stress zone to determine the optimal piling position, specifically includes the following steps:

[0115] Step 1: The vibration sensor collects vibration data from each detection hole within the high-stress zone and uploads the data to the host computer;

[0116] Step 2: The host computer compares and analyzes the vibration values ​​in each detection hole to determine the optimal piling location;

[0117] The optimal piling position is when the vibration value inside the detection hole is less than the third set threshold.

[0118] In the above technical solution, the seismic source is an electromagnetic shock source, which consists of an energy storage main control box and an electromagnetic shock hammer.

[0119] The electromagnetic shock source used was manufactured by Jining Huakuang Machinery Equipment Co., Ltd.

[0120] The energy storage main control box is model EHSS5000-EC, manufactured by Jining Huakuang Machinery Equipment Co., Ltd., and the electromagnetic impact hammer is model EHSS5000-S, manufactured by Jining Huakuang Machinery Equipment Co., Ltd.

[0121] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for positioning a sludge geo-architectural pile foundation, characterized in that, The method comprises the following steps: S10: determining a pile foundation pre-positioning target area and a seismic source placement area of the silt geology, wherein a first number of detection holes are arranged in the pile foundation pre-positioning target area, and a first number of geophones and a second number of vibration sensors are arranged in the first number of detection holes; S20: controlling the seismic source to send seismic waves from the seismic source placement area to the pile foundation pre-positioning target area; S30: collecting, by the geophones, a time at which a first received seismic wave in each detection hole is received and a wavelength of the first received seismic wave, i.e., a seismic wave first arrival time and a seismic wave first arrival wavelength, and uploading the collected time data to an upper computer; S40: determining, by the upper computer, a wave velocity and wavelength inversion map according to the first arrival time collected by the geophones and according to a seismic wave first arrival travel time tomography method; S50: determining a range of a high stress area in the silt geology in the wave velocity and wavelength inversion map according to a wave velocity of the seismic wave, a wavelength of the seismic wave, a corresponding relationship between the wave velocity and stress, and a corresponding relationship between the wavelength and stress; S60: comparing values of the vibration sensors in each detection hole in the high stress area range to determine an optimal pile driving position. The pile foundation pre-positioning target area is a silt geology area that needs to be detected for pile foundation positioning. The seismic source placement area is a position where the seismic source is placed. The wave velocity and wavelength inversion map is used to indicate a first target seismic wave velocity and a first target seismic wave wavelength of seismic wave propagation at each position in the pile foundation pre-positioning target area.

2. A method of positioning a sludge geo-architectural pile foundation according to claim 1, characterized in that, In step S10, the pile foundation pre-positioning target area is divided into a third number of grid units.

3. The method for positioning a sludge geo-architectural pile foundation according to claim 1, characterized in that, The detection hole positions are selected using a polar coordinate system, and the number of detection holes is N, wherein: The coordinates of the probe hole are P i (r i , θ i ) The value of i is 1, 2, 3, …, N. where i is the serial number of the detection hole, r i is the polar radius of the ith detection hole, θ i is the polar angle of the ith detection hole.

4. The method for positioning a sludge geo-architectural pile foundation according to claim 1, characterized in that, In step S20, the seismic source placement area includes four, namely a first seismic source placement area, a second seismic source placement area, a third seismic source placement area, and a fourth seismic source placement area.

5. A method of positioning a sludge geo-architectural pile foundation according to claim 4, characterized in that, In step S20, the seismic source sends seismic waves from the seismic source placement area to the pile foundation pre-positioning target area, which specifically comprises the following steps: Step one: the seismic source sends seismic waves from the first seismic source placement area to the pile foundation pre-positioning target area; Step two: after T1 time, the seismic source sends seismic waves from the second seismic source placement area to the pile foundation pre-positioning target area; Step three: after T2 time, the seismic source sends seismic waves from the third seismic source placement area to the pile foundation pre-positioning target area; Step four: after T3 time, the seismic source sends seismic waves from the fourth seismic source placement area to the pile foundation pre-positioning target area; Wherein T1 is the time interval from the seismic source sending seismic waves from the first seismic source placement area to the pile foundation pre-positioning target area to the last vibration sensor receiving vibration data; T2 is the time interval from the seismic source sending seismic waves from the second seismic source placement area to the pile foundation pre-positioning target area to the last vibration sensor receiving vibration data; T3 is the time interval from the seismic source sending seismic waves from the third seismic source placement area to the pile foundation pre-positioning target area to the last vibration sensor receiving vibration data.

6. The method for positioning a sludge geo-architectural pile foundation according to claim 1, characterized in that, In step S40, the wave velocity and wavelength inversion map is determined according to the seismic wave first arrival travel time tomography method, which specifically comprises the following steps: S401: Establish a seismic wave analysis model, wherein each grid cell in the pile foundation pre-positioning target area has a corresponding seismic wave velocity and seismic wave wavelength; S402: Perform multi-round iterative inversion calculation on the seismic wave analysis model according to the seismic wave velocity and the seismic wave first arrival time and the seismic wave wavelength and the seismic wave first arrival wavelength to update the initial seismic wave velocity in the seismic wave analysis model to the first target seismic wave velocity and update the initial seismic wave wavelength in the seismic wave analysis model to the first target seismic wave wavelength; S403: According to the first target seismic wave velocity and the first target seismic wave wavelength, obtain the velocity-wavelength inversion graph of the pile foundation pre-positioning target area.

7. The method for positioning a sludge geo-architectural pile foundation according to claim 1, characterized in that, In step S50, according to the first target seismic wave velocity and the corresponding relationship between velocity and stress, the range of the high stress area in the silt geology is determined in the velocity-wavelength inversion graph, which specifically includes the following steps: Step one: Obtain a reference velocity and a reference wavelength, wherein the reference velocity is the velocity of the seismic wave when the silt geology is adjacent to damage, and the reference wavelength is the wavelength of the seismic wave when the silt geology is adjacent to damage; Step two: According to the difference between the reference velocity and the second target seismic wave velocity of any grid cell and the difference between the reference wavelength and the second target seismic wave wavelength of any grid cell, determine the first velocity difference value and the first wavelength difference value of any grid cell, wherein the second target seismic wave velocity indicates the velocity of the seismic wave in any grid cell under the condition that the silt geology is not exploited, and the second target seismic wave wavelength indicates the wavelength of the seismic wave in any grid cell under the condition that the silt geology is not exploited; Step three: According to the difference between the first target seismic wave velocity and the second target seismic wave velocity, determine the second velocity difference value of any grid cell, and according to the difference between the first target seismic wave wavelength and the second target seismic wave wavelength, determine the second wavelength difference value of any grid cell; Step four: According to the ratio of the second velocity difference value to the first velocity difference value and the ratio of the second wavelength difference value to the first wavelength difference value, determine the velocity anomaly coefficient and the wavelength anomaly coefficient of any grid cell; Step five: In response to the velocity anomaly coefficient of any grid cell being greater than a second set threshold, and in response to the wavelength anomaly coefficient of any grid cell being greater than a second set threshold, determine that any grid cell is a high stress grid cell; Step six: According to each high stress grid cell, determine the range of the high stress area of the silt geology in the velocity-wavelength inversion graph.

8. A method of positioning a sludge geo-architectural pile foundation according to claim 6, characterized in that, In step S402, according to the initial velocity and the seismic wave first arrival time and the initial wavelength and the seismic wave first arrival wavelength, perform multi-round iterative inversion calculation on the seismic wave analysis model to update the initial velocity in the seismic wave analysis model to the first target seismic wave velocity and update the initial wavelength in the seismic wave analysis model to the first target seismic wave wavelength, which specifically includes the following steps: Step one: for any round of iterative inversion calculation, according to the initial wave velocity and initial wavelength adopted in this round, ray tracing is performed on each ray path to determine the ray path length of each ray path in this round, each grid cell passed by each ray path in this round, and the reference time length and reference wavelength required for seismic wave propagation along each ray path in this round; Step two: according to the seismic wave first arrival time and the reference time length and reference wavelength in this round, the travel time residual in this round and the loss function in this round are determined; If the value of the loss function in this round is greater than a first set threshold, the following steps are performed: Step three: according to the travel time residual in this round and the ray path length, the slowness update amount of any grid cell in this round is determined, wherein the slowness update amount in this round is used to indicate the update amount of the initial slowness in this round; the initial slowness in this round is the inverse of the initial wave velocity in this round; Step four: the initial wave velocity is updated to the inverse of the sum of the slowness update amount and the initial slowness, and the updated initial wave velocity is used as the initial wave velocity adopted in the next round; Step five: if the value of the loss function in this round is not greater than the first set threshold, the iterative inversion of the seismic wave analysis model is stopped, and the updated initial wave velocity is used as the target seismic wave velocity.

9. The method for positioning a sludge geo-architectural pile foundation according to claim 1, characterized in that, In step S60, the values of the vibration sensors in each detection hole in the high stress area range are compared to determine the best piling position, which specifically includes the following steps: Step one: the vibration sensors collect vibration data in each detection hole in the high stress area range, and upload the data to the upper computer; Step two: the upper computer compares and analyzes the vibration values in each detection hole to obtain the best piling position; Wherein, the vibration value in the detection hole is less than the third set threshold, which is the best piling position.

10. The method for positioning a sludge geo-architectural pile foundation according to claim 1, characterized in that, The seismic source is an electromagnetic impact seismic source, which is composed of an energy storage master control box and an electromagnetic impact hammer.

Images