一种利用Nyström离散与FFT快速褶积的散射地震波并行计算方法

利用数值方法解Lippermann-Schwinger (L-S)方程的主要困难在于系数矩阵存储和线性方程组求解.这主要是因为L-S方程的积分部分是一个空间褶积,在离散后将导致一个满秩矩阵,进而形成一个大型或超大型代数方程组.因此,在利用L-S解决地震波散射问题时,一般是利用散射级数法而非数值方法.然而,散射级数法的计算精度和收敛性强烈地依赖于速度扰动的强度,而克服这种依赖性的一个可能的途径就是对现有的数值方法进行改进或是建立新的数值求解方案.在这种思想指导下,首先对L-S方程进行改写,得到一个与原L-S方程等价的积分方程(等价L-S方程).然后,对等价L-S方程进行逐点归一化处理,并利用Nyström法对经归一化处理的等价L-S方程(归一化等价L-S方程)进行离散,并用FFT计算空间褶积.之所以这样选择是由于归一化等价L-S方程经Nyström法离散生成的系数阵为一个Toeplitz阵,可利用其Toeplitz性质降低存储空间;而FFT可以将矩矢空间褶积转化为乘积,且积分核部分只要计算一次即可.进一步,为节约正演计算时间,设计了进程级和线程级相结合的MPI+OpenMP并行模式.数值试验表明,与传统的积分方程数值算法相比,利用等价L-S方程、Nyström离散和FFT快速褶积的计算方案可极大地降低存储需求,进而在保证精度的同时提高计算效率.

     

内蒙古北山造山带时空结构与古亚洲洋演化

内蒙古北山造山带自北向南发育红石山－百合山、月牙山－洗肠井和帐房山－玉石山3条蛇绿构造混杂岩带。其中北部的2条蛇绿岩带揭示了北山造山带两阶段演化的历史：月牙山－洗肠井蛇绿混杂岩带洋壳形成于530~520 Ma,沿该带多处保存较完好的蛇绿岩洋壳残块（洋壳结构具向北变新特征）,与北侧的早古生代公婆泉岩浆弧（由南向北岛弧成熟度变高）共同指示了北山洋向北俯冲消减的过程,即490 Ma初始俯冲,450~440 Ma为俯冲峰期,430~420 Ma为同碰撞阶段,400 Ma的双峰式岩浆岩组合指示了北山洋的消亡和后造山伸展的过程;红石山－百合山蛇绿混杂岩带是发育在雀儿山－圆包山岛弧基础上的SSZ型蛇绿岩,弧后开裂洋壳的形成与南侧最早发育的岛弧岩浆作用年龄接近（340~320 Ma）,310~290 Ma俯冲峰期造成南侧白山岩浆弧大量的岩浆活动,早二叠世末期（275 Ma）的辉长岩和花岗岩侵位及早—中二叠世双堡塘组下部的角度不整合均反映了红石山洋盆的闭合。前人所划"石板井－小黄山蛇绿岩带"实为一条早古生代发育的深大断裂,沿带发育中基性侵入体及少量超基性岩,后期（志留纪末）叠加有较强的韧性剪切变形。中生代以来的走滑作用和逆冲推覆构造改造了古生代的构造格架,使红石山-百合山蛇绿混杂岩带向北左行切错了十余千米,北山南部的中—新元古界推覆至下古生界之上。对内蒙古北山造山带时空结构的厘定,有助于中亚造山带造山作用过程的理解及其对古生代地壳增生的深入研究,也对银额盆地晚古生代新层系油气资源勘查起到基础支撑。     

灰色GM(1,1)模型最优维数的研究

灰色GM(1,1)模型的维数对模型的预测精度有很大的影响,本文通过自编的FORTRAN程序,在合适的范围内对GM(1,1)模型的维数进行自动优化,通过对宁波市地面沉降数据的计算分析,得到了C和P的定性关系,确定了优化的评价标准,维数最优化后的模型能反映预测数列变化幅度的大小.     

从冷科1井下部地层的归属探讨柴达木盆地成盆时间

冷科1井是近年来柴达木盆地北缘最重要的一口科学探索井.因该井首次发现了巨厚的新地层,所以对该井的4 666~5 200 m地层划分存在较大争议.前人主要依据孢粉资料把冷科1井4 666~5 200 m层段归属为下侏罗统,也有人依据从该段发现的海相化石碎片认为该段为石炭系.本文从岩性组合、古生物特征、地震剖面、测井电性及水介质特征等多方面对该层段进行了分析,否定了该段为下侏罗统和石炭系的可能,认为该段为上三叠统.该研究改变了过去柴达木盆地始于早侏罗世的观点,指出柴达木盆地起始于晚三叠世零星分布的断陷沉积.     

Effects of temporal scales and space mismatches on the TRMM 3B42 v7 precipitation product in a remote mountainous area

