Experimental study of hydrothermal alteration of carbonate rocks by Na-F solutions under flow conditions

Na-F solutions with different pH were allowed to react under flow conditions with limestone, dolomite and siderite at 300 and 400°C at 500 bars, in order to understand the geochemical characteristics of the interaction between rock and F-bearing hydrothermal solutions. The mineralogical and chemical composition of these carbonate rocks was studied by chemical analysis, XRD and SEM/EDAX before and after reaction with Na-F solutions. Based on an approximate thermodynamic calculation the formation of alteration phases is discussed. The pH value at room conditions and the contents of K, Na, Ca, Mg, Fe, F, SiO₂ and total dissolved CO₂ in solution after interaction with carbonate rocks were determined. From solution chemistry data, it is found that the amount of Na and F in the reacted solutions seemed to have roughly similar fluctuations with time, and hence, carbonate might be an effective precipitating agent for Fe in hydrothermal solutions. The present experimental study demonstrates that hydrothermal metasomatism played a very important role in the formation of the Baiyun Ebo iron deposit.     

钾长石在碱性条件下的蚀变机制及其蚀变产物托贝莫来石的纳米结构研究

钾长石在碱性流体蚀变过程中会形成层状铝硅酸盐等矿物,前人对其反应机制和反应产物进行了研究,但缺乏微观尺度尤其是纳米尺度的探讨。因此,作者对钾长石在极端碱性条件下(190℃,24h,初始p H=12.4)的蚀变机制及其蚀变产物层状硅酸盐托贝莫来石的显微结构开展了纳米尺度的研究。X射线粉末衍射、扫描电子显微镜、能量色散光谱等观测结果显示,钾长石在碱性条件下水热蚀变所得到的产物主要为托贝莫来石、水钙铝榴石和方解石。高分辨率的透射电镜结果表明,在钾长石与次生矿物相的界面形成了纳米级的多孔非晶质层,且在空间上表现出结构的不连续性。界面溶解-再沉淀(CIDR)机制很好地解释了钾长石与次生矿物相界面的空间不连续性和非晶质层的形成。对蚀变产物中纤维状的托贝莫来石晶体进行显微结构表征,结果表明托贝莫来石的孔隙直径为0～160nm,平均孔径约为40nm;其构成的纳米孔隙和通道有利于增加周围流体中离子和气体的溶解度,并可能会影响局部化学平衡。这为层状铝硅酸盐作为自然界以及工业吸附材料和催化剂的更广泛应用提供了重要依据。     

渗流-化学溶解耦合作用下岩石单裂隙渗透特性研究

为揭示在渗流-化学溶解耦合作用下单裂隙渗透特性的变化规律,建立了描述二维渗流-化学溶解耦合作用的偏微分方程组,并利用COMSOL Multiphysics软件成功地求解该方程组。首先,模拟了文献[1]中的盐岩渗流-溶解耦合渗流试验结果,数值模拟结果与试验结果较为吻合,验证了数学模型的正确性和有效性。然后,利用分形理论生成了一个粗糙的裂隙面数字模型,着重分析了二维石灰岩粗糙裂隙面在水流、矿物溶解和输运过程中其渗透特性的变化规律。数值分析显示,（1）溶质浓度对裂隙面的溶解具有非常重要的作用,从而水流进口端的溶解厚度比出口端大得多。（2）裂隙的整体渗透性在初始时刻增加较慢,随着裂隙开度的增大和贯通,溶解速度会逐渐增大,是一个加速的过程。     

Hydrothermal alteration of basalt by seawater under seawater-dominated conditions

