Interpretation method and system for inner core pile-outer core pile non-coaxial angle of stiff composite pile

Abstract

The invention discloses an interpretation method and system for an inner core pile-outer core pile non-coaxial angle of a stiff composite pile, and belongs to the field of pile foundation detection. Comprising the steps that a site is leveled, a plurality of adjacent composite piles are selected, and the pile distance between two test piles, the thickness and the pile length of an inner core pile and an outer core pile and the like are measured through a measuring tape; an acoustic transmitter and an acoustic receiver of an ultrasonic instrument are placed in the inner core piles of the two composite piles respectively, monitoring points of acoustic transmitting and receiving signals are designed according to a sector scanning method, and the acoustic transmitting points and the acoustic receiving points are arranged at equal intervals and depths; according to the provided SIRT and a high-precision grid expansion algorithm, carrying out inversion calculation on a grid wave velocity cloud picture; and the non-coaxial angle of the inner core pile and the outer core pile is calculated according to a defined adaptation degree formula and an interpretation algorithm. The interpretation algorithm of the non-coaxial angle of the inner core pile and the outer core pile of the stiff composite pile is provided for the first time, intelligent detection of the pile forming quality of the stiff composite pile is achieved, and the method has wide application prospects in quality control of the stiff composite pile.

Description

Technical Field

[0001] This invention discloses a method and system for interpreting the non-coaxial angle between the inner core pile and the outer core pile of a rigid composite pile, which belongs to the field of pile foundation testing. Background Technology

[0002] In my country's coastal areas, the widespread distribution of soft soil is characterized by high water content, high compressibility, and low strength, making infrastructure built on such foundations prone to uneven settlement. To address this issue, reinforced composite piles, as an effective foundation reinforcement technology, leverage the high strength of pipe piles and the larger cross-section and high side friction of cement-soil piles to significantly control foundation settlement. Therefore, they are widely used in infrastructure construction in soft soil areas.

[0003] In recent years, with the large-scale application of reinforced composite piles in coastal areas such as Jiangsu, Zhejiang, Shanghai, and Guangdong, their excellent settlement control effect and high cost-effectiveness have gained widespread recognition in the engineering community. Simultaneously, with the continuous improvement of composite pile forming technology, construction processes, and related specifications, the application scope of reinforced composite piles has further expanded. However, in actual construction, due to the quality problems of cement-soil mixing pile forming and the difficulty in achieving perfect verticality when driving pipe piles into cement-soil piles, varying degrees of non-coaxiality exist between the inner and outer pile cores. This brings serious uncertainties to the pile foundation's bearing capacity and settlement control. In large-scale applications, this non-coaxiality is currently a significant cause of engineering disasters.

[0004] The non-coaxial angle refers to the angle between the axes of the inner core pile (pipe pile) and the outer core pile (cement-soil pile) and the normal axis of the ground plane. When the axes of the inner and outer piles are not completely coincident, it will affect the overall performance of the pile foundation and the settlement control effect. However, there is currently a lack of an effective detection method to quantify and measure this non-coaxial angle. Summary of the Invention

[0005] The technical problem to be solved by this invention is as follows: This invention proposes a method and system for interpreting the non-coaxial angle between the inner core pile and the outer core pile of a rigid composite pile. It adopts acoustic emission technology, array-type emission and acquisition deployment, and imaging interpretation algorithm for non-coaxial angle, which can accurately detect the non-coaxial characteristics of the inner and outer core piles in a rigid composite pile, so as to realize the test and evaluation of the non-coaxial angle between the inner core pile and the outer core pile of the composite pile.

[0006] To solve the above technical problems, the present invention adopts the following technical solution:

[0007] First, this invention proposes a method for interpreting the non-coaxial angle between the inner core pile and the outer core pile of a reinforced composite pile, the steps of which are as follows:

[0008] Step 1: Select any two adjacent composite piles A and B, and fill the inner core pile of the composite pile with water;

[0009] Step 2: Measure the thickness of the inner core pile and the outer core pile of the composite pile. and and pile spacing ;

[0010] Step 3: Place the acoustic transmitter and receiver in the inner core piles of pile A and pile B respectively, starting from the preset underground depth and arranging the points downwards at equal intervals.

