Estimates of Horizontal Groundwater Flow Velocities in Boreholes

Currently, monitoring tools can be deployed in observation boreholes to better assess groundwater flow, flux of dissolved contaminants and their mass discharge in an aquifer. The relationship between horizontal water velocity in observation boreholes and Darcy fluxes in the surrounding aquifer has been studied for natural flow conditions (i.e., no pumping). Interpretation of measurements taken with dilution tests, the colloidal borescope, the Heat Pulse Flowmeter, and other techniques require the conversion of observed borehole velocity u to aquifer Darcy flux q_∞ . This conversion is typically done through a proportionality factor α = u/q_∞ . In experimental studies as well as in theoretical developments, reported values of α vary almost three orders of magnitude (from 0.5 to 10). This large variability in reported values of α could be explained by: (1) unclear distinction between Darcy flux and water seepage velocity, (2) unclear definition of water velocity in the borehole, (3) effects of well screen and the presence of the measurement device itself on the observable velocities, and (4) hydraulic conditions in the borehole annulus. We address (1), (2) from a conceptual/theoretical perspective, and (3) by means of numerical simulations. We show that issue (1) in low porosity aquifers can yield to order-of-magnitude discrepancies in estimates of q_∞ ; (2) may result in discrepancies of up to 50%, and (3) can cause differences up to 20% of water velocity in the borehole void space compared to the theoretical case of an open borehole.     

Hyporheic exchange mechanism driven by flood wave

To study the hyporheic exchange driven by a single peak flood-induced water level fluctuation (i.e., flood wave), a method combining numerical simulation with theoretical derivation was proposed based on the Inbuk Stream, Korea, where flooding occurs frequently. The hyporheic exchanges induced by different flood waves were investigated by varying amplitude (A), duration (T), wave type parameter (r), and rising duration (t_p), which were adopted from the real-time stream stage fluctuations. Additionally, the idea of constant upstream flood volume (CUFV) condition for flood waves was put forward, and the effects of “Botan” (T/A) and peak number (N) on hyporheic exchange were studied. The results showed that the hyporheic exchange flux (q) was controlled by the water level h (sine-type) and its change rate v (cosine-type), and was proportional to the polynomial of them q ∝ (ω∙ h + v), where ω is the angular frequency of the flood wave. Based on this mechanism, the influence principles on hyporheic exchanges of the typical flood wave parameters (A, T, r and t_p) as well as T/A and N under CUFV condition were clarified. The main characteristic variables of hyporheic exchange, which were maximum aquifer storage and residence time, were positively correlated. They also had positive relations to the integral of the flood wave over time, which increased when the wave became higher, wider, rounder and less skewed. However, when CUFV condition was imposed, the residence time was positively correlated with T/A, whereas the maximum aquifer storage was negatively correlated with T/A. With the increase in N, water exchanged more frequently and some water returned to the stream early, leading to the slight decrease in maximum aquifer storage and residence time. These findings enriched the theory of hyporheic exchange driven by surface water fluctuation and be of great significance to enhance pollutant degradation in the hyporheic zone downstream of reservoirs.     

Effects of inertial and kinematic interaction on seismic behavior of pile with embedded foundation

Effects of inertial and kinematic forces on pile stresses are studied based on large shaking table tests on pile-structure models with a foundation embedded in dry and liquefiable sand deposits. The test results show that, if the natural period of the superstructure, T_b, is less than that of the ground, T_g, the ground displacement tends to be in phase with the inertial force from the superstructure, increasing the shear force transmitted to the pile. In contrast, if T_b is greater than T_g, the ground displacement tends to be out of phase with the inertial force, restraining the pile stress from increasing. With the effects of earth pressures on the embedded foundation and pile incorporated in, pseudo-static analysis is conducted to estimate maximum moment distribution in pile. It is assumed that the maximum moment is equal to the sum of the two stresses caused by the inertial and kinematic effects if T_b<T_g or the square root of the sum of the squares of the two if T_b>T_g. The estimated pile stresses are in good agreement with the observed ones regardless of the occurrence of soil liquefaction.     

The use of MRS in the determination of hydraulic transmissivity: The case of alluvial aquifers