Fresh mid-ocean ridge basalt glass and diabase have been reacted with seawater at 150–300°C, 500 bar, and water/rock mass ratios of 50, 62, and 125, using experimental apparatus which allowed on-line sampling of solution to monitor reaction progress. These experiments characterize reaction under what we have called “seawater-dominated” conditions of hydrothermal alteration.In an experiment at 300°C, basalt glass undergoing alteration removed nearly all Mg²⁺ from an amount of seawater 50 times its own mass. In the process, the glass was converted entirely to mixed-layer smectite-chlorite, anhydrite, and minor hematite. Removal of Mg from seawater occurred as a Mg(OH)₂ component incorporated into the secondary clay. This produced a precipitous drop in solution pH early in the experiment, accompanied by a dramatic increase in the concentrations of Fe, Mn, and Zn in solution. As Mg removal neared completion and the glass was hydrolyzed, pH rose again and heavy metal concentrations dropped.At water/rock ratios of 62 and 125 and 150–300°C, the mineral assemblage produced was similar to that at a water/rock ratio of 50. Solution chemistry, however, contrasted with the earlier experiment in that Mg concentrations in solution were greater and pH lower. This caused significant leaching of heavy metals. At 300°C nearly all of the Na, Ca, Cu, Zn, and CO₂ and most of the K, Ba, Sr, and Mn were leached from the silicates. H₂S, Al, Si, and possibly Co were also significantly mobilized, whereas V, Cr, and Ni were not. Little or no seawater sulfate was reduced.Although submarine hot spring solutions sampled to date along mid-ocean ridges clearly come from rock-dominated hydrothermal systems, evidence from ocean floor metabasalts and from heat flow studies indicates that seawater-dominated conditions of alteration prevail at least locally both in axial hightemperature systems and in ridge flank systems at lower temperatures.     

人工回灌条件下滨海砂质含水层渗透性的时空演变特征

人工回灌补给滨海含水层是海水入侵修复经常采用的工程措施,然而人工回灌条件下滨海砂质含水层经常出现渗透性显著降低的现象,其含水介质渗透性的时空演变规律仍不清楚。本文以青岛市大沽河下游咸水入侵区含水层砂样为研究对象,通过利用不同尺度的室内砂柱淡咸水驱替试验,对定水头和定流速条件下人工回灌咸淡水驱替过程中砂质含水层渗透性的时空演变特征进行了研究。结果表明：当淡水驱替咸水时含水介质的渗透性会发生显著变化,整个砂柱的渗透性先降低后回升,原因是回灌过程中砂柱中的黏土矿物发生释放、迁移、沉积;定水头条件下,在砂柱的前半段存在一个淘空区,渗透系数是初始渗透系数的1.6~2.0倍,砂柱其它部分渗透性是先降低再轻微回弹;定流速条件下,砂柱的渗透系数随时间的变化均呈现先降低后升高趋势,且砂柱各段离入水口的距离越远其渗透系数的值越小;砂柱黏粒含量的变化规律与渗透性的演化特征相吻合。研究成果可为人工回灌治理滨海含水层海水入侵提供一定的科学依据。     

Mineralogical alteration and associated permeability changes induced by a high-pH plume: Modeling of a granite core infiltration experiment

Within the framework of the HPF project (Hyperalkaline Plume in Fractured Rock) at the Grimsel Test Site (Switzerland), a small scale core infiltration experiment was performed at the University of Bern. A high-pH solution was continuously injected, under a constant pressure gradient, into a cylindrical core of granite containing a fracture. This high-pH solution was a synthetic version of solutions characteristic of early stages in the degradation of cement. The interaction between the rock and the solutions was reflected by significant changes in the composition of the injected solution and a decrease in the permeability of the rock. Changes in the mineralogy and porosity of the fault gouge filling the fracture were only minor. One-dimensional reactive transport modeling, using a modified version of the GIMRT code, was used to interpret the results of the experiment. Dispersive and advective solute transport, mineral reaction kinetics and a coupling between porosity and permeability changes were taken into account. In order to obtain a reasonable agreement between models and experimental results, reactive surface areas of the order of 10⁵ m²/m³ rock had to be used. These values are much smaller than the values measured for the fault gouge filling the fracture, which are in the order of 10⁶–10⁷ m²/m³ rock. However, the results could be improved by adding a small fraction of fine grained mineral, which could explain the high initial peaks in Al and Si concentration. With the inclusion of this fine grained fraction, the initial surface areas in the model were within the range of the measured specific surface areas of the fault gouge. The fact that the decrease in permeability was significant despite the minor changes in mineralogy, suggests that permeability may be controlled by changes in the structure of the rock (pore geometries) rather than by only the bulk volumetric (porosity) changes.     

