Seismic signals in geothermal areas of active volcanism: a case study from ‘La Fossa’, Vulcano (Italy)

Since the last eruption (1888–1890) volcanism at Vulcano, Aeolian Archipelago, southern Tyrrhenian Sea, has taken the form of persistent fumarolic activity. The gas-vapour phases of the geothermal systems are mainly discharged within two restricted areas about 1 km apart from each other, in the northern part of the island. These areas are La Fossa crater, and the beach fumaroles of the Baia di Levante. Fluids released at the two main fumarolic fields display quite different chemical and temperature characteristics, implying different origins. The local seismicity essentially takes the form of discrete shocks of shallow origin (depth1 km) at La Fossa, usually with energy < 10¹³ ergs. They are thought to be related to the uprise of pressurized hot gases and vapours discharged at the crater fumaroles. The present investigation points to the existence of two principal categories of seismic events (called M-shocks and N-shocks). These are short events (normally < 10 s). M-type shocks are thought to be due to resonance vibrations within the interior of the volcano, probably driven by the excitation of shock-waves within cavities deeply affected by deposition and alteration of self-sealant hydrothermal minerals. N-type events display features that resemble those of volcano-tectonic earthquakes, but have no recognizable S-phases. Here they are tentatively attributed to microfracturing of rocks which have been extensively hydrothermally altered. Results of the present study permit a preliminary conceptual model of the local shallow seismic processes in the framework of geochemical modelling of fumarolic activity and geological inferences from geothermal drilling.     

Estimation of statistical parameters for censored lognormal hydraulic conductivity measurements

Field tests of hydraulic conductivity (e.g., injection test, pumping test, etc.) in low permeability formations are subject to censoring due to the detection limit of the instruments used. An iterative method of estimating the mean and variance of hydraulic conductivity data with a presumed log-normal distribution function is presented. This method accounts for the data that are actually below the lower detection limit (called truncated data) and thus gives distribution parameters that are more representative for the underlying distribution. The proposed method is then tested on two simulated normally distributed random datasets having different variances. The results show that the means and variances estimated by the proposed method are very accurate. Finally, the method is used to estimate the mean and variance of hydraulic conductivity data from single hole water injection tests in a fractured geological formation.     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

地下水非稳定流的LTFAM算法

渗流域内应用拉普拉斯变换(LT)建立相应的有限分析(FAM)方程,顾及渗流域内地下水流的初始条件和边界条件,可在LT空间构成一个封闭的以水头像函数为变量的线性方程组。将此方程组所得的解,通过Stehfest数值反演公式,可归化为时间域的解(水头)。由于时间t被隐含在数值方程内,从而克服了传统数值法按时段(△t)逐步迭代的缺陷,提高了计算效率,也为用嵌入法建立地下水流管理模型提供了一条捷径。     

江西省安福地区地热水化学特征研究

为研究江西安福地区水文地球化学特征及控制因素, 本文采集了15组样品, 采用水化学、同位素分析等方法进行研究。结果表明: 研究区地热水以Na-HCO3型水为主, 地下水以HCO3-Na·Ca型为主, 地表水由西南向东北从HCO3-Na·Ca向HCO3-Ca型水演化。水化学组分演化过程主要受岩石风化作用控制, 地层封闭性较差, 水中的Na+、HCO– 3、Sr2+来源于硅酸岩风化溶解。由稳定同位素特征可知, 研究区地热水补给来源为大气降水。研究区热储温度为49.8~101.4 ℃, 地热水循环深度为1 502.6~1 513.6 m。热水在沿断裂带上升过程中与浅层冷水发生混合, 其混合比例为76.1%~87.5%。研究成果为安福地区水循环演化提供依据, 有利于地热资源的合理开采与保护。     

地铁施工对济南白泉泉群流量的影响

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用FEFLOW软件建立地下水流数值模型。以规划地铁M1号线为研究对象，分析了济南东站、梁王站、梁王东站分别施工及3个站同时施工4种情景下，采用施工降水或施工降水+人工回灌两种施工方式对白泉泉群流量的影响。结果表明：单独采用施工降水的施工方式使得白泉泉群流量衰减，其中3个站同时施工对泉流量的影响最大，泉流量最大衰减达5.48%；各站分别施工时，济南东站对泉流量影响最大，泉流量较未施工时减少了0.043×10⁴ m³/d。采用施工降水+人工回灌的施工方式，能够有效缓解泉流量的衰减，各车站施工时的泉流量衰减由仅施工降水时的2.26%~5.48%降低至0.08%~1.21%。岩溶强径流带有利于地下水形成优势径流，促进白泉泉群补给，一定程度上缓解因地铁施工引起的泉流量衰减。     

雄安新区地热地质模型探究：来自地球物理的证据

雄安新区内地热资源丰富，区内有牛驼镇地热田、容城地热田和高阳地热田，地热资源开发利用较早，但是对其深部热源机制仍未形成统一观点。为了研究雄安新区内地热田深部热源机制，在新区及外围进行了深反射地震和长周期大地电磁探测，对取得的同剖面的深反射地震和大地电磁数据进行处理和综合解释，探明了研究区从地表至莫霍面范围内地质构造和电性结构。下地壳结构在深反射地震剖面与大地电磁剖面上有很好的对应关系。电阻率低值区对应着在深反射地震剖面上存在一系列反射同相轴，且同相轴可以延续到莫霍面，电阻率高值区对应着在深反射地震剖面上无明显连续反射同相轴，尤其是在莫霍面之上呈现地震反射近似"空白区"。结合区域地热资料构建了研究区深部地热地质模型，对新区内深部地热机制进行了解释。该模型为"二元"生热模型，其热源包含两个部分，深部地幔热源和地壳放射性元素衰变生热。放射性元素衰变生热占地表热流的接近30%，而幔源热流在地表热流中的占比可达约70%。在牛驼镇下方，莫霍面以上，由于地幔热物质上涌造成下地壳上隆，幔源岩浆底侵作用于下地壳形成了局部热异常，该热异常具有低速高导的地球物理特征，认为是牛驼镇地热田和容城地热田的深部热源；以区域断裂为热通道，大地热流由深部向上传导、扩散到牛驼镇凸起和容城凸起顶部，对碳酸盐岩储水层进行加热，形成地热储层；上覆新近系沉积地层是良好的热盖层。     

Global climatic change – the latest scientific understanding

IPCCs statement in its 1995 report (IPCC 1996) that a human influence was discernible in global climate has been widely quoted but often misunderstood. The character of the evidence underpinning this detection statement is explained so that its strengths and weaknesses can be better understood and the subtleties of its message better appreciated. To demonstrate the close linkage between the government-approved summary and the underlying chapters of the IPCC report the detailed evolution of the detection statement from first draft through to the form finally approved by the IPCC is described.     

漳州热田的对流热流和传导热流的研究

漳州地热系统属对流型地热系统.漳州热田是我国东南沿海地区目前所见热田中温度最高（121.5℃）的一个.地表热流值的研究表明,热田中心具有最大的实测热流值（359mw/m²）.本文根据热田内152个钻孔的测温资料和56块岩石样品的热导率数据,采用三种不同方法计算出漳州热田及其邻近地区的大地热流值,并讨论了热流值分布的特点.