[0011] Step 4: Following the sector scanning method, raise the acoustic emitter sequentially, transmit signals to all acoustic receivers at each depth point, and record the arrival time data of the sound waves arranged in column vectors according to the depth from low to high.

[0012] Step 5: Based on the arrival time of the sound waves, divide the plane between the transmitter and receiver into multiple grids, construct a path matrix, and calculate the grid wave velocity using an inversion algorithm;

[0013] Step 6: Calculate the high-precision grid wave velocity matrix using a high-precision grid wave velocity inversion algorithm based on the grid wave velocity.

[0014] Step 7: Based on the high-precision grid wave velocity matrix, calculate the non-coaxial angle between the inner core pile and the outer core pile using the fit formula and interpretation algorithm.

[0015] Furthermore, in step 1 of the present invention, the inner core pile is filled with water as a coupling agent for acoustic emission and reception, and there is no obvious silt suspended in the water.

[0016] Furthermore, the measurement process in step 2 of this invention involves performing multiple measurements and taking the average value.

[0017] Furthermore, in the method proposed in this invention, step 5 uses an inversion algorithm and a joint algebraic reconstruction method to iteratively solve for the grid wave velocity. The iteration termination condition is that the number of iterations reaches a preset value or the standard deviation of the sound wave arrival time is less than a threshold. The specific steps are as follows: the plane between the transmitter and the receiver is divided into... There are grids, assuming there are... The arrival time of the experimental sound waves and the wave velocity in the grid are calculated using the following set of equations:

[0018] (1)

[0019] In the formula, For the ray in the first OK Path length in column grid For the first Line number The reciprocal of the wave velocity of the grid; For the first The arrival time of the test sound wave for the ray is n less than m.

[0020] Furthermore, in step 6 of the method proposed in this invention, the high-precision grid wave velocity inversion algorithm includes resetting the positions of the acoustic emission point and the acoustic reception point with a preset step size, calculating the arrival time of the expanded acoustic wave based on the grid wave velocity matrix, and then calculating the high-precision grid wave velocity matrix through the inversion algorithm. Specifically, the line segment of each ray in each grid is recorded as a vector, and code is written to traverse all vectors in the grid sequentially according to the ray order, where the magnitude of the vector is the path of the ray in the grid. .

[0021] Furthermore, a high-precision grid wave velocity matrix is ​​calculated using a high-precision grid wave velocity inversion algorithm, as detailed below:

[0022] S601. Calculate the initial value of the reciprocal of the wave velocity of the grid. Assuming the wave velocity of each ray is uniform, the arrival time of the sound wave of each ray is distributed among the grids according to the path length of each grid cell. The reciprocal of the initial wave velocity is obtained by dividing the total time of all rays passing through the grid cells by the total path length within the grid. The calculation formula is as follows:

[0023] (2)

[0024] S602, based on the initial value of the reciprocal of the calculated wave velocity. and the arrival time of the test sound waves The reciprocal of the wave speed is solved by calculating the system of equations in formula (1) using the following iterative algorithm. The iterative equation is as follows:

[0025] (3)

[0026] In the formula, The relaxation factor is selected from 1.0 to 1.2. For the first Line number List a grid in The reciprocal of the grid wave velocity after the next iteration; For the first Line number List a grid in The reciprocal of the grid wave velocity after the next iteration; For the first The arrival time of the test sound wave of the ray; For the first The first ray OK Path length in column grid; For the first Line number The number of rays passing through the column grid; the iteration termination condition of formula (3) is the number of iterations. or the standard deviation of the arrival time of sound waves during iteration. .

[0027] S603, the reciprocal of the grid wave velocity in formula (3) Substituting into the following formula yields the low-precision grid wave velocity. sum matrix :

[0028] (4)

[0029] S604. Based on the obtained grid wave velocity matrix With a step size of 5cm or 10cm, the positions of the sound emission point and the sound reception point are reset from a depth of 0m downwards. The arrival time of the amplified sound wave is calculated according to formula (5):

[0030] (5)

[0031] S605, based on the calculated Repeat steps 5 and 6 to calculate the high-precision grid wave velocity matrix.

[0032] .

[0033] Furthermore, the interpretation method for non-coaxial angles is as follows:

[0034] S701. Based on the grid wave velocity cloud map, it is determined that the inner core pile and outer core pile are not coaxial. When the inner core pile and outer core pile are not coaxial, there will be an overlap with the calculated high-precision grid wave velocity matrix. Within the overlapping area, there exists... Let there be a grid, and let the area of ​​each grid be... The area of ​​the overlapping region is The fit is then defined as:

[0035] (6)

[0036] In the formula: For compatibility, The larger value between the inner core pile or the outer core pile and the grid wave velocity. It is the smaller value between the inner core pile or the outer core pile and the grid wave velocity. The area of ​​the high-precision grid;

[0037] S702, Assuming the inner core pile generates a non-coaxial angle ,Will From 0 to 15°, with a step size of 0.1°, the area where the inner core pile overlaps with the high-precision grid when it is non-coaxial is searched. The non-coaxial angle corresponding to the highest fit is the actual non-coaxial angle of the inner core pile.

[0038] S703, using the calculation process of S701 and S702, obtains the non-coaxial angle of the outer core pile. .

[0039] Secondly, the present invention proposes a computer device, including a processor and a memory, wherein the memory stores a computer program, and the processor executes the computer program to implement the method steps proposed in the present invention.

[0040] Furthermore, the present invention also proposes a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the method steps proposed in the present invention.

[0041] Compared with the prior art, the present invention has the following technical effects:

[0042] (1) When detecting the non-coaxial angle of the inner and outer core piles of the composite pile, the present invention makes full use of the hollow structure of the inner core pile itself, without the need for additional drilling to install sonic logging tubes, which significantly reduces the construction cost and reduces the disturbance to the pile foundation structure. It is a true non-destructive testing.

[0043] (2) The non-coaxial angle interpretation method proposed in this invention fills the gap in the prior art. For the first time, this invention calculates the non-coaxial angle between the inner core pile and the outer core pile through acoustic transmission non-destructive testing and interpretation algorithm, providing an efficient non-coaxial detection method. It can provide technical support and theoretical basis for the quality control of pile foundations, fill the current technical gap, and has important practical application value and innovation. Attached Figure Description

[0044] Figure 1 A schematic diagram of a field test using acoustic transmission method for interpreting the non-coaxial angle between the inner core pile and the outer core pile of a rigid composite pile, provided as an example of the present invention.

[0045] In the diagram: A and B are test piles, respectively; 1 is the acoustic emitter; 2 is the hollow part of the inner core pile; 3 is the inner core pile; 4 is the cement-soil pile; 5 is the acoustic receiver; the white part represents the soil around the pile; the rays represent the path of acoustic wave emission in the fan-shaped scanning method.

[0046] Figure 2 This is a schematic diagram of the grid division and path length of the present invention.

[0047] Figure 3 This is a low-precision grid wave velocity cloud map for the present invention.

[0048] Figure 4This is a high-precision grid wave velocity cloud map of the present invention.

[0049] Figure 5 This is an interpretation of the non-coaxial angle of the inner core pile in this invention.

[0050] Figure 6 This is an interpretation of the non-coaxial angle of the outer core pile in this invention. Detailed Implementation

[0051] To provide a more specific description of the objectives and implementation schemes of the present invention, the following will describe the implementation schemes of the present invention in more detail and completely with reference to specific embodiments. The embodiments described are only some, not all, of the embodiments of the present invention. It should be understood that the specific embodiments described below are for illustrative purposes only and are not intended to limit the scope of the present invention. Other embodiments obtained by those skilled in the art based on the described embodiments of the present invention without creative effort are all within the scope of protection of the present invention.

[0052] This invention proposes a method and system for interpreting the non-coaxial angle between the inner core pile and the outer core pile of a rigid composite pile. Through theoretical analysis and algorithm optimization, it provides an efficient non-coaxial detection method, which can provide technical support and theoretical basis for the quality control of pile foundations, fill the current technical gap, and has important practical application value and innovation.

[0053] The present invention will be described in detail below with reference to the accompanying drawings. It should be noted that the simplified process diagrams used in the drawings are not intended to limit the actual product. The present invention provides a method and system for interpreting the non-coaxial angle of the inner core pile and outer core pile of a reinforced composite pile, as follows... Figure 1 As shown, A and B are composite piles, with a sound transmitter in pile A and a sound receiver in pile B; 1 is the sound transmitter; 2 is the hollow part of the inner core pile; 3 is the inner core pile; 4 is the cement-soil pile; 5 is the sound receiver; the white part represents the soil around the pile; the rays represent the path of sound wave emission in the fan-shaped scanning method. The sound transmitter and receiver used in the measurement are arranged on the hollow surface.

