2015—2019年辽宁省暴雨红色预警信号分布及其特征

利用2015—2019年辽宁省发布的暴雨红色预警信号和1605个自动站的分钟级降水资料,统计暴雨红色预警信号和短时大暴雨年际变化和时空分布,分析暴雨红色预警信号的高分布区、易发时段。结果表明：2015—2017年辽宁省暴雨红色预警信号发布站数逐年递增,最大值出现在2017年,发布站数为147个;2015—2018年预警信号准确率提升,提前时间略减少,最低值为2018年,提前时间为19 min;2019年比2018年暴雨红色预警信号发布站数减少59个,提前时间增加29 min;暴雨红色预警信号的空间分布为东南部地区多、中部地区少;暴雨红色预警信号多在夜间发布;在辽宁省发布的50%以上的暴雨红色预警信号中,降水量达到预警发布标准的时间滞后于最大雨强出现时间90 min,最大雨强出现时间为暴雨红色预警信号发布的重要指标。为了达到防灾减灾的服务效果,发布暴雨红色预警信号时,应充分考虑最大雨强出现时间、发布时机、短时大暴雨高发区及地形的影响。     

冰后期以来长江水下三角洲层序地层特征及沉积环境演化*

长江水下三角洲层序地层学研究有助于全面了解长江三角洲地层特征和沉积环境演化模式。通过对长江水下三角洲下切河谷区YD0901和YD0903孔岩心的详细沉积物粒度、特征元素比值(Cl/Ti和Zr/Rb)、沉积相对比分析,恢复了冰后期以来长江水下三角洲层序地层格架。研究区冰后期以来自下而上依次出现河流相、潮汐河流相、河口湾相、浅海相和三角洲相的沉积相序。末次冰期海平面下降,古长江形成下切河谷,古河间地发育硬黏土层,构成五级Ⅰ型层序界面。之后海平面回升,分别于15 cal ka BP和8.0 cal ka BP形成最大海退和最大海侵界面,水下三角洲区域最大海侵发生时间略滞后于平原区,约为7.5 cal ka BP。据此3个层序界面将冰后期地层划分为低位体系域、海侵体系域和高位体系域。钻孔岩心记录揭示了14.8 cal ka BP海侵到达研究区;14.8~13 cal ka BP期间,受MWP-1A冰融水事件影响海平面快速上升,海岸线向陆推进速率可达71.9,km/ka;海退期间各钻孔沉积速率较低,直至2 cal ka BP开始,沉积速率明显增加。     

基于有效邻域波场近似的起伏地表保幅高斯束偏移

随着我国陆上地震勘探向复杂地表探区的转移,高精度、适应性强的地震成像方法在地震资料的处理、解释及后续属性分析、储层预测中具有重要意义.本文基于有效邻域波场近似理论发展了一种成像精度更高且适用于复杂起伏地表条件的叠前保幅高斯束偏移方法.在传统水平地表高斯束偏移的基础上,本文根据中心射线附近有效邻域内高斯束表征的近似波场,导出了起伏地表条件下具有相对振幅保持的高斯束偏移公式,并给出了一种精度更高的旁轴射线传播角度计算方法.同现有的高斯束偏移方法相比,本文方法不仅考虑了起伏地表对高斯束走时的线性影响,而且首次引入了由地表高程差异和近地表速度变化引起的二次时差校正项和振幅校正项,使得成像结果更加准确可靠.两个典型模型算例验证了本文方法的正确性和有效性.     

Common‐shot Fresnel beam migration based on wave‐field approximation in effective vicinity under complex topographic conditions