Magnetic Resonance Sounding (MRS) is nowadays accepted as a new geophysical method that can be used for a reliable determination of the ground water content distribution in the top 150 m. A great effort has also been made in MRS development to deduce the hydraulic transmissivity, based on empiric relationships of the permeability with a factor F which is calculated with NMR parameters measured at laboratory scale. To use this relationship under field conditions a calibration coefficient C_T = T_pt / F has to be previously established, which demands the knowledge of the transmissivity T_pt evaluated in the pumping test. The transmissivity can then be calculated at any other site of the same aquifer using the relation T_mrs = C_TF. The C_T values reported suggest a certain relationship with the lithology, but with a great dispersion and contradictory results. MRS surveys carried out in alluvial aquifers in Spain have shown that the value of C_T evaluated at one site may not be valid at another place of the same aquifer, because of the great heterogeneity of this kind of geological environment. The demand of a pumping test at each site where a MRS is measured invalidates the method actually used for MRS transmissivity evaluation. More than 50 MRS have been used to propose a new methodology. The aquifers visited cover a great range of transmissivities (from 2 × 10^− 6 to 9 × 10^− 3 m²/s). The MRS signal amplitude varies between 20 and 1400 nV, the signal/noise ratio is in the range from 0.6 to 42, and the value of the decay time constant varies from 200 to 800 ms. It has been demonstrated that when the transmissivity increases, the value of F decreases, and C_T increases, except for certain groups of MRS taken at the same aquifer or part of one aquifer, for which F increases with T_pt, keeping C_T constant. A function C_T(F) of the type C_T = mF^− n has been obtained that allows the transmissivity evaluation without the need of T_pt. Considering that both values of transmissivity, T_pt and T_mrs, are subjected to deviations due to the experimental errors as well as due to evaluation errors, the prediction achieved by the proposed equation is rather good. To perform a better evaluation of the values of the coefficients m and n it is necessary to have a greater number of MR soundings of good quality and with a trustworthy inversion at locations where a really comparable and good performed pumping test is available, covering a sufficient range of transmissivities. Though the data we have used do not always fulfil these conditions, the result is promising. Once a trustable function is available, the forecast of the transmissivity using MRS will not need the existence of any pumping test in the area. The general extension of this methodology demands the availability of MRS taken at all kinds of geological and hydrogeological environments, which is impossible without the existence of a universal MRS data base.     

Denitrification in a Deep Basalt Aquifer: Implications for Aquifer Storage and Recovery

Aquifer storage and recovery (ASR) can provide a means of storing water for irrigation in agricultural areas where water availability is limited. A concern, however, is that the injected water may lead to a degradation of groundwater quality. In many agricultural areas, nitrate is a limiting factor. In the Umatilla Basin in north central Oregon, shallow alluvial groundwater with elevated nitrate‐nitrogen of <3 mg/L to >9 mg/L is injected into the Columbia River Basalt Group (CRBG), a transmissive confined aquifer(s) with low natural recharge rates. Once recovery of the injected water begins, however, NO₃‐N in the recovered water decreases quickly to <3 mg/L (Eaton et al. 2009), suggesting that NO₃‐N may not persist within the CRBG during ASR storage. In contrast to NO₃‐N, other constituents in the recovered water show little variation, inconsistent with migration or simple mixing as an explanation of the NO₃‐N decrease. Nitrogen isotopic ratios (δ¹⁵N) increase markedly, ranging from +3.5 to > +50, and correlate inversely with NO₃‐N concentrations. This variation occurs in <3 weeks and recovery of <10% of the originally injected volume. TOC is low in the basalt aquifer, averaging <1.5 mg/L, but high in the injected source water, averaging >3.0 mg/L. Similar to nitrate concentrations, TOC drops in the recovered water, consistent with this component contributing to the denitrification of nitrate during storage.     

The instability of the <Emphasis Type="Italic">M</Emphasis><Subscript>max</Subscript> parameter and an alternative to its using

The work presents statistical methods for estimating the distribution parameters of rare, strong earthquakes. Using the two main theorems of extreme value theory (EVT), the distribution of T-maximum (the maximum magnitude over the time period T). Two methods for estimating the parameters of this distribution are proposed using the Generalized Pareto Distribution (GPD) and the General Extreme Value Distribution (GEV). In addition, the that allow the determination of the distribution of the T-maximum for an arbitrary value of T are proposed. The approach being used clarifies the nature of the instability of the widely accepted M max parameter. In the work, instead of unstable values of the M max parameter, the robust parameter Q _T (q), the q level quantile for the distribution of the T-maximum, is proposed to be used. The described method has been applied to the Harvard Catalogue of Seismic Moments of 1977–2006 and to the Magnitude Catalogue for Fennoscandia in 1900–2005. Moreover, the estimates of parameters of the corresponding GPD and GEV distributions, in particular, the most interesting shape parameter and the values of the M _max and Q _T (q) parameters are given.     

