Calculating Auroral Oval Pattern by AE Index

The relationship between the auroral oval pattern, i.e., location, size, shape, and intensity, and the auroral electrojet activity index （AE index） is studied. It is found that the maximal auroral intensity is elliptically distributed, and the lengths of semimajor and semiminor axes are positively correlated to AE. The intensity along the normal of the auroral oval can be satisfyingly described by a Gaussian distribution, and the maximum and the full width at half maximum of the Gaussian distribution ave both positively correlated to AE. Based on these statistical results, a series of experimental formulas as a function of AE are developed to calculate the location, size, shape, and intensity of the auroral oval. These formulas are validated by the auroral images released by SWPC/NOAA.     

基于5分钟雷达资料降水估测模型的比较

For the Z-R relationship in radar-based rainfall estimation, the distribution of corresponding R values for a given Z value （or the corresponding Z value for a given R value） may be highly skewed. However, the traditional power-law model is physically deduced and fitted under the normal-distribution presumption of radar wave echoes associated with a rain rate value, and it may not be very appropriate. Considering this problem, the authors devised several generalized linear models with different forms and distribution presumptions to represent the Z-R relationship. Radar-reflectivity scans observed by a CINRAD/SC Doppler radar and 5-minute rainfall accumulation recorded by 10 ground gauges were used to fit these models. All data used in this study were collected during some large rainfalls of the period from 2005 to 2007. The radar and all gauges were installed in the catchment of the Yishu River, a branch of the Huaihe River in China. Three models based on normal distribution and a dBZ presumption of gamma distribution were fitted using maximum-likelihood techniques, which were resolved by genetic algorithms. Comparisons of estimated maximized likelihoods based on assumptions of gamma and normal distribution showed that all generalized linear models （GLMs） of presumed gamma distribution were better fitted than GLMs based on normal distribution. In a comparison of maximum-likelihood, the differences between these three models were small. Three error statistics were used to assess the agreement between radar estimated rainfall and gauge rainfall： relative bias （B）, root mean square error （RMSE）, and correlation coefficient （r）. The results showed that no one model was excellent in all criteria. On the whole, the GLM-based models gave smaller relative bias than the traditional power-law model. It is suggested that validations conducted in many previous works should have been made against a specific criterion but overlooked others.     

中国近海海平面变化趋势的对比分析

本文分析了近40年的中国近海验潮站资料海表面高度的线性变化趋势,并与卫星高度计资料进行了对比。通过对验潮站资料的分析发现,中国海域无论是近40年(1970~2013年)、还是近20年(1993~2013年)海平面均显著上升。各海区近20年的海平面上升有加速的趋势,且各时段上升速率大于全球平均海平面上升率。但是,受到海平面的年际和年代际变化的影响,近10年海平面上升趋势放缓。同时,本文也分析了不同季节海平面变化的趋势,北部海域秋季最大,冬季最小;南海海域春季最大,秋季最小。通过AVISO资料和验潮站资料的对比可以发现,AVISO资料在描述近20年海平面变化的线性趋势上与验潮站资料接近,较大的差异主要是由验潮站地表发生升降引起的。同时,通过对比也发现了用验潮站资料估算海域平均的海平面高度变化会有一定的误差,在黄海、渤海、东海海域验潮站估计的数值偏高,而在南海海域则偏低。     

赣南漂塘钨矿花岗岩成岩年龄与成矿年龄的精确测定

漂塘钨矿是赣南著名的钨锡多金属矿床,由漂塘本区、木梓园和大龙山三个矿段组成。多年以来,有关漂塘花岗岩成岩年龄,大多是利用邻近的西华山花岗岩的数据进行推断,而已有的年龄测试数据精度很难保证。对该矿床的成矿年龄则报道更少。本文利用单颗粒锆石U-Pb法对漂塘本区和木梓园的成矿花岗岩进行定年,获得了两条很好的U-Pb年龄谐和线,漂塘本区和木梓园成矿花岗岩的年龄分别为161.8±1Ma和 153.3±1. 9Ma。利用矿石中的白云母40Ar/39Ar同位素定年获得了漂塘本区成矿年龄为152±1.9Ma,利用矿石中的辉钼矿Re-Os同位素定年获得木梓园成矿年龄为151.1±8.5Ma。 以上年龄数据表明,矿床的成矿作用或与花岗岩侵位基本同时发生,或稍晚于花岗岩。赣南其他矿集区的一些代表性钨矿如淘锡坑、大吉山等的成矿与成岩也有这两种情况,可见这种成矿方式在赣南地区具有代表性。本文的年龄测定和近年来发表的一些年龄数据表明,赣南钨锡矿床属于华南燕山中期第一阶段陆壳重熔型花岗岩类大规模活动及其伴随的钨锡等金属大规模成矿作用的产物。     

基于MOS的广西北部湾沿海秋冬极大风速精细化预报

