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)时也存在低估,可能是由雨滴谱资料观测误差导致拟合的系数偏小或雷达硬件造成的观测偏差等造成的。     

Introducing a transferable fishing day management regime for Pacific Island countries

On 1 December 2007, eight ‘Small Island Developing States’ in the Western and Central Pacific Ocean implemented a management regime restricting the total number of days fished by tuna purse seine vessels within their waters, commonly referred to as the Vessel Day Scheme (VDS). The VDS is seen as one component of management arrangements to reduce fishing mortality on bigeye and yellowfin tuna, constrain fishing effort, and increase the rate of return from access fees by Distant Water Fishing Nations.     

稳健估计方法提高数据处理精度

介绍了近代平差理论的稳健估计方法,编制稳健估计方法的程序,并通过实例验证,与最小二乘估计进行比较,表明稳健估计在水准网粗差探测和平差计算中优于最小二乘估计方法,并且能够定位粗差,从而进行消除或者减弱,得到较为干净的观测值。因此,稳健估计方法应用于测量平差具有一定的抵抗粗差的能力,从而可以提高数据处理的精度。     

2022年1月8日青海门源6.9级地震机器学习地震预警震级估计与现地阈值报警的回溯验证

2022年1月8日1时45分青海省海北州门源县发生6.9级地震, 周边地区普遍有感, 并导致多条高铁线路临时停运. 本文利用这次地震获取的大量烈度计加速度记录, 基于正在进行系统研发的机器学习地震预警方法模块, 对地震预警震级估计与现地阈值报警进行了回溯验证. 结果表明: 在地震发生后3.1 s, 震级估计为6.5级, 且震级估计误差不受信噪比和震中距变化的影响, 随着首台触发后时间的增加, 震级估计逐渐接近实际震级. 对于现地地震动速度峰值PGV(Peak Ground Velocity)预测, 各个台站在P波到达后3 s时, 预测PGV与观测PGV呈现1:1线性关系, 随着P波到达后时间窗的增加, 预测PGV逐步接近观测PGV, 且PGV预测误差不受信噪比和震中距变化的影响. 现地台站仪器烈度阈值设置为Ⅵ度时, 报警成功、误报、漏报的百分比分别为99.53%、0%、0.47%, 平均预警时间为19.62 s, 且地震烈度Ⅵ度区内没有发生误报和漏报; 现地台站仪器烈度阈值设置为Ⅶ度时, 报警成功、误报、漏报的百分比分别为99.77%、0%、0.23%, 平均预警时间为9.69 s, 且地震烈度Ⅶ度区内没有发生误报和漏报. 此次回溯验证结果表明: 机器学习方法在这次地震中可以得到鲁棒的震级估计和现地阈值报警结果, 并为该方法的在线测试以及中国地震预警系统升级提供可行性依据; 其次, 在这次地震事件中, 烈度计可为预警提供额外的作用, 这也为烈度计在未来地震预警的研究和应用中提供了更多的可能性.

     

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.     

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

双偏振雷达差分传播相移的五种滤波方法对比分析

为了对双线偏振雷达的差分传播相移ΦDP进行滤波降噪处理,提高双线偏振雷达识别粒子、估测降水的能力。通过对C波段双线偏振雷达差分传播相移分别采用滑动平均、中值滤波、FIR滤波、卡尔曼滤波和小波分析5种不同方法地介绍,利用中国气象科学研究院灾害天气国家重点实验室移动式C波段双线偏振雷达(POLC)2013年外场试验观测资料,对比分析了经5种方法处理后的差分传播相移ΦDP径向数据,及其分别在降水估测中应用的效果。结果表明:采用FIR滤波方法、卡尔曼滤波方法和小波分析方法对ΦDP径向数据进行处理的效果比传统的滑动平均和中值滤波方法好;其中小波分析处理的ΦDP径向数据波动最小,再经过变距离法拟合的KDP数据的负值能够得到最有效地抑制,最后应用到R(KDP)降水估测中的精度最高。     

临界值可变的抗差估计等价权函数

抗差估计等价权函数一般由正态分布统计量构造,其临界值(或称准则)一般由实际经验确定。首先分析了正态分布统计量和学生化残差统计量的区别,然后分别讨论了基于这两种统计量构造的等价权函数的区别。研究表明,利用学生化残差统计量构造的等价权函数以及顾及误差显著性水平确定的临界值,不仅考虑了观测误差的大小,而且还可以顾及了实际观测的图形强度和多余观测数,可以克服人为确定临界值可能带来的参数估计的有效性和抗差性方面的风险。     

京津冀科技合作网络的演变特征及影响因素

构建层次分明、结构合理的科技合作网络对于推动京津冀协同创新具有非常重要的作用。本文从基于知识创新的科研合作网络和基于技术创新的技术合作网络的角度,使用合作论文和合作专利数据,采用社会网络分析法对2013—2018年京津冀地区地级及以上城市之间科技合作网络的演变特征进行了分析,并基于半参数估计和面板计量模型实证检验了影响京津冀78个城市对之间科技合作的主要因素。主要结论为：① 京津冀科研合作网络快速成长,呈现北京与天津双核引领,以京津为主轴、京保石为次轴的空间结构。② 京津冀逐渐形成了以北京为主中心,天津、石家庄为次中心,廊坊、保定、沧州为三级枢纽,其他城市为节点的技术合作网络结构。③ 空间距离越近的城市之间越容易产生科技合作,高铁开通带来的时空距离压缩可以减弱科技合作的空间衰减系数;相对于科研合作,技术合作对于空间距离更为敏感;技术邻近性是促进城市间科技合作的主要驱动因子,尤其是对技术合作的促进效应更为明显;经济发展差距越小的城市之间相对容易产生技术合作。     

变窗宽核加权估计下的变形趋势拟合方法

变形监测数据处理的方法有很多,但这些方法对数据量及数据的采集方式有特定的要求,或者计算过程复杂.针对这些问题,提出了基于变窗宽核加权估计的变形趋势拟合方法,即先用较大窗宽的核加权估计去拟合变形的整体趋势,再用较小窗宽的核加权估计去拟合残余变形量——局部趋势.并针对这一新方法,提出了一种新的窗宽计算方法,即时序间隔标准差窗宽.以某大坝某一监测点32期的高程变形拟合为例,比较了不同的窗宽以及不同变窗宽组合的核加权拟合效果.结果表明,采用时序间隔标准差窗宽的核加权拟合比经验窗宽的拟合精度高;而基于变窗宽的核加权拟合比前两者精度更高.