联合多点地质统计学与序贯高斯模拟的随机反演方法

岩相和储层物性参数是油藏表征的重要参数,地震反演是储层表征和油气藏勘探开发的重要手段.随机地震反演通常基于地质统计学理论,能够对不同类型的信息源进行综合,建立具有较高分辨率的储层模型,因而得到广泛关注.其中,概率扰动方法是一种高效的迭代随机反演策略,它能综合考虑多种约束信息,且只需要较少的迭代次数即可获得反演结果.在概率扰动的优化反演策略中,本文有效的联合多点地质统计学与序贯高斯模拟,并结合统计岩石物理理论实现随机反演.首先,通过多点地质统计学随机模拟,获得一系列等可能的岩相模型,扰动更新初始岩相模型后利用相控序贯高斯模拟建立多个储层物性参数模型;然后通过统计岩石物理理论,计算相应的弹性参数;最后,正演得到合成地震记录并与实际地震数据对比,通过概率扰动方法进行迭代,直到获得满足给定误差要求的反演结果.利用多点地质统计学,能够更好地表征储层空间特征.相控序贯高斯模拟的应用,能够有效反映不同岩相中储层物性参数的分布.提出的方法可在较少的迭代次数内同时获得具有较高分辨率的岩相和物性参数反演结果,模型测试和实际数据应用验证了方法的可行性和有效性.     

2013年南斯科舍海岭MW7.8地震的多点震源机制反演

2013年11月17日,在南极南奥克尼群岛北、南极板块与斯科舍板块之间发生了一次M_W7.8级地震（2013年南斯科舍海岭M_W7.8地震）,我们利用全球分布的长周期和宽频带地震记录反演确定了这次地震随时间和空间变化的震源机制,验证了提出的一种多点震源机制反演的新方法.首先利用长周期记录的W震相反演了这次地震的矩心矩张量解并利用体波提取了视震源时间函数,同时利用台阵反投影技术从宽频带记录中获得了这次地震的高频源的时空分布,然后基于矩心矩张量解、视震源时间函数以及高频源的时空分布,实现了采用新方法对2013年南斯科舍海岭M_W7.8地震的多点震源机制反演.矩心矩张量解表明,地震矩心在44.50°W/60.18°S,矩心深度19 km,半持续时间49 s,释放标量地震矩4.71×10²⁰ N·m,发震断层走向104°,倾角54°,滑动角8°.视震源时间函数清楚地揭示了地震矩随时间变化的方位依赖性,总体上可以将时间过程分为前60 s和后50 s两个阶段,但前60 s可细分为两次子事件.根据台阵反投影结果,这次地震为沿海沟从西到东的单侧破裂,破裂长度达311 km,可以分为5次子事件,能量释放的峰值点依次为13 s、30 s、51 s、64 s和84 s,平均破裂速度分别为0.6 km·s^-1、2.6 km·s^-1、2.3 km·s^-1、2.8 km·s^-1和3 km·s^-1.多点震源机制反演显示,5次子事件的矩震级分别为M_W7.57,M_W7.48,M_W6.80,M_W7.53和M_W7.08,半持续时间依次为21 s,17 s,6 s,16 s和8 s,走向分别为95°,105°,81°,98°和98°,倾角依次为57°,49°,86°,46°和64°,滑动角-9°,1°,-17°,13°和-4°.这些在震源机制、能量释放以及持续时间方面的变化都是当地构造和应力环境复杂性的反映.     

频率域海洋可控源电磁垂直各向异性三维反演

地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性,若仅对海洋可控源电磁(MCSEM)数据进行常规各向同性反演,有可能无法获得准确的反演解释结果,从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性(VTI)条件下频率域海洋可控源电磁数据三维反演算法.其中,三维正演采用基于二次场控制方程的交错网格有限体积法,并利用直接矩阵分解技术来求解离散所得的大型线性方程组,有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后,对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演,结果表明:各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释,而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布,展现出对海底电性各向异性结构更为优良的反演解释能力.     

北京市秋季城区和郊区大气边界层参数观测分析

应用2001年9月北京城区和郊区同步大气边界层观测资料,对大气边界层热力和动力参数进行了计算分析.结果表明:北京市秋季,逆温出现的时间城区滞后于郊区;逆温层高度城区大于郊区,200m通常为郊区的逆温层顶和城区的逆温层底部;逆温层强度城区弱于郊区;城市热岛强度为3℃;城区感热通量的输送大于郊区;城区下垫面粗糙度远大于郊区;郊区总辐射和紫外辐射的强度明显大于城区.     

南疆沙漠腹地大气边界层气象要素廓线分析