采用北部湾北部沿海6个自动站实况资料和欧洲中心细网格数值预报产品,根据模式输出统计法(MOS),对广西北部湾地区2012年—2016冬半年(9月到次年2月)日极大风速建立站点秋冬季不同预报时效的多元线性回归方程,并用2017冬半年的数据进行TS检验评分.结果表明,秋冬季日极大风速和代表引导气流的500hPa蒙古槽、代表...     

地表下沉系数选取的LWPSO-BP方法研究

针对BP神经网络自身收敛速度慢、容易陷入局部极小点的缺点,本文以线性下降惯性权重粒子群优化算法(LLWPSO)为前处理器,优化BP网络的权值和阈值,利用实测资料数据,建立LWPSO-BP的地表下沉系数预计模型,并同普通BP模型预计结果对比结果表明:LWPSO-BP神经网络不仅训练速度快,而且预测精度明显提高,该模型对地...     

CG-5重力仪格值系数修正与应用研究

在地表重力观测过程中,CG-5重力仪的格值系数会随着时间推移而发生变化.选用河北测区内2015—2017年相对重力观测数据,基于重力差法对两台编号为C859、C873的CG-5相对重力仪进行格值系数的误差计算及修正,研究其变化特征,并对修正前后不同时间尺度的重力场变化进行对比分析.结果表明:两台CG-5相对重力仪的格值系数在3年内存在不同程度的变化,均呈现整体上升的趋势,累计变化量分别为0.028％,0.034％;采用修正后的格值系数可得到点值精度更高的平差结果;不同时间尺度上的重力场变化相比修正前也出现明显差异,2017年河北大城ML3.6级地震发生在格值系数修正后0.5 a期重力变化零值线附近,这些结果表明了对CG-5重力仪进行格值系数修正的正确性和必要性.     

新的震级国家标准在陕西地震台网的应用

选取陕西地震台网及邻省共享台站记录到的2008年5月至2021年5月104次地震资料,应用新的震级国家标准GB 17740—2017《地震震级的规定》得到新国标震级。研究表明：①新的地方性震级M_L比原规范的M_L震级系统偏大0.3左右,这主要是二者的量规函数不同造成的;②宽频带面波震级M_{S （BB）}和矩震级M_W震级基本一致,无系统偏差,大多数地震震级偏差为-0.1—0.1;③宽频带面波震级M_{S （BB）}比面波震级M_{S （CENC）}整体偏小0.2,二者之间存在系统偏差;④绝大多数地震的面波震级M_{S （CENC）}比M_W偏大0.2—0.3,二者之间存在系统偏差。     

基于页岩线性岩石物理模型的TOC反演方法研究

总有机碳含量(TOC)是评价页岩有机质丰度和生烃能力的重要指标.常规TOC预测方法是通过对弹性参数体(例如密度反演体)进行线性拟合得到的,该方法具有地区经验性,并且缺乏有效的理论支撑.针对此问题,本文以岩石物理驱动的物性参数反演作为技术手段,首先对典型页岩的矿物组分以及孔隙流体特征进行分析,得出影响页岩弹性参数变化的四个主要因素分别是等效脆性矿物体积分数、黏土体积分数、TOC以及总孔隙度,进而可有效减少待反演参数个数,提高物性参数反演算法的可行性.在此基础上,将Grana提出的线性化岩石物理模型推广至页岩,建立线性化页岩岩石物理模型,并在贝叶斯理论框架下建立页岩物性参数反演,实现从弹性参数反演体直接提取TOC、总孔隙度以及等效脆性矿物体积分数等多个页岩物性参数,为后续的页岩气甜点(最佳的页岩气勘探与开发的区域或层位)评价提供可靠的数据支撑.模型数据和实际数据测试验证了该反演算法的有效性.     

三轴加载含水煤样波速与孔隙率关系实验研究

煤的孔隙率影响煤的吸附、渗透能力,且它也是衡量煤体内部结构及孔、裂隙发育状况的重要指标.而水力增透技术在瓦斯抽采中的作用效果常常受到水的影响,且瓦斯抽采量也会受到含水煤层的影响,为更客观地了解含水煤层状态下的孔隙率的变化情况,采用超声波透射法对赵固二矿二1煤层不同层理方向煤样在干燥及自然饱和吸水状态下进行三轴加载超声波测试试验,研究干燥及自然饱和水三轴加载状态下煤样超声波波速与孔隙率关系.研究结果表明:相同条件下,加载煤样饱和煤样波速大于干燥煤样,且加载煤样饱和状态下不同层理方向波速变化易区别.而加载煤样饱和煤样与干燥煤样孔隙率变化与波速变化相反,加载煤样饱和状态下不同层理方向孔隙率变化不易区别;通过分析加载条件下饱和煤样与干燥煤样波速与孔隙率关系,发现无论是干燥煤样还是饱和煤样在加载条件下都符合对数关系φ=-aIn(Vp)+b且a值与围压有关,b值与孔隙率变化量有关.通过研究含水煤样在不同加载条件下对煤体波速和孔隙率的影响以及波速和孔隙率之间的关系,为煤炭物性勘探技术提供有效试验数据与理论基础,也为煤炭物性预防安全开采提供理论基础.