地下工程渗流断层数值模拟的隐式复合材料单元法研究

结合地下工程渗流场三维有限元分析复杂性的实际,对断层的处理进行积极探讨。在拟连续介质模型的基础上,充分利用有限元法的特点,首次提出用隐式复合材料单元法对主要断层进行模拟。算例验证了该法的正确性和实用性。

Study on “triangle” water-inrush mode of strong water-guide collapse column

The water-inrush mechanism of strong water-guide collapse column in coal seam is studied based on the establishment of geological and mathematical models of "triangle" water-inrush mode. The geological background of Shuangliu mine is considered a prototype, similar simulation tests are adopted to analyze the water-inrush rules under this model, and the formation of water-guide channel and water-inrush process is investigated by examining the changes in rock resistivity. This work also uses the coupled cloud image derived from numerical simulation software to verify the results of simulation test. Results show that the numerical simulation of "triangle" water-inrush mode is consistent with the similar simulation. The "triangle" seepage area, which is located at the bottom of collapse columns and is connected to aquifer, is caused by the altered seepage direction and strengthened seepage actions after the overlapping of hydraulic transverse seepage in collapse column and hydraulic vertical seepage flow in aquifer. Under "triangle"water-inrush model, water-guide channel is formed by the communication between plastic failure zone of working face baseplate and"triangular" seepage area. Accordingly, the threatening water-inrush distance between working face and collapse column increases by 20 m compared with that of theoretical calculation.

Study on implicit composite element method for seepage analysis in underground engineering

Based on the basic principle of the finite element method, the implicit composite element method for numerical simulation of seepage in underground engineering is proposed. In the simulation, the faults and drainage holes are set implicitly in the model elements without adding additional elements. Elements containing fault or drainage-hole data are termed composite elements. Then, their information data in model could be obtained. By determining the osmotic transmission matrix of the composite elements, the permeability coefficient matrix is then obtained. The method was applied to the numerical simulation of the seepage field around the underground powerhouse of the Ganhe Pumping Station in Yunnan, China, using a compiled three-dimensional finite element method calculation program. The rock mass around the site includes two faults. The seepage field in the rock mass was analyzed at different stages of the engineering project. The results show that, before the excavation of the underground caverns, the rock mass seepage is affected by the faults and the groundwater permeated down along the tangential fault plane. After the excavation of the caverns during the operation period, the groundwater is basically drained away and the underground caverns are mostly above the groundwater level. Thus, the calculation results of the engineering example verify the implicit composite method for the simulation of faults and drainage holes. This method can well meet the calculation demands of practical engineering.

Percolation Phase Transitions from Second Order to First Order in Random Networks

We investigate a percolation process where an additional parameter q is used to interpolate between the classical Erd¨os–R′enyi(ER) network model and the smallest cluster(SC) model. This model becomes the ER network at q = 1, which is characterized by a robust second order phase transition. When q = 0, this model recovers to the SC model which exhibits a first order phase transition. To study how the percolation phase transition changes from second order to first order with the decrease of the value of q from 1 to 0, the numerical simulations study the final vanishing moment of the each existing cluster except the N-cluster in the percolation process. For the continuous phase transition,it is shown that the tail of the graph of the final vanishing moment has the characteristic of the convexity. While for the discontinuous phase transition, the graph of the final vanishing moment possesses the characteristic of the concavity.Just before the critical point, it is found that the ratio between the maximum of the sequential vanishing clusters sizes and the network size N can be used to decide the phase transition type. We show that when the ratio is larger than or equal to zero in the thermodynamic limit, the percolation phase transition is first or second order respectively. For our model, the numerical simulations indicate that there exists a tricritical point qcwhich is estimated to be between0.2 < qc< 0.25 separating the two phase transition types.

Simulation procedure of unconfined seepage with an inner seepage face in a heterogeneous field