[0054] This invention proposes a method for interpreting the non-coaxial angle between the inner core pile and the outer core pile of a reinforced composite pile, specifically including:

[0055] Step 1: Select any two adjacent composite piles A and B, and fill the inner core of the composite piles with water;

[0056] Step 2: Measure the thickness of the inner core pile and outer core pile of the reinforced composite pile using a measuring tape. and The pile spacing is ;

[0057] Step 3: The acoustic transmitter and receiver are placed in piles A and B respectively. The acoustic receiver and transmitter are arranged at equal depths of 0.2m to 1.0m, starting from 1m underground.

[0058] Step 4: Following the principle of the sector scanning method, raise the sound transmitters upwards sequentially. At each depth, transmit a signal to all receivers once until reaching the top of pile A. Then, store the arrival times of all sound waves from the receivers in the instrument.

[0059] Step 5: Grid-based wave velocity inversion step, specifically as follows: Divide the plane between the transmitter and receiver into grids. There are grids, assuming there are... The arrival time of the experimental sound waves. The wave velocity in the grid is calculated according to the following set of equations:

[0060] (1)

[0061] In the formula, For the ray in the first OK The path length in the column grid can be calculated once the grid division is determined; For the first Line number The reciprocal of the wave velocity of the column grid is the unknown that needs to be calculated; For the first The arrival time of the test sound wave of the ray can be determined through step 4;

[0062] Step 6: Calculate the path matrix based on the algorithm. The principle is as follows: denote the line segment of each ray in each grid as a vector, and write code to traverse all the vectors in the grid in the order of the rays, where the magnitude of the vector is the path of the ray in the grid. ;

[0063] Step 7: Calculate the initial value of the reciprocal of the wave velocity of the mesh. Assuming the wave velocity of each ray is uniform, the arrival time of the sound wave for each ray is distributed among the grid cells according to the path length of each grid cell. The reciprocal of the initial wave velocity is obtained by dividing the total time of all rays passing through the grid cells by the total path length within the grid. The calculation formula is as follows:

[0064] (2)

[0065] Step 8: Initial value of the reciprocal of the wave velocity calculated in Step 7 and the arrival time of the experimental sound waves obtained in step 4 The reciprocal of the wave speed is obtained by solving the system of equations in formula (1) using the following iterative algorithm. The iterative equation is as follows:

[0066] (3)

[0067] In the formula, The relaxation factor is set to 1.2.

[0068] For the first Line number List a grid in The reciprocal of the grid wave velocity after the next iteration;

[0069] For the first Line number List a grid in The reciprocal of the grid wave velocity after the next iteration;

[0070] For the first The arrival time of the test sound wave of the ray;

[0071] For the first The first ray OK Path length in column grid;

[0072] For the first Line number The number of rays passing through the column grid.

[0073] Step 9: Take the reciprocal of the grid wave velocity in formula (3) Substituting into the following formula yields the low-precision grid wave velocity. sum matrix :

[0074] (4)

[0075] Step 10: High-precision grid wave velocity inversion, the steps are as follows:

[0076] Based on the known low-precision grid wave velocity matrix (Formula 4), with a step size of 5cm or 10cm, reset the positions of the sound emission point and the sound reception point downward from the depth of 0m, and calculate the arrival time of the amplified sound wave according to Formula (5):

[0077] (5)

[0078] Step 11: Based on the already calculated Repeat steps 5-9 to calculate a high-precision grid wave velocity matrix. .

[0079] Step 12: Interpretation of non-coaxial angles, as follows:

[0080] Based on the mesh wave velocity cloud map, it was determined that the inner and outer core piles are not coaxial. When the inner and outer core piles are not coaxial, there will be an overlap with the calculated high-precision mesh wave velocity matrix. Within this overlapping region, there are... Let there be a grid, and let the area of ​​each grid be... The area of ​​the overlapping region is The fit is then defined as:

[0081] (6)

[0082] In the formula: For compatibility, The larger value between the inner core pile or the outer core pile and the grid wave velocity. It is the smaller value between the inner core pile or the outer core pile and the grid wave velocity. For high-precision grid area.

