致密砂岩的形成机制及其地质意义——以塔里木盆地英南2井为例

砂质岩的孔隙和喉道被网格状粘土矿物和次生加大矿物充填成微细孔喉状结构时,形成具有较高毛细管压力的致密砂岩。它与地层水发生水锁效应,可大大降低渗透率,成为致密砂岩盖层。当含水饱和度在50%以下时,束缚水饱和度比较低,致密砂岩储层可以产气;当含水饱和度在50%～90%区间时,具有较高的束缚水饱和度,相对渗透率非常低,它既不能产气也不能产水,反映为渗透率瓶颈区(具有盖层性质);当含水饱和度大于90%以上时,致密砂岩储层仅微量产水。塔里木盆地英南2井侏罗系气藏盖层由致密砂岩构成,不含水时的气体渗透率在(0.027～0.081)×10-3μm2,不能构成封堵;当含水饱和度达到60%以上时,相对渗透率几乎为零,构成有效盖层。     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

Zircon in amphibolites from Naxos,Aegean Sea,Greece: origin,significance and tectonic setting

We report U–Pb zircon ages of c. 700–550 Ma, 262–220 Ma, 47–38 Ma and 15–14 Ma from amphibolites on Naxos Island in the Aegean extensional province of Greece. The zircon has complex internal structures. Based on cathodoluminescence response, zoning and crosscutting relationships a minimum of four zircon growth stages are identified: inherited core, magmatic core, inner metamorphic (?) rim and an outer metamorphic rim. Trace element compositions of the amphibolites suggest igneous differentiation and crustal assimilation. Zircon solubility as a function of saturation temperatures, Zr content and melt composition indicates that the zircon did not originally crystallize in the mafic bodies but was inherited from felsic precursor rocks, and subsequently assimilated into the mafic intrusives during emplacement. Zircon inheritance is corroborated by the complex, xenocrystic nature of the zircon in one sample. Ages of c. 700–550 Ma and 262–220 Ma are assigned to inherited zircon. Available geochemical data suggest that the 15–14 Ma metamorphic rims grew in situ in the amphibolites, corresponding to a high‐grade metamorphic event at this time. However, the geochemical data cannot conclusively establish if the c. 40 Ma zircon rims also grew in situ, or whether they were inherited along with the xenocrystic cores. Two scenarios for emplacement of the mafic intrusives are discussed: (i) Intrusion during late‐Triassic to Jurassic ocean basin development of the Aegean realm, in which case the 40 Ma zircon rims would have grown in situ, and (ii) emplacement in the Miocene as a result mafic underplating during large‐scale extension. In this case, only the 15–14 Ma metamorphic outer rims would have formed in situ in the amphibolitic host rocks.     

Alteration of hydrological processes and streamflow with juniper (Juniperus virginiana) encroachment in a mesic grassland catchment

To understand the effect of woody plant encroachment on hydrological processes of mesic grasslands, we quantified infiltration capacity in situ, the temporal changes in soil water storage, and streamflow of a grassland catchment and a catchment heavily encroached by juniper (Juniperus virginiana, eastern redcedar) in previously cultivated, non‐karst substrate grasslands in north‐central Oklahoma for 3 years. The initial and steady‐state infiltration rates under the juniper canopy were nearly triple to that of the grassland catchment and were intermediate in the intercanopy spaces within the encroached catchment. Soil water content and soil water storage on the encroached catchment were generally lower than on the grassland catchment, especially when preceding the seasons of peak rainfall in spring and fall. Frequency and magnitude of streamflow events were reduced in the encroached catchment. Annual runoff coefficients for the encroached catchment averaged 2.1%, in contrast to 10.6% for the grassland catchment. Annual streamflow duration ranged from 80 to 250 h for the encroached catchment compared with 600 to 800 h for the grassland catchment. Our results showed that the encroachment of juniper into previously cultivated mesic grasslands fundamentally alters catchment hydrological function. Rapid transformation of mesic grassland to a woodland state with juniper encroachment, if not confined, has the potential to drastically reduce soil water, streamflow and flow duration of ephemeral streams in the Southern Great Plains. Copyright © 2013 John Wiley & Sons, Ltd.     

Micromechanical approach to swelling behavior of capillary-porous media with coupled physics

A detailed multiscale analysis is presented of the swelling phenomenon in unsaturated clay-rich materials in the linear regime through homogenization. Herein, the structural complexity of the material is formulated as a three-scale, triple porosity medium within which microstructural information is transmitted across the various scales, leading ultimately to an enriched stress-deformation relation at the macroscopic scale. As a side note, such derived relationship leads to a tensorial stress partitioning that is reminiscent of a Terzaghi-like effective stress measure. Otherwise, a major result that stands out from previous works is the explicit expression of swelling stress and capillary stress in terms of micromechanical interactions at the very fine scale down to the clay platelet level, along with capillary stress emerging due to interactions between fluid phases at the different scales, including surface tension, pore size, and morphology. More importantly, the swelling stress is correlated with the disjoining forces due to electrochemical effects of charged ions on clay minerals and van der Waals forces at the nanoscale. The resulting analytical expressions also elucidate the role of the various physics in the deformational behavior of clayey material. Finally, the capability of the proposed formulation in capturing salient behaviors of unsaturated expansive clays is illustrated through some numerical examples.     

