Parallel finite layer method for land subsidence and its homotopy parameter inversion

Abstract

Land subsidence is caused by the interaction between groundwater flow and soil deformation. At present, in the field of land subsidence calculation, there is an increasing tendency to use Biot's consolidation theory to precisely simulate the coupling processes of flow and deformation. However, there are some practical problems with large-scale numerical computation, such as large computational loads, low efficiency, and difficulty in obtaining parameters. The Finite Layer Method employs orthogonality of analytical functions to achieve decoupling among various series terms. And the method can be enhanced through parallel computing to optimize the efficiency of solving complex and extensive problems due to its decoupling property. Based on this, combined with the nonlinear homotopy method, a new approach for the reversal of land subsidence is proposed. The validity of the parallel program is confirmed through comparison with established solutions, and the impact of various factors on the efficiency of parallel computing is analyzed. The results indicate that the parallel approach enhances computational efficiency, while the homotopy inversion technique exhibits broad convergence characteristics.