[0083] Step 13: Assume that the inner core pile generates a non-coaxial angle. ,Will From 0 to 15°, with a step size of 0.1°, the search is performed on the overlapping area between the inner core pile and the high-precision mesh when they are not coaxial. The non-coaxial angle corresponding to the highest fit is the actual non-coaxial angle of the inner core pile. Non-coaxial angle of outer core pile The calculation process is the same.

[0084] Example 1: This example provides a method for interpreting the non-coaxial angle between the inner core pile and the outer core pile of a reinforced composite pile. The on-site acoustic transmission test is as follows: Figure 1 As shown, the specific steps include the following:

[0085] Step 1: Select any two adjacent composite piles A and B, and fill the inner core of the composite piles with water;

[0086] Step 2: Measure the thickness of the inner core pile and outer core pile of the reinforced composite pile using a measuring tape. and The pile spacing is ;

[0087] Step 3: The acoustic transmitter and receiver are placed in piles A and B respectively, with the acoustic receiver and transmitter arranged at equal depths of 1.0m, starting from 1m underground.

[0088] Step 4: Following the principle of the sector scanning method, raise the sound transmitters upwards sequentially. At each depth, transmit a signal to all receivers once until reaching the top of pile A. Then, store the arrival times of all sound waves from the receivers in the instrument.

[0089] The measurement depth in Example 1 was 14m, and the measured sound wave arrival time data... See Table 1.

[0090] Table 1. Sound wave arrival time in Example 1 experiment.

[0091]

[0092] Step 5: As Figure 2 As shown, a grid is created, dividing the plane between the transmitter and receiver into... There are grids, assuming there are... The arrival time of the experimental sound waves. The wave velocity in the grid is calculated according to the following set of equations:

[0093] (1)

[0094] In the formula, For the ray in the first OK The path length in the column grid can be calculated once the grid division is determined; For the first Line number The reciprocal of the wave velocity of the column grid is the unknown that needs to be calculated; For the first The arrival time of the test sound wave of the ray can be determined through step 4;

[0095] Step 6: Calculate the path matrix based on the algorithm. The principle is as follows: denote the line segment of each ray in each grid as a vector, and write code to traverse all the vectors in the grid in the order of the rays, where the magnitude of the vector is the path of the ray in the grid. ;

[0096] Step 7: Calculate the initial value of the reciprocal of the wave velocity of the mesh. Assuming the wave velocity of each ray is uniform, the arrival time of the sound wave for each ray is distributed across the grid according to the path length of each grid cell. The reciprocal of the initial wave velocity is obtained by dividing the total time of all rays passing through the grid cell by the total path length within the grid. The calculation formula is as follows:

[0097] (2)

[0098] The initial value of the reciprocal of the grid wave velocity During the calculation, if no rays pass through the grid, the reciprocal of its wave velocity is tentatively set to... .

[0099] Step 8: Initial value of the reciprocal of the wave velocity calculated in Step 7 and the arrival time of the experimental sound waves obtained in step 4 The reciprocal of the wave speed is obtained by solving the system of equations in formula (1) using the following iterative algorithm. The iterative equation is as follows:

[0100] (3)

[0101] In the formula, The relaxation factor is set to 1.2. For the first Line number List a grid in The reciprocal of the grid wave velocity after the next iteration; For the first Line number List a grid in The reciprocal of the grid wave velocity after the next iteration; For the first The arrival time of the test sound wave of the ray; For the first The first ray OK Path length in column grid; For the first Line number The number of rays passing through the column grid, and the condition for termination of the iteration in formula (3) is the number of iterations. Or iterative standard deviation of sound wave arrival time .

[0102] Step 9: As Figure 3 As shown, the reciprocal of the grid wave velocity in formula (3) Substituting into the following formula yields the low-precision grid wave velocity. sum matrix :

[0103] (4)

[0104] The steps for high-precision grid wave velocity inversion are as follows:

[0105] Step 10: Based on the known low-precision grid wave velocity matrix (Formula 4), with a step size of 10cm, reset the positions of the sound emission point and the sound reception point downward from the depth of 0m, and calculate the arrival time of the amplified sound wave according to Formula (5):

[0106] (5)

[0107] Step 11: Based on the already calculated Repeat steps 5-9 to calculate a high-precision grid wave velocity matrix. . Figure 4 This is a high-precision grid wave velocity cloud map of the present invention.