Constitutive modelling of fine-grained gassy soil: A composite approach

Fine-grained marine sediments containing large undissolved gas bubbles are widely distributed around the world. Presence of the bubbles could degrade the undrained shear strength (s_u ) of the soil, when the gas pressure u_g is relatively high as compared with the effective stress in the saturated soil matrix. Meanwhile, the addition of bubbles may also increase s_u when the difference between u_g and pore water pressure u_w becomes smaller than the water entry value, causing partial water drainage from the saturated matrix into the bubbles (bubble flooding) during globally undrained shearing. A new constitutive model for describing the two competing effects on the stress-strain relationship of fine-grained gassy soil is proposed within the framework of critical state soil mechanics. The gassy soil is considered as a three-phase composite material with compressible cavities, which allows water entry from the saturated matrix. Bubble flooding is modelled by introducing an additional positive volumetric strain increment of the saturated clay matrix, which is dependent on the difference between pore gas and pore water pressure based on experimental observations. A modified hardening law based on that of the modified Cam clay model is employed, which in conjunction with the expression for bubble flooding, can describe both the detrimental and beneficial effects of gas bubbles on soil strength and plastic hardening in shear. Only two extra parameters in addition to those in the modified Cam clay model are used. It is shown that the key features of the stress-strain relationship of three fine-grained gassy soils can be reproduced satisfactorily.     

便携式Li-K分析仪的研制及其在锂辉石中锂的分析应用

我国锂钾主要矿产资源大多分布在西部偏远地区,其勘查找矿或综合利用迫切需要现场快速分析技术的支持。本文介绍了自行开发的便携式Li-K分析仪的主要性能及其在锂辉石中锂的分析应用。Li-K分析仪是一种基于大气压的液体阴极辉光放电光谱仪,以待测液体为放电阴极,实现了样品中Li、K等元素的原子化和激发。该仪器以光纤传导CCD光谱仪作为检测器,其波长范围为345~1015 nm,分辨率3 nm;以Li 670.78 nm和K 769.90 nm分析谱线,在仪器最佳工作条件下测定Li、K的精密度(RSD,n=14)均低于2%,检出限为0.03μg/m L,检测范围0.1~10μg/m L。研究表明,不同酸度和酸的类型对谱线强度影响较大,且存在显著的样品基体效应,标准曲线法分析锂辉石中Li的结果偏差高达267%;而标准加入法可克服基体效应的影响,获得与ICP-OES一致性较好的分析结果,为现场开展固体样品中Li的测定奠定了基础。     

伊通盆地鹿乡断陷C区块含油饱和度与含油性分析

综合研究伊通盆地鹿乡断陷双二段储层的岩性、含油性和含油饱和度,确定了适用于本区的含油饱和度计算方法,并与实际含油性进行了对比,分析了沉积和构造对含油性的影响。研究区双二段含油储层的岩性以粉砂岩、细砂岩和砂砾岩为主,含油级别为荧光级以上。根据阿尔奇公式计算出的含油饱和度大部分与实际的含油性相符合,并将含油饱和度70%、50%和38%值作为本区气层、油层、油水层和水层的区分标准。本区储层的含油性受到沉积和构造两方面的控制,沉积作用主要控制储层的岩石结构和物性,构造控制油气的聚集,即在构造高部位,实际含油饱和度比计算结果要高,而在构造低部位,实际含油饱和度比计算结果要低。     

高效液相色谱-离子阱质谱法测定海洋微藻藻粉中的8种脂溶性毒素

建立了利用高效液相色谱-电喷雾离子阱质谱(HPLC-ESI-IT-MS)测定海洋微藻藻粉中8种典型脂溶性毒素的分析方法。藻粉样品经超声细胞破碎后,采用超声波辅助提取法对藻毒素进行提取,用HPLC-ESI-IT-MS多反应离子监测(MRM)模式对各种毒素(包括大田软海绵酸(OA)、鳍藻毒素1(DTX-1)、扇贝毒素2(PTX-2)、虾夷扇贝毒素(YTX)、原多甲藻酸1(AZA1)原多甲藻酸2(AZA2)、罗环内酯毒素(SPX),米氏裸甲藻毒素(GYM))进行测定。8种脂溶性藻毒素均在线性范围内线性关系良好(R2均在0.991以上),检出限均介于0.085～1.315 pg之间,加标回收率在88.5%～111.4%之间,方法重复性相对标准偏差(RSD)在4.82%～10.17%范围。应用该方法对利玛原甲藻干藻粉中的毒素进行了测定,分析结果良好,说明本方法是海洋微藻藻粉中脂溶性藻毒素测定的有效方法。