Investigation of gas content of organic-rich shale:A case study from Lower Permian shale in southern North China Basin,central China

Measuring gas content is an essential step in estimating the commerciality of gas reserves. In this study,eight shale core samples from the Mouye-1 well were measured using a homemade patented gas desorption apparatus to determine their gas contents. Due to the air contamination that is introduced into the desorption canister, a mathematical method was devised to correct the gas quantity and quality.Compared to the chemical compositions of desorbed gas, the chemical compositions of residual gas are somewhat different. In residual gas, carbon dioxide and nitrogen record a slight increase, and propane is first observed. This phenomenon may be related to the exposure time during the transportation of shale samples from the drilling site to the laboratory, as well as the differences in the mass, size and adsorptivity of different gas molecules. In addition to a series of conventional methods, including the USBM direct method and the Amoco Curve Fit(ACF) method, which were used here for lost gas content estimation, a Modified Curve Fit(MCF) method, based on the 'bidisperse' diffusion model, was established to estimate lost gas content. By fitting the ACF and MCF models to gas desorption data, we determined that the MCF method could reasonably describe the gas desorption data over the entire time period, whereas the ACF method failed. The failure of the ACF method to describe the gas desorption process may be related to its restrictive assumption of a single pore size within shale samples. In comparison to the indirect method, this study demonstrates that none of the three methods studied in this investigation(USBM, ACF and MCF) could individually estimate the lost gas contents of all shale samples and that the proportion of free gas relative to total gas has a significant effect on the estimation accuracy of the selected method. When the ratio of free gas to total gas is lower than 45%, the USBM method is the best for estimating the lost gas content, whereas when the ratio ranges from 45% to 75% or is more than 75%, the ACF and MCF methods, are the best options respectively.     

黔东松桃地区燕山期构造特征研究——以高地超大型锰矿床为例

黔东松桃地区构造格局主要是受到燕山运动影响形成。燕山早期,区内处于SEE—NWW向区域挤压收缩,形成了NNE—SSW向褶皱和断裂。下寒武统杷榔组—变马冲组地层和南华系下统大塘坡组地层作为两套滑脱层,在构造变形过程中起到滑脱作用,滑脱层内褶皱变形强烈,以大量紧闭、顶厚不协调褶皱为主要构造样式,而其上、下岩层则以断块或宽缓褶皱变形为主,以软弱层为界面,其上、下岩层发生逆冲滑脱作用,形成了区内隔槽式褶皱的构造样式。燕山晚期,区内处于NNE—SSW挤压为主,形成了区内NWW向小型褶皱和断裂,同时燕山晚期构造横跨叠加在早期变形之上,导致早期NNE向褶皱枢纽呈波状起伏,轴迹弯曲,部分地区形成了构造穹窿和构造盆地,部分早期断层呈左行走滑性质。在燕山早、晚两期构造作用下,高地矿区含锰岩系作为软弱滑脱层,向猴子坳向斜的核部进行滑脱聚集,使其在向斜的核部富集,翼部相对减薄。同时,由于犁式正断层的存在,也影响了锰矿层的保存与破坏。     

融雪洪水跳跃变异对其参数估计不确定性影响分析

利用Pettitt非参数检验法和Mann-Kendall非参数趋势检验法,分析年最大洪峰流量序列的非一致性,确定序列的变异形式,采用“分解-合成”理论对其进行一致性修正,得到过去、现状两种条件下年最大洪峰流量序列,根据贝叶斯理论对序列一致性修正前后参数不确定性进行估计,并对其预报区间优良性进行评价。研究结果表明：年最大洪峰流量序列变异点发生在1993年,序列整体上升趋势不显著,在1957-1993年子序列呈显著下降趋势,而1994-2006年子序列变化趋势不显著,跳跃变异为序列主要变异形式;给出了实测、还原及还现序列参数后验分布估计值及95%置信区间,将其结合优化适线法进行P-Ⅲ型频率分析,得到修正前后设计频率年最大洪峰流量预报区间估计值;还原、还现序列与实测序列相比,预报区间覆盖率均提高24%,平均带宽分别减少39.59%、23.17%,平均偏移幅度分别减少28.45%、11.39%。通过对非一致性年最大洪峰流量序列还原/还现计算,可减小参数估计不确定性对其计算产生的影响,从而提高预报区间的可靠性。     

Satellite estimates and subpixel variability of rainfall in a semi-arid grassland

