Application of Chemical Dynamics Method in Solving HydroGeological Parameters

With the increasing of coal mining depth, the coal seam floor is threatened more and more seriously with the limestone aquifer of Taiyuan Formation. In order to make the coal seam floor water inrush accident to be effectively prevented, it is very important to ascertain the hydro-geological condition of the research area and adequately gain hydro-geological parameters. Based on the theory of chemical dynamics, the chemical dynamical equation of the mineral was obtained. The chemical dynamical equation of the mineral and the Darcy’s Law were combined to derive the expression of permeability coefficient and transmissibility coefficient that was indicated by hydro-chemistry index. Based on the data of water quality analysis in the research area, the hydro-geology parameter of the aquifer of Taiyuan Formation was calculated. And it was compared with what the hydro-geology parameter was determined by the data of pumping test. The results show that the permeability coefficient (K) of the research area is 1.7512 m/d, and transmissibility coefficient (T) is 42.0282 m/d. The calculated results of chemical dynamics and pumping test results are little different between them. And it has certain practical significance in the mine water control. It provides an effective method for obtaining hydro-geological parameters.


Introduction
With the increasing demand for coal resources in China, the work of coal mining is gradually developed into the deep mining depth and complex geological conditions [1] [2] [3] [4].The continuous increase of mining depth has become the main method of mining.Along with the continuous increase of mining depth, the Mining under Pressure has become the main method of mining.And the threat of coal seam mining will gradually aggravate from the depth of the aquifer of Taiyuan Formation.The safety of coal seam floor is very important [1] [2] [3].The main reason is the complex geological and hydro-geological conditions of aquifer and insufficient understanding of hydro-geological conditions.
Meanwhile, it is difficult to obtain the basic hydro-geological parameters.Therefore, it is of great practical significance to identify the aquifer properties and hydro-geological parameters of the coal seam floor before mining, which is of the most importance for the prevention and control of water in coal mines [3] [4] [5] [6].
For a long time, many methods were usually used in the process of obtaining hydro-geological parameters.For example, pumping test, slug test, water pressure test and so on [4] [5].According to Zhang H.Y. et al. [6], taking the South deep mining area of the Chensilou coal as the research object, the hydro-geological parameters of limestone aquifer were calculated by unsteady flow dewatering test.According to Ju X.M. et al. [7], the pumping tests and slug tests at the same time in the Laosui river of Huaibei in Anhui province were carried out, and the results of the two methods were compared and analyzed.According to Hu W.Y. et al. [8], a simple method of calculating hydro-geological parameters was introduced by the observed dynamic data of ground water level for a long time.These methods of obtaining hydro-geological parameters could be summarized as hydrodynamic methods [4]- [9].In recent years, experts and scholars have put forward the method of chemical dynamics to obtain hydro-geological parameters [10]- [15] and the quantitative evaluation of hydro-geological conditions in the study area is realized.
On this basis, this paper took a mine in Anhui province as an example.The expression of hydro-geological parameters was derived by using the method of Chemical Dynamics.And the hydro-geology parameter of the limestone aquifer of Taiyuan Formation was calculated.The data of pumping test in the research area was collected and arranged.The calculated results of chemical dynamics and pumping test results were compared to analyze the reason for the deviation.

Determination of Ion Activity
In the ground water, part of the ions does not work in the reaction process due to the interaction of various ions.This leads to a difference between the effective concentration of the ions and the true concentration of the ions in the process of chemical reaction [10] [11] [12] [14].Therefore, it is necessary to correct the measured concentration in the ground water.The corrected concentration is known as the Ion Activity and is represented by α i .The Relationship between Ion Activity (α i ) and the true concentration of the ions (c i ) is as follows [10]: where, γ i is the activity coefficient.
At this time the theory of Debye and Hükel is used to explain difference between the ground water solution and ideal solution.And the formula for calculating the activity coefficient of an ion solution by mathematical method is as follows [10] [11]: where, A and B are constants that depend on the temperature of the water.i a is constant that depends on the ionic radius.I is the ionic strength of the solution.
The ionic strength (I) is a measure of the electrostatic field that formed by various ions in the solution.And the formula for calculating the ionic strength (I) is as follows: where, m i is the substance concentration of the ion.z i is the Valence number of ion.M is the kinds of ion.
The theory of Debye and Hükel was applied to the temperature of 0 to 60.The values of A and B are as follows.
According to the data in Table 1, the curve fitting of T with A and B was made.The curve fitting is as follows (Figure 1).

Determination of Saturation Index
The Saturation index is a parameter to determine the state of water and minerals, and is an important means to evaluate the formation and development of karst.It is expressed by the symbol of SI.When SI > 0, it means that the mineral is in a super-saturation state.When SI = 0, it means that the mineral is in a saturation state.When SI < 0, it means that the mineral is in a dissolved condition.The multi mineral system refers to the dissolution of calcite, dolomite, and gypsum in the CO 2 -H 2 O system.Under these circumstances, the saturation index is calculated by the mineral composition of ion activity and solubility product constant.The saturation index of Calcite, Dolomite and Gypsum are as follows: where, α i is the ion activity of i, K 1 is the first-order dissociation constant of carbonic acid, K 2 is the second-order dissociation constant of carbonic acid, K c is the dissociation constant of Calcite, K d is the dissociation constant of Dolomite, K g is the dissociation constant of Gypsum.