[0108] The method for interpreting non-coaxial angles is as follows:

[0109] Step 12: Based on the grid wave velocity cloud map, determine if the inner and outer core piles are non-coaxial. When the inner and outer core piles are non-coaxial, there will be overlap with the calculated high-precision grid wave velocity matrix. Within this overlapping area... Let there be a grid, and let the area of ​​each grid be... The area of ​​the overlapping region is The fit is then defined as:

[0110] (6)

[0111] In the formula: For compatibility, The larger value between the inner core pile or the outer core pile and the grid wave velocity. It is the smaller value between the inner core pile or the outer core pile and the grid wave velocity. For high-precision grid area.

[0112] Step 13: As Figure 5 and Figure 6 As shown, it is assumed that the inner core pile generates a non-coaxial angle. ,Will From 0 to 15°, with a step size of 0.1°, the non-coaxial angle corresponding to the highest fit when the core pile tilts and coincides with the high-precision grid wave velocity is the actual non-coaxial angle of the core pile. Non-coaxial angle of outer core pile The calculation process is the same. According to the steps described in Examples 1-4, the non-coaxial angles of the inner core pile and the outer core pile are calculated as follows: and .

[0113] Example 2: This example proposes an electronic system, including: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, the instructions being executed by the at least one processor to enable the at least one processor to perform the method steps of the present invention.

[0114] Example 3: This example proposes a computer-readable storage medium storing a computer program thereon. When the computer program is executed by a processor, it implements the steps of the method described in this invention, which will not be repeated here.

[0115] It should be noted that the processing flow of embodiments 2-3 corresponds to the specific steps of the method provided in the embodiments of the present invention, and has the corresponding functional modules and beneficial effects of the method. Technical details not described in detail in this embodiment can be found in the method provided in the embodiments of the present invention.

[0116] The program code used to implement the methods of this application may be written in any combination of one or more programming languages. This program code may be provided to a processor or controller of a general-purpose computer, special-purpose computer, or other programmable data processing device, such that when executed by the processor or controller, the functions / operations specified in the flowcharts and / or block diagrams are implemented. The program code may be executed entirely on a machine, partially on a machine, as a standalone software package partially on a machine and partially on a remote machine, or entirely on a remote machine or server.

[0117] In the context of this application, a machine-readable medium can be a tangible medium that may contain or store a program for use by or in conjunction with an instruction execution system, apparatus, or device. A machine-readable medium can be a machine-readable signal medium or a machine-readable storage medium. Machine-readable media can be, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination of the foregoing. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fibers, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.

[0118] The above embodiments are merely illustrative of the technical concept of the present invention and should not be construed as limiting the scope of protection of the present invention. Any modifications made to the technical solutions based on the technical concept proposed in this invention shall fall within the scope of protection of this invention.

Claims

1. A method for interpreting the non-coaxial angle between the inner core pile and the outer core pile of a reinforced composite pile, characterized in that, The steps are as follows: Step 1: Select any two adjacent composite piles A and B, and fill the inner core pile of the composite pile with water; Step 2: Measure the thickness of the inner core pile, the thickness of the outer core pile, and the pile spacing of the composite pile; Step 3: Place the acoustic transmitter and receiver in the inner core piles of pile A and pile B respectively, starting from the preset underground depth and arranging the points downwards at equal intervals. Step 4: Following the sector scanning method, raise the acoustic emitter sequentially, transmit signals to all acoustic receivers at each depth point, and record the arrival time data of the sound waves arranged in column vectors according to the depth from low to high. Step 5: Based on the arrival time of the sound waves, divide the plane between the transmitter and receiver into multiple grids, construct a path matrix, and calculate the grid wave velocity using an inversion algorithm; Step 6: Calculate the high-precision grid wave velocity matrix using a high-precision grid wave velocity inversion algorithm based on the grid wave velocity. Step 7: Based on the high-precision grid wave velocity matrix, calculate the non-coaxial angle between the inner core pile and the outer core pile using the fit formula and interpretation algorithm.

2. The method according to claim 1, characterized in that, In step 1, the inner core pile is filled with water as a coupling agent for acoustic emission and reception, and there is no obvious silt suspended in the water.

3. The method according to claim 1, characterized in that, Step 2 involves taking multiple measurements and averaging them.