The accurate measurement of precipitation is essential to understanding regional hydrological processes and hydrological cycling. Quantification of precipitation over remote regions such as the Tibetan Plateau is highly unreliable because of the scarcity of rain gauges. The objective of this study is to evaluate the performance of the satellite precipitation product of tropical rainfall measuring mission (TRMM) 3B42 v7 at daily, weekly, monthly, and seasonal scales. Comparison between TRMM grid precipitation and point‐based rain gauge precipitation was conducted using nearest neighbour and bilinear weighted interpolation methods. The results showed that the TRMM product could not capture daily precipitation well due to some rainfall events being missed at short time scales but provided reasonably good precipitation data at weekly, monthly, and seasonal scales. TRMM tended to underestimate the precipitation of small rainfall events (less than 1 mm/day), while it overestimated the precipitation of large rainfall events (greater than 20 mm/day). Consequently, TRMM showed better performance in the summer monsoon season than in the winter season. Through comparison, it was also found that the bilinear weighted interpolation method performs better than the nearest neighbour method in TRMM precipitation extraction.     

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.     

Suspended sediment and total dissolved solid yield patterns at the headwaters of Urumqi River,northwestern China: a comparison between glacial and non‐glacial catchments

In order to understand the differences in the suspended sediment and total dissolved solid (TDS) yield patterns between the glacial and non‐glacial catchments at the headwaters of Urumqi River, northwestern China, water samples were collected from a glacier catchment and an empty cirque catchment within the region, during three melting seasons from 2006 to 2008. These samples were analyzed to estimate suspended sediment and TDS concentrations, fluxes and erosion rates in the two adjoining catchments. There were remarked differences in suspended sediment and TDS yield patterns between the two catchments. Suspended sediment concentrations were controlled mainly by the sediment source, whereas TDS concentrations were primarily related to the hydrologic interaction with soil minerals. Generally, the glacial catchment had much higher suspended sediment and TDS yields, together with higher denudation rates, than the non‐glacial catchment. Overall, glacial catchment was mainly dominated by physical denudation process, whereas the non‐glacial catchment was jointly influenced by physical and chemical denudation processes. The observed differences in material delivery patterns were mainly controlled by the runoff source and the glacial processes. The melting periods of glacier and snow were typically the most important time for the suspended sediment and TDS yields. Meanwhile, episodic precipitation events could generate disproportionately large yields. Subglacial hydrology dynamics, glaciers pluck and grind processes could affect erodibility, and the large quantities of dust stored on the glacier surface provided additional sources for suspended sediment transport in the glacial catchment. These mechanisms imply that, in response to climate change, the catchment behaviour will be modified significantly in this region, in terms of material flux. Copyright © 2013 John Wiley & Sons, Ltd.     

基于空间技术实测数据的板块构造运动研究

将ITRF2008解算结果 GPS、SLR、VLBI实测站速度信息和北美CORS网络、中国陆态网络、欧洲CORS网络站速度信息进行基准统一,联合解算IERS组织划定的全球14个主要大陆板块的欧拉矢量,对原有的板块运动参数进行精化;提出了基于最小二乘法的海量台站选择方法,实验结果表明,该台站选取方法能够有效提高海量台站的筛选效率和可靠性。稳定台站的数量越多,分布越均匀,越能够体现板块运动的现势性。     

高分二号卫星影像融合方法探析

图像融合技术是卫星影像处理过程中的一个重要步骤,它可以将全色影像的高空间分辨率特点与多光谱影像的高光谱分辨率特点有机地结合起来,使得融合后的影像更利于后续解译。由于融合处理方法的原理以及获取影像的传感器特性不同,不同的卫星影像适合于不同的融合方法。针对高分二号卫星影像的特点,本文采用Gram-Schmidt融合法、FIHS变换方法、HPF方法、Pansharp方法进行试验。实验结果表明:HPF方法所得融合结果在光谱特性保持方法性能最优,而Pansharp方法融合结果在信息量的保持以及清晰度方面效果更优,因此,在实际应用中,应根据具体需求进行方法选择。     

基于等高线隔断线确定DEM适宜分辨率的方法研究

分辨率是等高线数据插值生成格网DEM数据过程中需要考虑的一个重要参数,适宜的分辨率既能将等高线高程信息尽可能保留,又不会因插值造成数据冗余及误差。基于对格网DEM生成方式的理解,提出了等高线隔断线的相关概念,分析了基于地形图矢量化的等高线数据插值生成格网DEM数据,由不同长度等高线隔断线频率统计来确定格网适宜分辨率参考值的方法。同时,不同长度等高线隔断线的长度与其频率相乘得到标准长度,以及等高线隔断线中心点生成的密度分析图,均可在频率特征决定适宜分辨率参考值时起到一定的参考作用。通过样区实验验证,实现简便,可度量性强,与其他相关研究具有一致性,对确定适宜格网分辨率参考值具有一定的借鉴意义。