3D高阶抛物Radon变换在不规则地震数据保幅重建中的应用

3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换（NHOPRT）地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而获得最接近完整规则数据的高阶抛物Radon变换域系数.在用SVD求解反演方程过程中,利用高阶抛物Radon变换算子在频率域为指数函数,具有线性可分解特性,将二维空间的高阶抛物Radon变换算子分解为两个独立的一维空间变换算子,减小了变换算子的矩阵大小,从而很大程度地提高了计算效率.理论模型和实际地震数据重建测试证明了本文方法的有效性以及实用性.     

酒泉、嘉峪关地域城镇系统的空间规模结构研究

受干旱区自然地理、人口发展与分布等特殊条件的制约,酒泉、嘉峪关地域城镇系统的空间规模结构呈现“中小市镇、多中心”的特征。针对这种状况,空间规模结构的研究和规划,应注重建构多中心的地域城镇系统,带动地域整体发展要促进相邻城市间、城市与境内或周边单位间的空间关联及与近郊集镇的融合,推动城乡一体化应突出酒泉-嘉峪关复合城市的中心地位,强化敦煌市的规模和职能,保证玉门市区持续发展并尽快发展玉门镇,积极发展其他小城镇。     

非均匀介质中交错网格高阶有限差分数值模拟

地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。     

节气气候变化与上消化道出血的关系研究

通过对2000-2005年龙岩市第一医院上消化道出血（UGB）患者入院时的节气归属及当时气象资料的收集,分析不同节气气候变化与UGB的关系,运用最优子集回归方法建立关系模式进行综合分析,得出影响上消化道出血发病趋势的关键气象要素是节气最高气温、最低气温、平均水汽压;通过节气气候因素对UGB影响的认识,提出相应的预防保健措施。     

1961~2010年黄河中下游地区24节气气候变化特征分析

24节气是我国古代劳动人民的独特创造。全面掌握"24节气"的气候变化规律,不但有利于指导农事生产,提高气象服务质量和水平,而且在为人类预防和治疗疾病方面也有重要意义。在全球变暖的气候背景下,统计分析了1961~2010年黄河中下游地区24节气的气温、湿度、风速等6个气象要素的变化特征,得到以下结论:黄河中下游地区随节气变换气候变化显著,大暑、小暑节气高温高湿,小寒、大寒节气寒冷干燥,清明节气寒温反复大风将至,霜降节气天气渐凉秋燥加剧等。50年内,春季型节气(平均、最高、最低)气温显著升高,冬季型节气最低气温升高显著。气压随节气变化特征与气温大致相反,夏、秋季节气有升压趋势。相对湿度与降水均呈减少趋势,以秋季型节气减小趋势最明显。春季风速最大,夏、秋季风速最小,所有节气风速均呈减小趋势,冬夏季节气日照时间呈缩短趋势。     

明嘉靖六年《大统历》历书的气朔推步精度

历书是历法的民用本.考察了明嘉靖六年(1527年)<大统历>历书,根据历书所颁朔日和节气时刻,以及先期恢复的<授时历>和<大统历>推步术,并利用现代天文历算结果,可考校得出:(1)历书所载朔日和节气日期、时刻确由<大统历>推得;(2)<大统历>与<授时历>推步结果存在差异,推朔平均差33.35 min或34.93 min,推气平均差86.40 min,且两者均有差1日现象;(3)<大统历>与现代推算合朔时刻的平均误差为24.50 min或44.03 min,与现代推算节气时刻的差平均为1.25 d,若扣除太阳中心差则节气推步的误差将减少到86.40 min.     

On the recent warming of the southeastern Bering Sea shelf

During the last decade, the southeastern Bering Sea shelf has undergone a warming of 3 °C that is closely associated with a marked decrease of sea ice over the area. This shift in the physical environment of the shelf can be attributed to a combination of mechanisms, including the presence over the eastern Bering Sea shelf of a relatively mild air mass during the winter, especially from 2000 to 2005; a shorter ice season caused by a later fall transition and/or an earlier spring transition; increased flow through Unimak Pass during winter, which introduces warm Gulf of Alaska water onto the southeastern shelf; and the feedback mechanism whereby warmer ocean temperatures during the summer delay the southward advection of sea ice during winter. While the relative importance of these four mechanisms is difficult to quantify, it is evident that for sea ice to form, cold arctic winds must cool the water column. Sea ice is then formed in the polynyas during periods of cold north winds, and this ice is advected southward over the eastern shelf. The other three mechanisms can modify ice formation and melt, and hence its extent. In combination, these four mechanisms have served to temporally and spatially limit ice during the 5-year period (2001–2005). Warming of the eastern Bering Sea shelf could have profound influences on the ecosystem of the Bering Sea—from modification of the timing of the spring phytoplankton bloom to the northward advance of subarctic species and the northward retreat of arctic species.     

天然地震数值模拟与实例对比

为了解决长时间正演模拟的不稳定性问题, 实现天然地震波场长时间数值模拟.在高阶有限差分数值模拟的基础上, 给出了多轴向完全匹配吸收边界(M-PML)二维划分方式, 讨论了M-PML吸收边界转换系数P的取值对其吸收能力与截断误差对数值模拟有效信号的影响.实现了二维弹性波高阶有限差分的长时间数值模拟.采用2013年12月16日湖北省巴东县地震及余震信息, 对该地震数据进行了数值模拟, 并与武汉、秭归两台站数据进行了分析和讨论, 验证该正演模拟方法的正确性和有效性.为天然地震波场传播规律、成像和震相识别及震源定位等研究提供了基础.      

2002年12月25日乌恰5.8级地震的长、中、短期跟踪预测和反思

在乌恰地区地震活动呈现长期平静异常的背景下,通过多种方法对2002年12月25日乌恰5．8级地震进行了1-3年的中长期预测。应用“模糊地震活动度”和“空区参数法”两种方法进行了1年和半年的中期预测,结合对该地区小震活动动态的分析与震例分析,尤其是预警地震的分析判断,提出了3个月的短期预测意见,较好地预测了这次地震。最后对预测情况进行了正反两个方面的总结和反思,提出了七点看法。     

An Adaptive Nonhydrostatic Atmospheric Dynamical Core Using a Multi-Moment Constrained Finite Volume Method

An adaptive 2 D nonhydrostatic dynamical core is proposed by using the multi-moment constrained finite-volume(MCV) scheme and the Berger-Oliger adaptive mesh refinement(AMR) algorithm. The MCV scheme takes several pointwise values within each computational cell as the predicted variables to build high-order schemes based on single-cell reconstruction. Two types of moments, such as the volume-integrated average(VIA) and point value(PV), are defined as constraint conditions to derive the updating formulations of the unknowns, and the constraint condition on VIA guarantees the rigorous conservation of the proposed model. In this study, the MCV scheme is implemented on a height-based, terrainfollowing grid with variable resolution to solve the nonhydrostatic governing equations of atmospheric dynamics. The AMR grid of Berger-Oliger consists of several groups of blocks with different resolutions, where the MCV model developed on a fixed structured mesh can be used directly. Numerical formulations are designed to implement the coarsefine interpolation and the flux correction for properly exchanging the solution information among different blocks. Widely used benchmark tests are carried out to evaluate the proposed model. The numerical experiments on uniform and AMR grids indicate that the adaptive model has promising potential for improving computational efficiency without losing accuracy.