A new solution of transient confined–unconfined flow driven by a pumping well

A new solution of transient confined–unconfined flow driven by a pumping well is developed and compared to previous approximate solutions of Moench and Prickett [Moench AF, Prickett TA. Radial flow in an infinite aquifer undergoing conversion from artesian to water table conditions. Water Resour Res 1972;8:494–9] and Hu and Chen [Hu L, Chen C. Analytical methods for transient flow to a well in a confined–unconfined aquifer. Ground Water 2008;46(4):642–6]. The problem is rewritten in dimensionless form with the Boltzmann transform. The nonlinear equation for flow in the unconfined zone is solved with the Runge–Kutta method. Position of the conversion interface is determined with an iteration scheme. This study shows that the confined–unconfined flow depends on three dimensionless parameters that represent the confined–unconfined storativity ratio (a_D), the ratio of the initial hydraulic head over the aquifer thickness (f_i), and the dimensionless pumping rate (q_D). The rate of expansion of the unconfined zone increases with q_D, but decreases with a_D and f_i. Differences between the two previous approximate solutions and the new solution of this study are observable in the estimated position of the conversion interface and the drawdown–time curves. The new solution can be applied to estimate the time for confined–unconfined conversion to occur (critical conversion time), and the time when the pumping well becomes dry (critical drying time). The critical conversion time is found to be very sensitive to the initial hydraulic head. The critical drying time is often much larger than the critical conversion time and may never be observed during a finite pumping period.     

Effect of pore geometry on Gassmann fluid substitution

Although there is no assumption of pore geometry in derivation of Gassmann's equation, the pore geometry is in close relation with hygroscopic water content and pore fluid communication between the micropores and the macropores. The hygroscopic water content in common reservoir rocks is small, and its effect on elastic properties is ignored in the Gassmann theory. However, the volume of hygroscopic water can be significant in shaly rocks or rocks made of fine particles; therefore, its effect on the elastic properties may be important. If the pore fluids in microspores cannot reach pressure equilibrium with the macropore system, assumption of the Gassmann theory is violated. Therefore, due to pore structure complexity, there may be a significant part of the pore fluids that do not satisfy the assumption of the Gassmann theory. We recommend that this part of pore fluids be accounted for within the solid rock frame and effective porosity be used in Gassmann's equation for fluid substitution. Integrated study of ultrasonic laboratory measurement data, petrographic data, mercury injection capillary pressure data, and nuclear magnetic resonance T₂ data confirms rationality of using effective porosity for Gassmann fluid substitution. The effective porosity for Gassmann's equation should be frequency dependent. Knowing the pore geometry, if an empirical correlation between frequency and the threshold pore‐throat radius or nuclear magnetic resonance T₂ could be set up, Gassmann's equation can be applicable to data measured at different frequencies. Without information of the pore geometry, the irreducible water saturation can be used to estimate the effective porosity.     

Sensitivity analysis of a no-crossflow model for the transient flowmeter test

 Logarithmic sensitivities and plausible relative errors are studied in a simple no-crossflow model of a transient flowmeter test (TFMT). This model is identical to the model of a constant-rate pumping test conducted on a fully penetrating well with wellbore storage, surrounded by a thick skin zone, and situated in a homogeneous confined aquifer. The sensitivities of wellbore drawdown and wellface flowrate to aquifer and skin parameters are independent of the pumping rate. However, the plausible relative errors in the aquifer and skin parameters estimated from drawdown and wellface flowrate data can be proportionally decreased by increasing the pumping rate. The plausible relative errors vary by many orders of magnitude from the beginning of the TFMT. The practically important flowrate and drawdown measurements in this test, for which the plausible relative errors vary by less than one order of magnitude from the minimum plausible relative errors, can begin approximately when the dimensionless wellface flowrate exceeds q _D =q/Q≈0.4. During most of this stage of the test, the plausible relative errors in aquifer hydraulic conductivity (K _a ) are generally an order of magnitude smaller than those in aquifer specific storativity. The plausible relative errors in the skin hydraulic conductivity (K _s ) are generally larger than the plausible relative errors in the aquifer specific storativity when the thick skin is normal (K _s >K _a ) and smaller when the thick skin is damaged (K _s <K _a ). The specific storativity of the skin zone would be so biased that one should not even attempt to estimate it from the TFMT. We acknowledge Wiebe H. van der Molen for recommending the De Hoog algorithm and sharing his code. This research was partially supported by the US Geological Survey, USGS Agreement #1434-HQ-96-GR-02689 and North Carolina Water Resources Research Institute, WRRI Project #70165.     