Effect of temperature on the dissolution and thermal alteration of combined amino acids fixed in natural sediment under simulated hydrothermal conditions

Seafloor sediment containing biogenic amino acids was heated with NaCl solutions at 50–200 °C for 240 h to investigate the dissolution process of amino acids and evaluate their stabilities under hydrothermal conditions. Dissolved amino acids in the combined phase (dissolved combined amino acids, DCAAs) and free phase (dissolved free amino acids, DFAAs) were rapidly released into the solution during heating. The amount of DCAAs in the solutions was 4–9 times higher than the amount of DFAAs at each temperature. When heated at ⩽ 100 °C, most of the total dissolved hydrolyzable amino acids (TDHAAs) were in the combined form (DCAAs/TDHAAs ratios > 0.9). The compositions of the DCAAs in solutions heated at ⩽ 100 °C were similar to that of the total hydrolyzable amino acids (THAAs) of the initial sediment, indicating that the DCAAs, which are derived from organisms and biodebris in the sediment, are barely altered during the hydrothermal reaction at these temperatures. On the other hand, the DCAAs/TDHAAs ratios were 0.72 and 0.57 at 150 and 200 °C, respectively, and the compositions of the DCAAs at 150 and 200 °C were significantly different from that of the initial THAAs. In addition, non-protein amino acids (β-alanine and γ-aminobutyric acid), which are sensitive biochemical indicators of the diagenetic alteration of natural organic matter, drastically increased to 80.9% of the DCAAs after heating at 200 °C. These results suggests that DCAAs are thermally unstable in the hydrothermal solutions at ⩾ 150 °C. These DCAA would be transformed into thermally stable geo-polymers such as humic-like substances and hydrolyzable kerogens.     

Mineralogical and geochemical characteristics of the hydrothermal illite from Hoam granite,South Korea: Implications for episodic fluid injections in the hydrothermal alteration system