卫星估雨精度的不确定性受到当地降雨类型和像元内降雨非均匀性影响,而结合这两个关键因素开展半干旱草原卫星估雨的研究有限.2009年夏,我们在中国锡林郭勒半干旱草原用多部微雨雷达和雨量计构建了9 km卫星像元降雨观测网,观测了像元内降雨非均匀性(空间变异系数CV),并评估了卫星估雨精度.结果表明:(1)CV值受像元内平均降雨量,降雨类型,降雨云面积及移向等影响,如高Cv值的降雨过程大多为平均降雨量小,对流性降雨过程,降雨云边缘像元CV值较高;(2)TRMM 3B42V7卫星估雨产品适用性较好,CMORPH和PERSIANN次之,但TRMM 3B42V7易在半干旱草原湖泊处高估降雨.     

A CRPS-Based Spatial Technique for the Verification of Ensemble Precipitation Forecasts

Traditional precipitation skill scores are affected by the well-known“double penalty”problem caused by the slight spatial or temporal mismatches between forecasts and observations. The fuzzy (neighborhood) method has been proposed for deterministic simulations and shown some ability to solve this problem. The increasing resolution of ensemble forecasts of precipitation means that they now have similar problems as deterministic forecasts. We developed an ensemble precipitation verification skill score, i.e., the Spatial Continuous Ranked Probability Score (SCRPS), and used it to extend spatial verification from deterministic into ensemble forecasts. The SCRPS is a spatial technique based on the Continuous Ranked Probability Score (CRPS) and the fuzzy method. A fast binomial random variation generator was used to obtain random indexes based on the climatological mean observed frequency, which were then used in the reference score to calculate the skill score of the SCRPS. The verification results obtained using daily forecast products from the ECMWF ensemble forecasts and quantitative precipitation estimation products from the OPERA datasets during June-August 2018 shows that the spatial score is not affected by the number of ensemble forecast members and that a consistent assessment can be obtained. The score can reflect the performance of ensemble forecasts in modeling precipitation and thus can be widely used.     

粤港澳大湾区地下水资源特征及开发潜力

粤港澳大湾区是中国开放程度最高、经济活力最强的区域之一,在国家发展大局中具有重要的战略地位,大湾区未来的发展离不开水资源的支撑和良好的水生态环境。近年来,随着大湾区经济的快速发展,人口激增,需水量上升,水资源环境问题也日益突出,水安全保障程度不足;地下水是水资源的重要组成部分,具有水量稳定、水质较好的特点,可作为重要的应急备用水源。本文从地下水资源着手,系统梳理了大湾区水资源环境条件、地下水资源状况、特征和开发利用潜力,并提出了应急后备水源地建议,得到以下认识:(1)地下水可划分为松散岩类孔隙水、碳酸盐岩岩溶水、基岩裂隙水三大类,其中松散岩类孔隙水和基岩裂隙水分布最广;(2)湾区内地下水水化学类型较为复杂,丘陵山区以HCO_3-Na型、HCO_3-Ca型及HCO_3-Na+Ca型为主,冲积平原及山间盆地以HCO_3+Cl-Na型及HCO_3+Cl-Na+Ca型为主,三角洲地区以Cl-Na型微咸-咸水为主;(3)西江、北江及东江干流构成湾区内地下水排泄的总渠道,各支流为地下水的局部排泄基准面,地下水动态变化具季节性特征;(4)地下水整体水质较好,Ⅰ-Ⅲ类水占比高达66.25%,从丘陵山区到三角洲平原,水质呈变差趋势,尤其在广州、江门、中山、东莞等城市周边水质较差,超标因子主要为氨氮、氯化物、氟化物、硫酸盐等,三角洲地区发育大量"铁质水"和"氨氮水",水质性缺水问题突出;(5)地下水开发利用程度很低,东莞及中山等城市基本未开发利用地下水,在各类地下水中,碳酸盐岩岩溶裂隙水具有规模开采的开发利用潜力;(6)综合分析相关资料,提出将广花盆地等10处富水块段作为应急水源地备选,经初步计算每年可为大湾区提供约4.18亿m~3的应急水源保障。为应对突发性水质污染及极端干旱气候等大规模供水危机,保障粤港澳大湾区用水安全,促进大湾区高质量发展,建议加强大湾区的基础水文地质调查工作,掌握地下水的水位、水质、水量的动态变化特征,精准计算可用于应急备用开采的地下水储存量。     

花岗岩风化料弃土快速堆填过程中不排水抗剪强度评估