利用塔中80m观测塔梯度系统采集的2006年8月、10月和2007年1月、4月的风、温度、湿度资料,结合气象站的同步气象资料,对南疆沙漠腹地近地层四季的晴天平均风速、温度、湿度廓线分布特征进行分析。结果表明,晴天平均风速白天随高度升高增加缓慢,夜间较快,低层风速白天比夜间大,高层则白天比夜间小,春夏季风速较大;四季平均温度廓线表现为夜间辐射型、早上过渡型、白天日射型和傍晚过渡型等四种类型,早、晚过渡时间四季各有不同,日最低、最高温度出现时间四季则相差不大;冬季夜间比湿随高度升高而增大,整个80m近地层表现为逆湿状态,其他季节逆湿一般出现在0．5—1m、1～2m、32—47m、63—80m等4个层次上,各逆湿层出现的时间各季节有所差异。     

湖北省冻雨天气特征及其判定方法研究

利用常规地面气象站多时次气象要素和天气现象观测资料,以及武汉、宜昌和恩施的探空资料,对1998—2014年湖北地区冻雨天气的主要环流背景进行了分析,归纳了冻雨发生发展过程的气象要素和大气层结特征,给出了基于正负能量面积的冻雨发生定量化判断方法。结果表明:乌山阻高型和贝湖阻高型是导致湖北省出现冻雨天气的两种主要天气形势;2005年和2月是冻雨主要发生的年份和月份;共存在3种典型的气温层结演变特征,其中"暖雨-冻雨-固态"的气温层结演变最易导致持续时间长、影响范围大的冻雨天气;基于正负能量面积的冻雨判定方法为:3层层结时,正能量面积(A_(SP))小于80℃·hPa,负能量面积(A_(SN))小于400℃·hPa,且正负能量关系为(5.71A_(SP)-257.14)≤A_(SN)≤(6.25A_(SP)+200);2层层结时,ASP为350~650℃·hPa,A_(SN)为200~400℃·hPa,且正负能量关系为(A_(SP)-350)≤A_(SN)≤(ASP-100)。     

上海一次龙卷风过程分析

在分析环流背景的基础上结合风廓线仪探测资料对“990 90 6”龙卷风作了较为详细的物理剖析。分析发现 ,途经上海市的 990 9号热带风暴的残留涡旋由于西风槽的切入 ,动量迅速增大 ,形成中尺度强对流辐合云团 ,其中生成了 4个龙卷风暴。龙卷风暴在成熟阶段其环境气团的低空暖气柱逆温结构十分明显 ,垂直测风突然逆转了 36 0°。沙氏指数 SI小于 - 6℃ ,气层极不稳定。     

绥中一次暴雪伴雷电天气过程的成因分析

为了探讨绥中一次暴雪伴雷电天气过程的成因,利用常规观测资料、NCEP每6h间隔的1°×1°的再分析资料和营口多普勒雷达的资料,分析此过程的天气形势特点、高低空急流的作用、雷达回波的特征及反映动力、热力和水汽条件的相关物理量场的特征。结果发现:雷电发生在对流层中层的西南风急流和底层偏东风均处在最强的时刻,当对流云团发展到-20℃温度层时,温差起电产生雷电;雷电发生在低层850hPa附近存在的逆温层消失之后,同时配合低层水汽的辐合,产生了暴雪天气;雷电和强降雪发生在大气底层南风和北风转换的过程中,强降雪的时间与冷空气扩散加强的时间比较一致,当冷空气扩散到整个大气底层时强降雪结束;引起雷电和强降雪的对流不稳定层结主要处在对流层中层,并为上升运动的发生提供了动力和热力条件,促使雷电发生和强降雪的维持。     

青岛一次沙尘污染事件的气象条件特征

2008年5月28日至6月1日青岛市区出现了持续5 d的大气污染,其中29-30日为重度污染,其余3日为轻度污染。本文基于青岛常规地面和高空观测资料以及环境监测站SO2、NO2、PM10监测资料,利用资料分析和中尺度数值模拟的方法,分析造成青岛此次持续多日的大气污染的污染源、大气环流和气象要素特征。分析结果表明：此次污染过程主要是外来沙尘引起的PM10污染;持续较强的逆温层结以及近地面层弱南北风频繁交替出现造成沙尘在近地层往复、积聚,最终导致连续多日的空气污染。     

基于探空资料的江西典型冻雨天气过程垂直结构分析

利用探空资料和地面观测资料,对1959—2008年发生在江西的典型冻雨天气过程的大气垂直结构进行了分析。结果表明:冻雨的形成与大气饱和层、冰晶层、暖层、逆温层、冷层及特征物理量等因素的关系密切。对流层中下层的大气呈饱和或准饱和的逆温状态,饱和层顶位于暖层之上是冻雨发生的主要层结特征。冰晶层底部的平均高度为675hPa,平均厚度为2309m;暖层底部的平均高度为834hPa,平均温度为4.1℃,平均厚度为1765m;逆温层的平均温差为7.5℃;冷层的平均厚度为1668m,平均最低温度为-3.9℃;平均地面最低温度为-1.2℃;冻雨发生时大气整层比湿积分指数须增大为1000以上,干暖盖强度指数须降到-5℃以下。