Gaussian beam migration is a versatile imaging method for geologically complex land areas, which overcomes the limitation of Kirchhoff migration in imaging multiple arrivals and has no steep‐dip limits of one‐way wave‐equation migration. However, its imaging accuracy depends on the geometry of Gaussian beam that is determined by the initial parameter of dynamic ray tracing. As a result, its applications in exploration areas with strong variations in topography and near‐surface velocity are limited. Combined with the concept of Fresnel zone and the theory of wave‐field approximation in effective vicinity, we present a more robust common‐shot Fresnel beam imaging method for complex topographic land areas in this paper. Compared with the conventional Gaussian beam migration for irregular topography, our method improves the beam geometry by limiting its effective half‐width with Fresnel zone radius. Moreover, through a quadratic travel‐time correction and an amplitude correction that is based on the wave‐field approximation in effective vicinity, it gives an accurate method for plane‐wave decomposition at complex topography, which produces good imaging results in both shallow and deep zones. Trials of two typical models and its application in field data demonstrated the validity and robustness of our method.     

Effect of pore geometry on Gassmann fluid substitution

Although there is no assumption of pore geometry in derivation of Gassmann's equation, the pore geometry is in close relation with hygroscopic water content and pore fluid communication between the micropores and the macropores. The hygroscopic water content in common reservoir rocks is small, and its effect on elastic properties is ignored in the Gassmann theory. However, the volume of hygroscopic water can be significant in shaly rocks or rocks made of fine particles; therefore, its effect on the elastic properties may be important. If the pore fluids in microspores cannot reach pressure equilibrium with the macropore system, assumption of the Gassmann theory is violated. Therefore, due to pore structure complexity, there may be a significant part of the pore fluids that do not satisfy the assumption of the Gassmann theory. We recommend that this part of pore fluids be accounted for within the solid rock frame and effective porosity be used in Gassmann's equation for fluid substitution. Integrated study of ultrasonic laboratory measurement data, petrographic data, mercury injection capillary pressure data, and nuclear magnetic resonance T₂ data confirms rationality of using effective porosity for Gassmann fluid substitution. The effective porosity for Gassmann's equation should be frequency dependent. Knowing the pore geometry, if an empirical correlation between frequency and the threshold pore‐throat radius or nuclear magnetic resonance T₂ could be set up, Gassmann's equation can be applicable to data measured at different frequencies. Without information of the pore geometry, the irreducible water saturation can be used to estimate the effective porosity.     

末次冰盛期时吉兰泰盐湖的湖泊状态与古气候特征

通过现代季风边缘区的吉兰泰盐湖钻孔JLT11-A孔沉积岩芯的矿物分析,结合地层盘星藻的含量探讨末次冰盛期(LGM)时湖泊的状态和古气候特征。结果显示:在LGM时吉兰泰湖泊沉积矿物主要是石英、长石为主的碎屑岩沉积,含量在85%左右,显示出陆源碎屑矿物的高输入状态,可能指示区域寒冷干旱的环境;其次是以方解石为主的碳酸盐的沉积,含量约为10%;氯化物为主的石盐类矿物一般不足5%,但持续存在,指示湖泊仍然有较高的盐度,因此地层中的淡水藻类盘星藻可能是由河流输入。由于陆源碎屑矿物输入强烈,矿物组合可能难以直接指示吉兰泰盐湖湖水状态。区域的干冷的气候与大多数的古气候记录一致,而与新疆西部的冷湿的环境不同。对比邻近区域的古气候研究结果发现,本区域在LGM时段夏季降水相对于冰消期偏多,而相对于MIS3阶段晚期偏少,整体夏季风减弱。吉兰泰盐湖末次冰盛期到末次冰消期以来矿物组合的变化表明,湖泊环境可能受到夏季太阳辐射、全球与区域温度变化以及夏季风强弱变化的影响。     

最大准则优化技术在贴体网格中的应用

贴体网格在地质数值模拟中具有广阔的应用前景,为解决贴体网格生成时边界离散问题,提出了最大长度准则和最大面积准则,把曲线逼近和曲面网格优化问题转化为数学优化问题,为求解该问题,提出了改进的单粒子优化算法.试验表明,最大长度准则和最大面积准则的优化效果好于常规方法;以改进的单粒子优化算法求解该问题时,计算效率是智能单粒子优化算法的30倍左右（节点量为200）,从而实现最大长度准则和最大面积准则在贴体网格生成中的应用.针对最大面积准则优化曲面网格不能控制网格步长的情况,提出了限定步长的网格优化算法,使网格步长合理化,并通过实例验证了该算法的有效性.研究成果提供了生成贴体网格时边界优化准则和求解方法,对今后复杂边界的贴体网格生成具有重要意义.     