沿海沿江城市地下开挖产生的工程渣土含泥量大、含水率高且较松散,主要运往渣土场进行堆填处置。由于产量巨大而处置场地有限,许多渣土场在运营过程中存在堆填速度快、缺乏排水设施、超高超库容堆填等问题,容易引发堆填体失稳事故。目前对于非饱和工程渣土堆体在快速堆填过程中的失稳机制认识尚不清晰,尤其是对这一过程中的高饱和度工程渣土强度增长规律缺乏足够的认识。以深圳红坳渣土场填料——花岗岩风化料(CDG)填土为研究对象,对不同初始饱和度土样进行三轴不排水不排气等向压缩和剪切试验,结果表明:非饱和CDG填土不排水抗剪强度随围压增大呈非线性增长,增长速率与试样初始饱和度密切相关;当土样压缩后的饱和度超过0.7,不排水强度随围压的增长速率显著降低。基于有效固结应力法的原理,结合Hilf孔压公式和修正剑桥模型,提出了一种工程渣土不排水抗剪强度估算方法,并通过与试验结果对比验证了该方法对初始饱和度高于0.6的CDG填土的适用性。利用该方法确定的不排水强度cu与正应力σn的关系可应用于高饱和度工程渣土快速堆填中的稳定性分析。     

CINRAD-SA双偏振雷达资料在降水估测中的应用初探

对基于水平反射率ZH和差分传播相移率K_(DP)的降水估测综合法R(C)进行了改进,并对广州S波段双偏振雷达2016年2次飑线和2次台风降水过程的Φ_(DP)使用小波分析进行滤波处理,在此基础上使用变距最小二乘法拟合得到K_(DP)的值。分别使用R(C)和R(Z_H)法对2次飑线和2次台风降水过程进行降水估算,将估算结果和雨量计小时雨量进行了对比,并将两种方法的评估结果进行了对比。结果表明:(1)对于飑线类型降水,R(C)法对5 mm·h~(-1)以上的降水估测精度要好于R(Z_H)法,且降水率越大,R(C)法优势越明显,当降水率≥20 mm·h~(-1)时,两次过程R(C)法比R(Z_H)法的平均相对误差(RE)降低了17. 2%,平均绝对误差(AE)减少了1.89 mm,平均均方根误差(RMSE)减少了1.66 mm;(2)对于台风类型降水,R(C)法对5 mm·h~(-1)以上的降水估测精度也好于R(Z_H)法,当降水率≥20 mm·h~(-1)时,两次过程R(C)法比R(Z_H)法的平均RE降低了33. 19%,平均AE减少了3. 95 mm,平均RMSE减少了4.05 mm;(3)对于飑线和台风两种类型降水R(C)法都明显改善了降水率较大时的R(Z_H)法低估问题,但R(C)法在降水率10 mm·h~(-1)时也存在低估,可能是由雨滴谱资料观测误差导致拟合的系数偏小或雷达硬件造成的观测偏差等造成的。     

Assessing uncertainties in a conceptual water balance model using Bayesian methodology / Estimation bayésienne des incertitudes au sein d’une modélisation conceptuelle de bilan hydrologique

Abstract

Abstract The aim of this study was to estimate the uncertainties in the streamflow simulated by a rainfall–runoff model. Two sources of uncertainties in hydrological modelling were considered: the uncertainties in model parameters and those in model structure. The uncertainties were calculated by Bayesian statistics, and the Metropolis-Hastings algorithm was used to simulate the posterior parameter distribution. The parameter uncertainty calculated by the Metropolis-Hastings algorithm was compared to maximum likelihood estimates which assume that both the parameters and model residuals are normally distributed. The study was performed using the model WASMOD on 25 basins in central Sweden. Confidence intervals in the simulated discharge due to the parameter uncertainty and the total uncertainty were calculated. The results indicate that (a) the Metropolis-Hastings algorithm and the maximum likelihood method give almost identical estimates concerning the parameter uncertainty, and (b) the uncertainties in the simulated streamflow due to the parameter uncertainty are less important than uncertainties originating from other sources for this simple model with fewer parameters.     

On-line estimation of nonlinear physical systems

Recursive algorithms for estimating states of nonlinear physical systems are presented. Orthogonality properties are rediscovered and the associated polynomials are used to linearize state and observation models of the underlying random processes. This requires some key hypotheses regarding the structure of these processes, which may then take account of a wide range of applications. The latter include streamflow forecasting, flood estimation, environmental protection, earthquake engineering, and mine planning. The proposed estimation algorithm may be compared favorably to Taylor series-type filters, nonlinear filters which approximate the probability density by Edgeworth or Gram-Charlier series, as well as to conventional statistical linearization-type estimators. Moreover, the method has several advantages over nonrecursive estimators like disjunctive kriging. To link theory with practice, some numerical results for a simulated system are presented, in which responses from the proposed and extended Kalman algorithms are compared.