An inner seepage face phenomenon is given and a numerical simulation procedure has been developed.It may appear at the interface of two materials when an unconfined seepage flows from a porous media to a coarser porous media with a higher permeability.Inaccuracy and divergent problems may arise both in a saturated-only and in a variably saturated analysis while an inner seepage face is not simulated with a special procedure.The position of the seepage face is determined during the nonlinear iteration process and the flux of the inner seepage face nodes is transferred to the downstream side nodes.Validity and efficiency of the procedure are illustrated by the simulation of two dimensional steady state seepage examples of heterogeneous zoned dams which is usually used to validate algorithms.An analysis of a three-dimensional earth core rockfill dam is also presented here.The procedure can also be applied to general transient seepage problems.

Threshold pressure gradient of fluid flow through multi-porous media in low and extra-low permeability reservoirs

After analyzing many studies of fluid flow theory of multi-porous media in low and extra-low permeability reservoirs and the numerical simulation of non-Darcy flow, we found that a negative flow rate occurs in the existing non-Darcy flow equation, which is unreasonable. We believe that the existing equation can only be considered as a discriminant to judging Darcy flow or non-Darcy flow, and cannot be taken as a fluid flow governing equation of multi-porous media. Our analysis of the experimental results shows that the threshold pressure gradient(TPG) of low and extra-low permeability reservoirs is excessively high, and does not conform to fluid flow through multi-porous media in the actual reservoir situation. Therefore, we present a reasonable TPG ranging from 0.006 to 0.04 MPa/m at the well depth of 1500 m and oil drainage distance of 500 m. The results of our study also indicate that the non-Darcy flow phenomenon will disappear when the TPG reaches a certain value. In addition, the TPG or non-Darcy flow in low and extra-low permeability reservoirs does not need to be considered in the productivity prediction and reservoir numerical simulation. At present, the black oil model or dual-porous media is suitable for simulating low and extra-low permeability reservoirs.

Slip-enhanced Reverse Electrodialytic Power Generation in Ion-selective Nanochannels

Power generation by reverse electrodialysis in ion-selective nanochannels is numerically investigated. Especially,in the present study, the influence of hydrodynamic slip at the surface of nanochannels is investigated. The current-potential characteristics of the nanochannels are calculated by solving several governing equations:Nernst-Planck equation for the ionic concentrations, the Poisson equation for the electric potential, and the Navier-Stokes equation for the diffusioosmotic flow. Hydrodynamic slip is applied as the boundary condition at the surface of nanochannels. As the slip length increases, the diffusioosmotic flow velocity and electrical conductance of ions increase because the friction at the surface of nanochannels decreases. It is shown that the power generation is enhanced by 44% with a moderate 100nm slip length by using a nanochannel with 10nm height.

Numerical investigation on critical length of impermeable plate below underwater pipeline under steady current

This study proposes a new approach in which an impermeable plate is placed under the pipeline to prevent the local scour around the pipeline.In order to understand the performance of this approach,the finite volume method(FVM) and volume of fluid(VOF) method are adopted to simulate the flow field around the pipeline.The pressure distribution along the sandy bed surface is obtained by considering the variation of water surface.Furthermore,the effects of water depth,unidirectional and bidirectional impermeable plates on pressure difference are discussed.The seepage flow field of sandy bed near underwater pipeline is numerically simulated using the laminar and porous media model.On this basis,the effect of the impermeable plate length on hydraulic gradient is investigated and the critical length of impermeable plate is obtained.The simulated results show that when the water depth is smaller than 5.00D(D is the diameter of pipeline),the effect of the water depth on the pressure difference is remarkable.The pressure differences between two endpoints of both the unidirectional and bidirectional plates decrease with the increase of the plate length.The variations of the pressure differences for both the unidirectional and bidirectional plates are similar.With the increase of plate length,the hydraulic gradient decreases and the piping at the seepage exit is avoided effectively as long as it reaches a certain length.Such a critical length of the plate decreases with the increase of the water depth.When water depth is larger than 4.00D,the effect of the water depth on the critical length is small.For the same water depth,the critical length of impermeable plate increases with the increase of the dimensionless flow parameter.Numerical simulation results are in good agreement with the available experimental measurements.