Groundwater investigation on sand dunes area in southern part of Vietnam by magnetic resonance sounding

In the last five years, magnetic resonance sounding (MRS), as a non-invasive geophysical method, has emerged as a new technique for ground water investigation in Vietnam. In this paper, we present the general theoretical basis of this method together with acquisition, processing, and interpretation of the MRS data. We show a case study of MRS surveys in sand dunes area in order to characterize aquifers situated in the southern part of Vietnam. From the interpretation of MRS soundings we delimited an aquifer layer in the subsurface with strong lateral variations for which we determined the depth at 44 m and water content between 3% and 9.5%. The longitudinal relaxation constant T*₁ is about 250 m s, while the transverse relaxation T*₂ is between 150–200 m s. That indicates fine to medium grain size and thus low to medium hydraulic permeability. These results are confirmed by the observations from the well LK1 between 45 to 70 m. The results of other MRS measurements showed the presence of a low water bearing aquifer and were confirmed by the observations in two other wells.     

CO2 storage in the high Arctic: efficient modelling of pre‐stack depth‐migrated seismic sections for survey planning

The sequestration of CO₂ in subsurface reservoirs constitutes an immediate counter‐measure to reduce anthropogenic emissions of CO₂, now recognized by international scientific panels to be the single most critical factor driving the observed global climatic warming. To ensure and verify the safe geological containment of CO₂ underground, monitoring of the CO₂ site is critical. In the high Arctic, environmental considerations are paramount and human impact through, for instance, active seismic surveys, has to be minimized. Efficient seismic modelling is a powerful tool to test the detectability and imaging capability prior to acquisition and thus improve the characterization of CO₂ storage sites, taking both geological setting and seismic acquisition set‐up into account. The unique method presented here avoids the costly generation of large synthetic data sets by employing point spread functions to directly generate pre‐stack depth‐migrated seismic images. We test both a local‐target approach using an analytical filter assuming an average velocity and a full‐field approach accounting for the spatial variability of point spread functions. We assume a hypothetical CO₂ plume emplaced in a sloping aquifer inspired by the conditions found at the University of Svalbard CO₂ lab close to Longyearbyen, Svalbard, Norway, constituting an unconventional reservoir–cap rock system. Using the local‐target approach, we find that even the low‐to‐moderate values of porosity (5%–18%) measured in the reservoir should be sufficient to induce significant change in seismic response when CO₂ is injected. The sensitivity of the seismic response to changes in CO₂ saturation, however, is limited once a relatively low saturation threshold of 5% is exceeded. Depending on the illumination angle provided by the seismic survey, the quality of the images of five hypothetical CO₂ plumes of varying volume differs depending on the steepness of their flanks. When comparing the resolution of two orthogonal 2D surveys to a 3D survey, we discover that the images of the 2D surveys contain significant artefacts, the CO₂‐brine contact is misplaced and an additional reflector is introduced due to the projection of the point spread function of the unresolvable plane onto the imaging plane. All of these could easily lead to a misinterpretation of the behaviour of the injected CO₂. Our workflow allows for testing the influence of geological heterogeneities in the target aquifer (igneous intrusions, faults, pervasive fracture networks) by utilizing increasingly complex and more realistic geological models as input as more information on the subsurface becomes available.     

Binary Mixtures of Permanganate and Chlorinated Volatile Organic Compounds in Groundwater Samples: Sample Preservation and Analysis

