Application of Linear Mean-Square Estimation in Ocean Engineering

The attempt to obtain long-term observed data around some sea areas we concern is usually very hard or even impossible in practical offshore and ocean engineering situations. In this paper, by means of linear mean-square estimation method, a new way to extend short-term data to long-term ones is developed. The long-term data about concerning sea areas can be constructed via a series of long-term data obtained from neighbor oceanographic stations, through relevance analysis of different data series. It is effective to cover the insufficiency of time series prediction method’s overdependence upon the length of data series, as well as the limitation of variable numbers adopted in multiple linear regression model. The storm surge data collected from three oceanographic stations located in Shandong Peninsula are taken as examples to analyze the number-selection effect of reference oceanographic stations (adjacent to the concerning sea area) and the correlation coefficients between sea sites which are selected for reference and for engineering projects construction respectively. By comparing the N-year return-period values which are calculated from observed raw data and processed data which are extended from finite data series by means of the linear mean-square estimation method, one can draw a conclusion that this method can give considerably good estimation in practical ocean engineering, in spite of different extreme value distributions about raw and processed data.     

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.     

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

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

内蒙古翁公敖包航磁异常的推断与解释

文章对内蒙古四子王旗—土默特右旗中北部近年发现的翁公敖包航磁异常进行推断与解释,在对找矿有望的航磁异常进行筛选后,对翁公敖包航磁异常进行了野外地质-地球物理查证,确定异常为铁矿所致;采用向上延拓、垂向导数、斜导数及小波变换等方法对航磁异常进行数据转换处理,以获取铁矿体的横向及纵向分布信息,取得了较好的效果。     

基于小波分析与Mann-Kendall法的岩溶大泉动态研究

研究地下水动态是认识地下水资源的有效手段。根据1956－2013年济南岩溶泉域大气降水及地下水水位动态监测资料,采用小波分析法、Mann-Kendall趋势检验、突变检验法研究了58个水文年泉水位对大气降水的响应,可以看出:（1）大气降水和泉水位呈现出多尺度的变化特征,长时间尺度上两者的变化周期基本相同,变化周期为16年和12年,说明大气降水对泉水位有直接影响;（2）在1956—2013年,济南泉域地下水水位具有0.65 m·（10a）-1的年际显著下降趋势,但降水具有12.65 mm·（10a）-1的不显著上升趋势,说明在人为因素影响下泉水动态的影响因素的权重发生了变化;（3）大气降水在1999年发生突变,1999年之后年降水为增加趋势;而地下水水位突变年份为1967年,1967年以后水位持续降低,2004年以后水位快速上升,泉水位未来趋势应与降水保持一致,呈上升趋势,说明大气降水并非泉水动态的唯一影响因素;（4）通过建立不同时段的多元回归模型,表明近58年来地下水水位的主要影响因素由大气降水到人工开采之间的转换,同时验证了小波分析和Mann-Kendall法研究地下水动态的适宜性和可靠性,也为济南市的保泉提供了参考依据。      

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.     

磁法三维建模在龙王沟磁铁矿勘探中的应用

以龙王沟磁铁矿高精度磁测数据为基础,采用小波多尺度分解方法,对化极后的磁法数据进行多尺度分解,通过功率谱法计算小波多尺度分解后的各阶细节异常的场源深度。结合地形、地质及钻孔资料,采用2.5D人机交互反演对磁测剖面上的磁铁矿体进行拟合,利用Golden Software公司开发的Voxler软件,在虚拟三维环境下实现二维小波多尺度分解数据的可视化,建立三维地质体模型,直观的预测磁铁矿体的空间范围和分布形态,实现地质、物探和地上、地下信息一体化集成展示,进行专业分析评价和辅助决策,可以有效的提高地质找矿效果。     

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

沿海沿江城市地下开挖产生的工程渣土含泥量大、含水率高且较松散,主要运往渣土场进行堆填处置。由于产量巨大而处置场地有限,许多渣土场在运营过程中存在堆填速度快、缺乏排水设施、超高超库容堆填等问题,容易引发堆填体失稳事故。目前对于非饱和工程渣土堆体在快速堆填过程中的失稳机制认识尚不清晰,尤其是对这一过程中的高饱和度工程渣土强度增长规律缺乏足够的认识。以深圳红坳渣土场填料——花岗岩风化料(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)时也存在低估,可能是由雨滴谱资料观测误差导致拟合的系数偏小或雷达硬件造成的观测偏差等造成的。