The Darcy's Law Is Expressed by Chemical Indexes
If the component of i is generated through many minerals, the equation of the mineral dissolution process can be superimposed.In some monographs [10] [12] [16] [17] [18], according to the theory of chemical dynamics, the chemical dynamics equation for the dissolution of minerals was derived.And the equation is as follows: ( ) The Equation ( 4) was deformed as follows: ( ) In the process of ground water movement, the ds i /dt can be considered as the infiltration velocity (V ni ) on the groundwater flow line.If V ni is substituted into Equation ( 5), it will get the following relationship: ( ) ( ) ( ) where, j is the number of dissolved mineral, i is the component of ion from mineral dissolution, V ni is the infiltration velocity of the component of i on the groundwater flow line, α ij is the ion activity of the component of i from the mineral of j, S i is the flow distance of the component of i along with the stream line, v ij is the stoichiometric number of the component of i that dissolve from the mineral of j, k j is the chemical dynamics constant for the dissolution of mineral of j.
Because of the ion activity (α i ) and the saturation index (β j ) are the functions of time and space, the Equation ( 7) can be integrated.It will get the following relationship: ( ) The Equation ( 8) can be simplified by using the trapezoid formula.The Equation (8) was deformed as follows: ( ) where, the ion activity of the component of i from A to B, β jA is the saturation index of the mineral of j at A, β jB is the saturation index of the mineral of j at B.

Expression of Hydro-Geological Parameters
If Equation ( 9) is substituted into the expression of the Darcy's Law that is expressed by Hydraulic index, it will get the following relationship: ( ) where, K AB is the permeability coefficient, h ∆ is the height deviation between A and B. M is thickness of the limestone aquifer of Taiyuan Formation.

Expression of Chemical Dynamics Constant
There are two methods for obtaining chemical dynamics constants, which are laboratory test and field test.The chemical dynamics constants were obtained by taking field experiments in this paper.And the equation is as follows: ( ) where, K is the permeability coefficient that was obtained by hydraulic index.
According to the results of pumping test in the study area, the K is equal to 1.3913 m/d.
On the basis of the chemical dynamics constants was obtained by predecessors, the chemical dynamics constants of Calcite, Dolomite and Gypsum were calculated respectively.The corresponding permeability coefficient (K) and transmissibility coefficient (T) were obtained by using the chemical dynamics constants of Calcite, Dolomite and Gypsum.Then the average of permeability coefficient (K) and transmissibility coefficient (T) were calculated.
The chemical dynamics constants of Calcite (c), Dolomite (d) and Gypsum (g) are expressed as follows: ( ) ( )

The Calculation of Hydro-Geological Parameters by Chemical Indexes
The research area is located in the south of the mine field.The main aquifers include the first to fourth aquifer of Cenozoic Erathem, the sandstone aquifers in coal measures and the limestone aquifer of Taiyuan Formation.The aquifer in the research area is shown in Figure 2. The limestone aquifer of Taiyuan Formation is one of the major hidden dangers of mine filling water.And it is of great significance to research the limestone aquifer of Taiyuan Formation.
A total of 10 water samples were taken at 1 to 5 # sampling points in the research area, and there are 2 water samples at each sampling point.The 10 water samples were used for water quality analysis.The distributions of sampling points are shown in Figure 3.The results of water quality analysis are shown in Table 2.
According to the geological report of the research area, the thickness of the limestone aquifer of Taiyuan Formation is 24 meters.So, the chemical dynamics constant could be calculated by Equation (14), Equation (15) and Equation (16).
The results are shown in Table 3.
Table 3.The chemical dynamics constant.
No.  The contour map of the permeability coefficient in the research area is shown as follows.

Evaluation of Hydro-Geological Conditions
It can be seen from Figure 4 that the permeability coefficient of the limestone aquifer of Taiyuan Formation in the whole research area has little change in general.According to the results of the hydro-geological parameters in

Conclusions
1) The results show that the limestone aquifer of Taiyuan Formation has strong water permeability and medium water conductivity, and the hydro-geological parameters of the whole research area have not changed much.2) The little difference between the hydro-geological parameters was obtained by the chemical dynamics method and was obtained from the pumping test.It shows that it is feasible to evaluate hydro-geological condition from the point of view of chemical dynamics.And the chemical dynamics has a certain application prospect.

Figure 1 .
Figure 1.The curve fitting of T with A and B.

3 )
In order to improve the accuracy of the calculation of hydro-geological parameters by chemical dynamics method, the following points should be made in the following work: a) More water samples are collected for water quality analysis; b) The sampling points should be evenly distributed in time and space; c) The retrieved water samples should be immediately analyzed for water quality to avoid evaporation of water and affect the results of water quality analysis; d) Water samples were collected for water quality analysis to reduce test errors.

Table 1 .
The value of A and B in different temperature.
The Value of A and B

Table 4 .
The calculation of hydro-geological parameters