Groundwater samples collected at sites where in situ chemical oxidation (ISCO) has been deployed may contain binary mixtures of groundwater contaminants and permanganate (MnO₄^–), an oxidant injected into the subsurface to destroy the contaminant. Commingling of the oxidant and contaminant in aqueous samples may negatively impact the quality of the sample as well as the analytical instruments used to quantify contaminant concentrations. In this study, binary mixtures comprised of (1) a multicomponent standard with permanganate and (2) groundwater samples collected at two ISCO field sites were preserved with ascorbic acid. Ascorbic acid reacts rapidly with the MnO₄^– and limits the reaction between MnO₄^– and the organic compounds in the mixture. Consequently, most of the compounds in the multicomponent standard were within the control limit for quality assurance. However, despite timely efforts to preserve the samples, the rapid reaction between permanganate and contaminant caused the concentration of several sensitive compounds to fall significantly below the lower control limit. Concentrations of volatile organic compounds in the field‐preserved binary mixture groundwater samples were greater than in samples refrigerated in the field and preserved upon arrival at the laboratory, indicating the time‐dependency and benefit of field preservation. The molar ratio of ascorbic acid required to neutralize KMnO₄ was 1.64 (mol ascorbic acid/mol KMnO₄); this provided a baseline to estimate the volume of ascorbic acid stock solution and/or the weight of crystalline ascorbic acid required to neutralize MnO₄^–. Excess ascorbic acid did not negatively impact the quality of the aqueous samples, or analytical instruments, used in the analyses.     

Inelastic displacement demands in steel structures and their relationship with earthquake frequency content parameters

This paper deals with the estimation of peak inelastic displacements of SDOF systems, representative of typical steel structures, under constant relative strength scenarios. Mean inelastic deformation demands on bilinear systems (simulating moment resisting frames) are considered as the basis for comparative purposes. Additional SDOF models representing partially‐restrained and concentrically‐braced (CB) frames are introduced and employed to assess the influence of different force‐displacement relationships on peak inelastic displacement ratios. The studies presented in this paper illustrate that the ratio between the overall yield strength and the strength during pinching intervals is the main factor governing the inelastic deformations of partially‐restrained models and leading to significant differences when compared with predictions based on bilinear structures, especially in the short‐period range. It is also shown that the response of CB systems can differ significantly from other pinching models when subjected to low or moderate levels of seismic demand, highlighting the necessity of employing dedicated models for studying the response of CB structures. Particular attention is also given to the influence of a number of scalar parameters that characterise the frequency content of the ground motion on the estimated peak displacement ratios. The relative merits of using the average spectral period T_aver, mean period T_m, predominant period T_g, characteristic period T_c and smoothed spectral predominant period T_o of the earthquake ground motion, are assessed. This paper demonstrates that the predominant period, defined as the period at which the input energy is maximum throughout the period range, is the most suitable frequency content scalar parameter for reducing the variability in displacement estimations. Finally, noniterative equivalent linearisation expressions based on the secant period and equivalent damping ratios are presented and verified for the prediction of peak deformation demands in steel structures. Copyright © 2011 John Wiley & Sons, Ltd.     

In situ Mixing of Organic Matter Decreases Hydraulic Conductivity of Denitrification Walls in Sand Aquifers

In a previous study, a denitrification wall was constructed in a sand aquifer using sawdust as the carbon substrate. Ground water bypassed around this sawdust wall due to reduced hydraulic conductivity. We investigated potential reasons for this by testing two new walls and conducting laboratory studies. The first wall was constructed by mixing aquifer material in situ without substrate addition to investigate the effects of the construction technique (mixed wall). A second, biochip wall, was constructed using coarse wood chips to determine the effect of size of the particles in the amendment on hydraulic conductivity. The aquifer hydraulic conductivity was 35.4 m/d, while in the mixed wall it was 2.8 m/d and in the biochip wall 3.4 m/d. This indicated that the mixing of the aquifer sands below the water table allowed the particles to re-sort themselves into a matrix with a significantly lower hydraulic conductivity than the process that originally formed the aquifer. The addition of a coarser substrate in the biochip wall significantly increased total porosity and decreased bulk density, but hydraulic conductivity remained low compared to the aquifer. Laboratory cores of aquifer sand mixed under dry and wet conditions mimicked the reduction in hydraulic conductivity observed in the field within the mixed wall. The addition of sawdust to the laboratory cores resulted in a significantly higher hydraulic conductivity when mixed dry compared to cores mixed wet. This reduction in the hydraulic conductivity of the sand/sawdust cores mixed under saturated conditions repeated what occurred in the field in the original sawdust wall. This indicated that laboratory investigations can be a useful tool to highlight potential reductions in field hydraulic conductivities that may occur when differing materials are mixed under field conditions.     

DETERMINATION OF THE EFFECTIVE POROSITY FOR PARTIALLY PENETRATED AQUIFERS