Two stages of illite mineralization are recognized in the hydrothermal alteration zone of the Hoam granite. These illites are formed as a result of pervasive alteration by re-equilibration with high water/rock in a brittle environment below <2 km; the mineralization timing is middle Oligocene (26–27 Ma), coinciding with the timing of crustal deformation related to the opening of the East Sea (Sea of Japan). The mineralogical and geochemical characteristics of the clearly distinguished illites at each site indicate that they were mineralized from different fluid injections in distinct geological environments. Illites at the site-1 alteration zone are characterized by high-K content [K_0.84 per O₁₀(OH)₂], 2M₁ polytype of 99 %, hexagonal plate shape, and coexistence with pyrite. These observations indicate that the illites were formed in a slow cooling system (>250 °C), high fS₂, and a relatively acidic environment. The pseudomorphic replacement combined with matrix-filling texture indicates that the illites at the site-1 alteration zone recorded the changes in fluid conditions from low to high water-rock ratio. In contrast, the illites at the site-2 alteration zone show the coexistence of polytypes (2M₁, 1M, and 1M_d), high-K illites [(K_0.83 per O₁₀(OH)₂]/low-K illites [K_0.63 per O₁₀(OH)₂], platy/hairy shapes, and presence of magnetite. Furthermore, this alteration zone no longer exhibits primary textures because of pervasive alteration induced by the dissolution-precipitation process. These results indicate that they were formed in a rapid cooling system and were continuously under conditions of high water-rock ratio, as well as in a less acidic and fS₂ environment than that observed at the site-1 alteration zone. The behavior of trace elements for each illite primarily depends on the constituents of the hydrothermal fluid, which reflect different degrees of fluid evolution. The enrichment of high field strength elements (Nb and Ta), large ion lithophile elements (B, Be, and Cs), rare earth elements, and actinide elements (U and Th) in illite at the site-2 alteration zone shows that these elements formed by a more evolved fluid than that of the illite at the site-1 alteration zone. In addition, negative Ce anomalies at the site-2 alteration zone indicate that these crystallized in a reducing environment. Considering the mineralogical and geochemical properties of illites at the site-1 and site-2 alteration zones, the illite mineralization in the Hoam granite was likely generated by at least two episodes of hydrothermal illite mineralization, which originated from episodic injections of fluids, rather than continuously evolved fluids.     

西藏羊八井地热田水热蚀变的时空演化特征

西藏亚东—谷露裂谷中—北段的羊八井地热田是我国著名的高温地热田,研究其水热蚀变的时空演化有助于更好地认识藏南地热的发育特征。通过对羊八井地热田及其水热蚀变岩开展地表调查、显微特征与X衍射分析等工作,总结了其主要蚀变类型特征,划分出黄褐色蚀变中心带、灰白色中强蚀变带、灰白色中等蚀变带和浅灰白色弱蚀变带4个不同的水热蚀变带,并区分出红褐色—黄褐色蚀变期、灰白色蚀变期和淡黄色—灰色蚀变期共3期蚀变。研究结果揭示,羊八井地热田高温地热活动中心一直在北区硫磺沟区域,其水热蚀变活动主要受亚东—谷露裂谷内部的活动断裂构造控制,并与断裂构造活动具同步性;地热水的排泄方式早期为沿北东向断裂构造直接排泄,晚期为经浅层第四系径流后再排泄,由直接排泄向间接排泄转变;中高温地热水的排泄区由北区硫磺沟地区向南区藏布曲迁移。根据研究结果推断,硫磺沟区域的北东向断裂与北西向断裂交汇区可作为羊八井热田北区深部地热勘查的主要方向。     

Hydromechanical properties and alteration of natural fracture surfaces in the Soultz granite (Bas-Rhin, France)

Natural fractures are characterized by rough surfaces and complex fluid flows. A large distribution of apertures (residual voids) within their walls and the presence of contact points (in situ normal loads) produce heterogeneous flows (channeling). The resulting permeabilities, porosities or fluid–rock exchange surfaces cannot be realistically modeled by parallel and smooth plate models. Four natural fractures are sampled at different depths and degrees of alteration in the Soultz sandstone and granite (EPS1 drillhole, Soultz-sous-Forêts, Bas-Rhin, France). The fracture surfaces are measured with mechanical profilometry and maps of asperity heights (XYZ). Resulting local apertures (XYe) are then calculated. A statistical study of the surface profiles (XZ) show that the fractures are more or less rough and tortuous according to the types of alteration. Altered samples are characterized by smoother surfaces of fractures. Such differences imply that (i) the average fracture aperture is not representative for the whole fracture and that (ii) the different local apertures should be integrated in hydraulic and mechanical models. A hydraulic model (finite difference calculations) of fluid flow, taking into account the elastic closure (Hertz contact theory) of fractures with depth, is used. Maps of contact points and relative local loads within the fracture planes are compared to flow maps. They show different channeling of fluid flows. Strongly altered fractures are characterized by homogeneous fluxes despite the presence of numerous contact zones during the closure of fracture. By contrast, fresh fractures develop, increasing fluid flow channels with depth.Fracture closure (increasing normal stress) does not systematically increase the channeling of fluid flow. There is evidence for a general smoothing out of the irregularities of the fracture walls due to precipitation of secondary minerals, indicating that the cubic law can be commonly valid, also at great crustal depth but this validity depends on the degree of fracture alteration. Mineralogical and geochemical observations, thus, should be taken into account to perform more accurate permeability calculations and models of fluid circulation in fracture networks.     

The reaction kinetics of alanine and glycine under hydrothermal conditions

Experimental data on the hydrothermal reaction kinetics of α-alanine, glycine, and β-alanine were acquired using a custom-built spectrophotometric reaction cell which permits in situ observation under hydrothermal conditions. Quantitative kinetic information, including rate constants, concentration versus time profiles, and calculations of the individual component spectra, was obtained from the data using a self-modeling chemometric approach based on factor analysis which treats the rate expressions simultaneously as a system of differential algebraic equations (DAE) of index 1. Experimental data collected at 120-165 °C and 20 bar indicates that aqueous α-alanine, glycine and β-alanine will preferentially undergo dimerization and subsequent cyclization when heated in an inert reactor. The results presented here lend further support to the roles of temperature, exposed reactive surfaces, and matrix additives in the reaction kinetics of the structurally simple amino acids examined in this study.     

Origin and granite alteration effects of hydrothermal fluid: isotopic evidence from fluorite veins, Co. Galway, Ireland

The Costelloe Murvey Granite is a chemically evolved, high heat production, leucocratic component of the 400 Ma old Galway Granite batholith and is host to hydrothermal fluorite-quartz-calcite veins. A previously reported clinopyroxene ⁴⁰Ar-³⁹Ar age of 231±4 Ma obtained from a pre-mineralization dolerite dyke is reinterpreted as dating this mineralization. The hydrothermal fluid extensively altered its granite wallrocks, leading to lower Sm and Nd and higher Rb concentrations in altered granite, disturbing both its Rb-Sr and Sm-Nd isotopic systems. The ⁸⁷Sr/⁸⁶Sr ratio of the hydrothermal fluid from which fluorite and calcite precipitated ranged from 0.7101 to 0.7139. These ratios are very much lower than in the Costelloe Murvey Granite at the time of mineralization, precluding the granite as a source for more than 2% of the hydrothermal Sr. The initial ¹⁴³Nd/¹⁴⁴Nd ratio varies between fluorite in different veins due to Nd derivation from local wallrocks, and between fluorite of petrographically distinct growth phases within a single hand specimen, highlighting the difficulty of Sm-Nd isochron dating of fluorite in cases where there are multiple sources of hydrothermal Nd. It is proposed that fluorite and calcite precipitated where hot, dilute fluids rising through the granite mixed with cooler, more saline fluids of basinal origin migrating through Lower Carboniferous limestone which then overlay the granite. Received: 3 August 1995 / Accepted: 11 April 1996     

压剪作用下裂隙扩展过程渗流与应力耦合模型研究

裂隙的剪胀特性及扩展演化规律对岩体的渗流特性具有重要影响。为揭示裂隙剪胀及扩展演化对岩体渗流的影响,基于残余强度提出了一种能较好描述岩石硬化-软化特性的全剪切本构关系;结合剪切变形与裂隙开度的关系,利用最小势能原理和立方定律,建立了压剪作用下考虑裂隙剪胀特性的渗流应力耦合模型;假定压剪作用下裂隙发生Ⅰ型扩展,提出了伴有翼型裂纹的渗流模型,该模型不仅考虑了岩石的剪胀特性,更反映了裂隙扩展过程渗流的演化规律。对不同裂隙粗糙度的剪切应力-位移曲线进行分析,全剪切本构模型表现出较高的拟合精度。在剪切应力-位移关系基础上,通过剪切渗流试验数据对压剪作用下渗流模型进行验证,结果表明,该模型能较好地描述岩体剪胀阶段渗透性演化规律。利用等效裂隙简化裂隙网络,并通过试验数据进行验证,证实了裂隙扩展过程渗流-应力耦合模型的准确性与适用性。     

沁水盆地煤样静水压力下渗透率实验及模型分析

为研究沁水盆地煤样渗透率演化规律,构建了煤样渗透率测定的瞬态压力脉冲法实验装置,使用N2和CO2在实验室开展了3种试验条件的渗透率测定,应用Connell模型对实验结果进行分析,并讨论了模型预测值和实验值之间差别的原因。结果表明,（1）在恒定孔压变围压条件下渗透率随有效应力增大而减小;在等有效应力条件下,渗透率随孔压增大而减少;在恒定围压变孔压条件下,随孔压增大,渗透率呈先减小后变大的趋势。（2）运用Connell模型预测的恒定围压变孔压条件的渗透率值大于实验值,原因可能是由于裂隙压缩性系数和吸附应变系数存在估计误差。通过开展实验室渗透率实验和模型分析,对指导实验室内二氧化碳封存和气体驱替实验及其模拟研究具有借鉴意义。     

Leaching and differential recrystallization of metamict zircon under experimental hydrothermal conditions

We report results of hydrothermal experiments on four alluvial zircons from Sri Lanka, which cover a wide range of radiation damage, at 450 °C and 1.3 kbar for 744 h with 2 M CaCl₂ solution as reactive fluid. After the hydrothermal treatment, the most metamict samples show micrometer-thick reaction rims, which surround apparently unreacted zircon, as revealed by cathodoluminescence (CL) and Nomarski differential interference contrast (NDIC) images. These rims have sharp, curved, and transgressive boundaries with unreacted zircon and are, in some cases, spread out along cracks. The thickness of reaction rims increases with increasing cumulated !-dosage of the starting materials. The reaction rims are strongly enriched in Ca (up to 7000 ppm) and a water species and depleted in radiogenic Pb, Zr, and Si, as revealed by electron microprobe analyses. A significant Th loss from the reaction rims was detected in the case of the most metamict sample, whereas U remained in the structure. FT-infrared spectrometry and X-ray diffraction measurements revealed that the bulk run products were recrystallized. Using micro-Raman spectrometry, we were able to demonstrate that differential recrystallization took place. The reaction rims are strongly recrystallized, whereas the unreacted grain interiors underwent only minor recrystallization. Recrystallization of the rims is accompanied by an enhancement of the integral CL intensity. It is suggested that recrystallization in the reaction rims was catalyzed by water infiltration and ion exchange and prevented significant congruent zircon dissolution under the given experimental conditions. Previous zircon studies have shown that (1) a transgressive morphology, (2) a reduced Th-U ratio, and (3) an enhanced CL emission are also characteristics of rims in zircons from high-grade metamorphic rocks. Based on these similarities between natural and experimentally produced rims, it is suggested that leaching-catalyzed recrystallization is an important alteration process in zircon under wet geological conditions and can account for many complex core-rim structures found in natural zircons. Furthermore, the strong enrichment of Ca in the reaction rims supports previous assumptions that high Ca concentrations in natural zircons are of secondary origin. It is suggested that lower U-Pb concordia intercept ages obtained from single-phase zircons with high Ca contents date a leaching event.     

Temporal relationship between granite cooling and hydrothermal uranium mineralization at Dalongshan in China: a combined radiometric and oxygen isotopic study

A combined study using multi-radiometric dating and oxygen isotopic geothermometry was carried out for Mesozoic quartz syenite, alkali-feldspar granite and associated hydrothermal uranium mineralization at Dalongshan in the Middle-Lower Yangtze valley of east-central China. Radiometric dating of the quartz syenite yields a whole-rock Rb–Sr isochron age of 135.6±4.3 Ma, a zircon U–Pb isochron age of 132.9±2.2 Ma, and K–Ar ages of 126±2, 118±3 and 94±4 Ma for hornblende, biotite and orthoclase, respectively. The alkali-feldspar granite yields a whole-rock Rb–Sr isochron age of 117.3±3.3 Ma, a zircon U–Pb isochron age of 114.7±2.1 Ma, and K–Ar ages of 112±2, 109±3 and 88±4 Ma for hornblende, biotite and orthoclase, respectively. Oxygen isotope thermometry for both granites gives temperatures of 685 to 720, 555 to 580, 435 to 460 and 320 to 330 °C, for hornblende, magnetite, biotite and orthoclase respectively, when paired with quartz. The systematic differences among the ages by the different techniques on the different minerals are used to reconstruct the cooling history of the granite. The results yield rapid cooling rates of 27.4 to 58.6 °C/Ma from 800 to 300 °C in the early stage, but slow cooling rates of 6.3 to 7.2 °C/Ma from 300 to 150 °C in the late stage. The regular sequence of oxygen isotope temperatures for the different quartz–mineral pairs demonstrates that diffusion is a dominant factor controlling the closure of both radiometric and O isotopic systems during granite cooling. Pitchblende U–Pb isochron dating yields an uranium mineralization age of 106.4±2.9 Ma, which is younger than the age of the granite emplacement and thus considerably postdates the time of magma crystallization, but is close to the closure time of the K–Ar system in the biotite. This points to a close relationship between granite cooling and ore-forming process. It appears that hydrothermal mineralization took place in the stage of slow cooling of the granite, whereas the rapid cooling of the granite was concurrent with the migration of hydrothermal fluids along fault structures. Therefore, the activity of the ore-forming hydrothermal system is temporally dictated by the cooling rates of the granite and may lag about 25 to 30 Ma behind the crystallization timing of associated granite.     

The Hydrothermal Na—Ca Alteration at the Marginal Part of the Florina Granite and the Associated Magnetite—Hematite Bands With Thorium and Uranium Mineralization

The Oxia mineralized granite is the product of differentiation in the external parts of the Florina magmatic mass. Acidic hydrothermal solutions either of magmatic or of meteoric origin reacted with the upper tectonically fractured parts of the Florina granite and became enriched in iron, thorium, uranium, zircon and rare-earth elements. The most abundant alteration minerals are sericite and quartz, while the minerals of the mineralization bands include magnetite, hematite, thorite, monazite and zircon. The outer parts of the Oxia granite made it easy the percolation of hydrothermal solutions from the deeper heater to the upper cooler parts of the granite which acted as a hot spot.     

煤矿采动条件下煤层气储层模型及应用

煤矿区煤层气开发受煤矿采动影响,为实现煤矿采动条件下煤层气井产能数值模拟,抽象概化了采动条件下煤层气开发的地质模型,构建了采动条件下煤层气储层的数学模型,并通过对CBM-SIM软件二次开发,实现了采动条件下煤层气储层的数值模拟。在建模和数值解算软件开发中,基于采动离层量变化曲线公式构建了采动条件下储层渗透率变化曲线公式,利用采动井水位变化规律构建了储层漏失水量变化公式,利用时间卡机制解决了煤储层渗透率及漏失水量的动态求解和循环迭代过程中作为系数和边界条件的调用赋值,实现了渗透率随采动影响的动态变化、储层水漏失降压和储层产气的耦合解算。应用开发的软件对淮南矿区某矿采动条件下煤层气抽采井生产数据进行历史拟合和产量模拟应用,预测煤层气产量曲线与实际生产曲线基本一致,判定系数达到0.92。     

单轴压缩下饱水花岗岩破裂过程声发射频谱特征试验研究

通过开展饱水花岗岩单轴压缩声发射试验,采用快速傅里叶变换提取声发射信号主频及次主频,进一步提出了综合表征声发射信号频谱特征的参数:主频比F ,即主频与次主频的比值,研究岩石破裂过程声发射频谱特征。研究结果表明:饱水花岗岩破裂过程声发射主频、次主频频带数量呈现“3→6→3”先增大后减小的规律,与主频分布特征相比,次主频低频区间占比减小,中、高频区间占比增大,次主频总体平均值比主频高约5 kHz。声发射主频比尸分布区间在弹性阶段达到最大:0 < F < 8 ,且F≠1,根据主频比与1 的关系,可将声发射信号划分为2 类:A 类信号(0 < F < 1 ),其主频小于次主频;B 类信号(1 < F < 8 ),其主频大于次主频。总体来看,饱水花岗岩破裂过程声发射A 类与B 类信号的数量关系约为3:2, A 类信号数量相对较多,特别是塑性阶段A 类信号占绝对优势,其数量为B 类信号的3.0?3.5倍。饱水花岗岩破裂过程声发射A 、B类信号事件率近似同步变化,二者均出现明显的平静期前兆特征,在实际应用时,可以对岩石破裂过程数量庞大且计算繁琐的声发射信号进行甄别,选择A、B 两类信号中的一种,进而实现对岩石破坏前兆信息的快速有效识别。