4. The method according to claim 1, characterized in that, Step 5 uses an inversion algorithm and a joint algebraic reconstruction method to iteratively solve for the grid wave velocity. The iteration terminates when the number of iterations reaches a preset value or the standard deviation of the sound wave arrival time is less than a threshold. The specific steps are as follows: Divide the plane between the transmitter and receiver into... There are grids, assuming there are... The arrival time of the experimental sound waves and the wave velocity in the grid are calculated using the following set of equations: (1) In the formula, For the ray in the first OK Path length in column grid For the first Line number The reciprocal of the wave velocity of the grid; For the first The arrival time of the test sound wave for the ray is n less than m.

5. The method according to claim 1, characterized in that, In step 6, the high-precision grid wave velocity inversion algorithm includes resetting the positions of the acoustic emission point and the acoustic reception point with a preset step size, calculating the arrival time of the expanded acoustic wave based on the grid wave velocity matrix, and then calculating the high-precision grid wave velocity matrix through the inversion algorithm. Specifically, the line segment of each ray in each grid is recorded as a vector, and code is written to traverse all vectors in the grid sequentially according to the ray order. The magnitude of each vector is the path of the ray in the grid. .

6. The method according to claim 5, characterized in that, The high-precision grid wave velocity matrix is ​​calculated using a high-precision grid wave velocity inversion algorithm, as detailed below: S601. Calculate the initial value of the reciprocal of the wave velocity of the grid. Assuming the wave velocity of each ray is uniform, the arrival time of the sound wave of each ray is distributed among the grids according to the path length of each grid cell. The reciprocal of the initial wave velocity is obtained by dividing the total time of all rays passing through the grid cells by the total path length within the grid. The calculation formula is as follows: (2) S602, based on the initial value of the reciprocal of the calculated wave velocity. and the arrival time of the test sound waves The reciprocal of the wave speed is solved by calculating the system of equations in formula (1) using the following iterative algorithm. The iterative equation is as follows: (3) In the formula, The relaxation factor is selected from 1.0 to 1.

2. For the first Line number List a grid in The reciprocal of the grid wave velocity after the next iteration; For the first Line number List a grid in The reciprocal of the grid wave velocity after the next iteration; For the first The arrival time of the test sound wave of the ray; For the first The first ray OK Path length in column grid; For the first Line number The number of rays passing through the column grid; S603, the reciprocal of the grid wave velocity in formula (3) Substituting into the following formula yields the low-precision grid wave velocity. sum matrix : (4) S604. Based on the obtained grid wave velocity matrix With a step size of 5cm or 10cm, the positions of the sound emission point and the sound reception point are reset from a depth of 0m downwards. The arrival time of the amplified sound wave is calculated according to formula (5): (5) S605, based on the calculated Repeat steps 5 and 6 to calculate the high-precision grid wave velocity matrix. 。 7. The method according to claim 1 or 6, characterized in that, The interpretation method for non-coaxial angles is as follows: S701. Based on the grid wave velocity cloud map, it is determined that the inner core pile and outer core pile are not coaxial. When the inner core pile and outer core pile are not coaxial, there will be an overlap with the calculated high-precision grid wave velocity matrix. Within the overlapping area, there exists... Let there be a grid, and let the area of ​​each grid be... The area of ​​the overlapping region is The fit is then defined as: (6) In the formula: For compatibility, The larger value between the inner core pile or the outer core pile and the grid wave velocity. It is the smaller value between the inner core pile or the outer core pile and the grid wave velocity. The area of ​​the high-precision grid; S702, Assuming the inner core pile generates a non-coaxial angle ,Will From 0 to 15°, with a step size of 0.1°, the area where the inner core pile overlaps with the high-precision grid when it is non-coaxial is searched. The non-coaxial angle corresponding to the highest fit is the actual non-coaxial angle of the inner core pile. S703, using the calculation process of S701 and S702, obtains the non-coaxial angle of the outer core pile. .

8. The method according to claim 6, characterized in that: The iteration of formula (3) terminates when the number of iterations has been reached. or the standard deviation of the arrival time of sound waves during iteration. .

9. A computer device comprising a processor and a memory, the memory storing a computer program, characterized in that, When the processor executes the computer program, it implements the method as described in any one of claims 1 to 8.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the method as described in any one of claims 1 to 8.

Images