Abstract

The effective porosity θ _e for partially penetrated aquifers was determined. The model basin sandy aquifer available in the Centre was used. The values obtained for θ _e were in good agreement with the adopted values.     

Accuracy in measurements and the temperature and volume dependence of thermoelastic parameters

To obtain the temperatureT and volumeV (or pressureP) dependence of the Anderson-Grüneisen parameter _T , measurements with high sensitivity are required. We show two examples:P, V, T measurements of NaCl done with the piston cylinder and elasticity measurements of MgO using a resonance method. In both cases, the sensitivity of the measurements leads to results that provide information about _T (,T), where V/V ₀ andV ₀ is the volume at zero pressure. We demonstrate that determination of _T leads to understanding of the volume and temperature dependence ofq=( ln / lnV) _T over a broadV, T range, where is the Grüneisen ratio.     

Solute dynamics across the stream‐to‐riparian continuum under different flood waves

To enhance the understanding of solute dynamics within the stream‐to‐riparian continuum during flood event‐driven water fluctuation (i.e., flood wave), a variable saturated groundwater flow and solute transport model were developed and calibrated against in situ measurements of the Inbuk stream, Korea, where seasonal flooding prevails. The solute dynamics were further investigated for flood waves (varying by amplitude [A], duration [T], roundness [r], and skewness [t_p]) that were parameterised by real‐time stream stage fluctuations. We found that the solute transferred faster and farther in the riparian zone, especially within the phreatic zone, above which in the variable saturated zone the concentration required a significantly longer time, particularly at higher altitudes, to return to the initial state. By comparison, solute transferred shallowly in the streambed where the solute plume exhibited an exponential growth trend from the centre to the bank. The dynamic changes of solute flux and mass along the stream–aquifer interface and stream concentration were linked to the shape of flood wave. As the flood wave became higher (A↗), wider (T↗), rounder (r↘), and less skewed (t_p↗), the maximum solute storage in aquifer increased. Maximum stream concentration (C_str?max) not only presented a positive linear relationship with A or t_p but also showed a negative logarithmic trend with increasing T or r. The sensitivity of C_{str_max} to A was approximately two times that of t_p, and between these values, the r was slightly more sensitive than T. C_str?max linearly increased as hydraulic conductivity increased and logarithmically increased as longitudinal dispersivity increased. The former relationship was more sensitive than the latter.     

Transient Divergent Flow and Transport in an Infinite Anisotropic Porous Formation

When seeking to predict plume geometry resulting from fluid injection through partially penetrating wells, it is common to assume a steady-state spherically diverging flow field. In reality, the flow field is transient. The steady-flow assumption is likely to cause overestimation of injection plume radius since the accommodation of fluid by increases in porosity and fluid density is ignored. In this paper, a transient solution is developed, resulting in a nonlinear ordinary differential equation expressing plume radius as a function of time. It is shown that the problem can be fully described by one type curve. A critical time, t_c, is identified at which the percentage error of the steady-state flow solution compared to the fully dynamic problem is less than 1%. Only for large injection rates and low permeabilities, does t_c become greater than 1 h. Nevertheless, an improved approximate solution is obtained by a simple linearization procedure. The critical time, t_c for the new approximate solution is 0.3% of that required for the steady-state flow solution.     

The pore structure of tight limestone—Jurassic Ziliujing Formation,Central Sichuan Basin,China

The pore structure of the tight limestone in the Daanzhai Member of the Ziliujing Formation, Jurassic System, in central Sichuan Basin, China, is complex but essential to the exploration and development of tight oil. The pore structure of the tight limestone is studied by using scanning electron microscopy (SEM), nitrogen adsorption, high-pressure mercury intrusion, and nuclear magnetic resonance (NMR). The experimental results suggest that the pores are mainly slit pores and mesopores and macropores contribute to the pore volume and specific surface. The displacement pressure, average pore size, and homogeneity coefficient correlate with porosity and permeability and can be used to evaluate the pore structure. The full pore-size distribution was obtained by combining nitrogen adsorption and high-pressure mercury intrusion. We find that the limestone mainly contains mesopores with diameter of 2?50 nm. The T₂ distribution was converted into pore-size distribution, well matching the full pore-size distribution. The relation between T₂ and pore size obeys a power law and the geometric mean of T₂ correlates with the pore structure and can be used in the pore structure evaluation.