基于相对质心震中的地震破裂方向性测定方法研究:以2008年云南盈江MS6.0地震为例

地震破裂方向性参数包含了断层几何形态、破裂长度等信息,对于地震灾害评估、孕震机理研究具有重要意义.但是基于点源近似的震源矩中心张量机制解(CMT)或断层面解(fault plane solution)只能给出两个节面,无法确定破裂断层面.本文提出一种基于质心与起始震中的差异确定发震断层的方法,该方法利用P波到时来测定主震与参考地震间相对起始震中,同时利用波形反演震源机制过程中的主震时移与参考地震时移之差测定相对质心震中,在假定单侧破裂的情形下,根据时移之差随方位角的变化推断破裂断层面.本文使用该方法研究了2008年云南盈江M_S6.0走滑型地震,发现其破裂方向与其他方法结果一致.由于该方法基于相对起始震中和质心震中,可有效削弱速度结构模型不准确以及地震绝对位置误差带来的影响,应可适用于其他类似类型6级左右地震的破裂方向性研究,但仍需进一步工作对其进行验证.     

接地导线源电磁场全域有效趋肤深度

为了更好地推进广域电磁法的发展和应用,本文以接地长导线源为例,研究了可控源电磁场全场域的有效趋肤深度.利用频域电偶极源在均匀半空间产生电磁场的闭合表达式,计算了不同电磁场分量定义的有效趋肤深度,并讨论了在不同频率、不同偏移距、不同电导率情况下,有效趋肤深度的变化特性.根据有效趋肤深度随偏移距的变化特征以及与平面波趋肤深度之间的关系,利用多项式拟合的办法在五个不同频率范围内给出了适用于全场域的有效趋肤深度快速估算公式.研究结果表明：不同分量定义的有效趋肤深度是不同的,但是它们随偏移距、频率、电导率等参数的变化趋势是类似的,而且在一定范围内都趋近于平面波趋肤深度.理论模型的研究表明,有效趋肤深度可以作为测量参数选取和数据解释工作的参考依据.     

中亚地震目录震级转换及其完整性分析

为编制能够应用于地震危险性分析的中亚地区统一震级标度为矩震级的地震目录,从国际地震中心（ISC）下载得到该地区1907-2012年的地震数据,该数据包含众多机构不同震级标度的地震记录.以MW,GCMT为参考机构震级标度,并用最小二乘法拟合了其他机构震级标度与MW,GCMT之间的转换关系.在挑选某次地震事件的唯一震级记录时,以对应转换关系的相关系数大且剩余标准差小为准则进行筛选.对于少量没有震级转换关系或者拟合优度过差的地震记录,使用间接转换关系或者全球转换关系予以补充转换.震级转换后,用时空窗法删除前余震,并考虑构造环境和地震活动水平的空间差异性将研究区划分成5个子区域.采用地震记录时间累积曲线法、最大曲率法（MAXC）和拟合优度检验法（GFT）综合分析各个分区的最小完整震级（MC）,并在此基础上用极大似然法拟合相应的地震活动性参数.结果表明,每种完整性分析方法各具一定的优缺点,但采用综合分析的方法能够得出最佳的MC.地震记录时间累积曲线法能分析出高质量地震目录的起始时间,以作为后两种方法的基础,但容易受到地震活动水平随时间波动的影响.由于研究区目录质量较差,最大曲率法误差过大以至于只能作为其他方法的补充,拟合优度检验法的GFT参数也普遍只能达到80％左右.GFT最大值点可能并不对应MC,但是MC通常都在GFT极大值点取得.5个分区1964-2012年的MC普遍在Mw4.8左右,b值在1.136-1.514之间波动.