复杂地形下黄河流域月平均气温分布式模拟

通过对常规气象站月平均气温资料的分析,发现在影响气温的众多因素中,海拔高度、太阳总辐射、地表长波有效辐射对气温具有显著影响.建立了月平均气温的物理经验统计模型,结合复杂地形下太阳总辐射分布式模拟结果,提出了依托常规地面气象观测资料实现复杂地形下月平均气温分布式模拟的方法,生成了黄河流域1km×1km分辨率月平均气温、月平均最高气温、月平均最低气温的空间分布图.分析表明,模拟结果能较好地反映气温的宏观分布趋势和局地分布特征.交叉验证结果表明,模型具有很好的稳定性,各月平均气温、月平均最高气温、月平均最低气温的模拟误差平均为0.19~0.35℃;加密站验证和个例年验证表明,模型具有良好的空间维和时间维模拟能力.提出的月平均气温分布式模型立足于常规地面气象观测资料,不依赖于山地野外考察资料,可以方便地在广大地区推广应用.     

起伏地形下黄河流域太阳散射辐射分布式模拟研究

针对天文、大气、宏观地理与局地地形等因子对起伏地形下太阳散射辐射的复杂作用,将影响起伏地形下太阳散射辐射的天空因素与地面因素分开处理.通过基于数字高程模型（DEM）数据的起伏地形下天文辐射模型和地形开阔度模型,综合考虑地面因素对散射辐射的影响;基于常规地面气象站观测资料建立的水平面散射辐射模型考虑天空因素对散射辐射的影响;依据各向异性散射机理,建立了起伏地形下太阳散射辐射分布式计算模型,探索出一条利用DEM数据和常规气象观测资料实现山区太阳散射辐射定量模拟的技术路线.以1 km×1 km分辨率的DEM数据作为地形的综合反映,实现了起伏地形下黄河流域1 km×1 km分辨率的太阳散射辐射分布式模拟.     

复杂地形三维直流电阻率有限元数值模拟

系统地论述了用有限单元法研究复杂地形条件下三维直流电阻率的正演计算技术.首先给出了三维构造中点源电场的边值问题以及相应的变分问题;然后利用有限单元法求解变分问题,采用四面体单元对研究区域进行剖分,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后,考虑到节约计算时间,利用对称超松弛顸条件共轭梯度迭代算法求解大型线性方程组,得到了各节点的电位值,进而计算出地表的视电阻率.通过理论模型的计算检验了算法的可行性之后,给出了几种常见纯地形异常的数值模拟结果和一个组合模型的计算结果,其研究工作为研究三维直流电阻率反演奠定了基础.     

起伏地形下天文辐射分布式估算模型

基于数字高程模型（DEM）,建立了起伏地形下天文辐射分布式估算模型.模型全面考虑了地形因子对天文辐射的影响,只需DEM数据作为输入项,适用于遥感图像处理、地理信息系统等数据处理平台.以1km×1km分辨率的DEM数据作为地形的综合反映,计算了我国全年各月天文辐射的空间分布.结果表明：我国年天文辐射总量有明显的纬向分布特点,随着纬度的降低,年天文辐射总量由北向南增加;由于受坡向、坡度和地形遮蔽因子影响,山区天文辐射表现出非地带性分布特征.本文所提供的我国天文辐射数据产品,可作为基础地理数据供相关研究应用.     

复杂地形条件下四面体剖分电阻率三维有限元数值模拟

地形起伏对电阻率影响非常大,复杂地形条件下的电阻率异常解释一直受到很多的关注.本文提出了一种四面体网格交叉剖分方法,剖分的网格相互交错,使剖分的网格具有多样性,能较好地模拟复杂地形情况下的地电模型.同时,从点电源场满足的方程出发,本文推导了三维复杂地形条件下有限元数值模拟算法,并编制了计算程序.算例表明,本方法行之有效,计算精度高.     

考虑地形起伏和障碍体破裂的汶川地震强地面运动数值模拟

2008年5月12日四川汶川地区发生M_W7.9地震,震中位置103.4°E,31.06°N.这次地震造成了以汶川、映秀为中心及其周边地域建筑物的严重破坏和人员的重大伤亡,且因为高山等地形复杂区域抢险救灾的艰巨性,为及时救援造成很大干扰.为更好理解地形因素对于强地面数值模拟结果的影响,建立了包含地形起伏影响及去除地形影响的两类模型.同时,依据震源破裂过程运动学反演结果,建立了包含障碍体破裂过程的震源滑动模型,实现断层分段、空间倾角以及滑移角的动态设定.基于动力学的地震动模拟方法,通过对地震波传播过程的数值计算和后处理分析,模拟由地震激发的区域强地面运动过程.结果显示：（1）强震动台站的断层距对地形效应具有放大或抑制作用,距离断层破裂带越近,地形效应越明显,反之,距离越远,则地形效应越微弱;（2）因为地形高差与障碍体的影响,地震造成的峰值可能出现在震中区域之外;（3）考虑地形影响模型的地表峰值速度（PGV）区域位于汶川与北川附近;而未考虑地形影响模型的PGV区域位于灌县—江油断层的后半段,处清平、安县附近;对汶川地震近实时强地面运动波场的模拟、峰值图谱的圈定及未来大地震强地面运动特征的预测都有重要指示意义.

     

三维复杂山谷地形SV波垂直输入地震反应分析

本文基于显式有限元法研究了地震波垂直入射时三维复杂山谷地形对地震地面运动的影响,在数值分析中应用了三维化二维的解法和黏弹性人工边界的处理方法,实现了地震波垂直输入下三维复杂场地地震动数值模拟,并验证了该方法的合理性.以四川桃坪地区一山谷地形作为研究对象,基于地表高程数据分别建立了二维和三维场地模型,对比研究表明：在复杂地形情况下考虑二、三维模型时具有明显差异,三维模型能更真实地反映地形变化对地震动的影响,复杂地形条件下有必要考虑三维实际场地模型.本文对边界自由场的处理方法也可用于处理三维复杂场地地震动斜入射问题,为三维复杂地形场地地震效应研究提供参考.

     

三维地形直流电阻率有限元法模拟

基于稳定电流场的基本方程、三维区域满足的边值问题以及相应的变分问题,研究了三维起伏地形条件下电阻率的有限单元数值模拟算法. 离散积分区域时,以三棱柱为最小研究单元,推导了含有地形特征信息的三线性插值型函数以及单元刚度矩阵. 采用变带宽、一维数组方式只存储稀疏刚度矩阵中非零元素,能够节约内存. 利用Cholesky分解法只分解一次大型稀疏矩阵,通过回代可以求出方程组的全部解,当求解有多个供电点的测深问题时可以缩短计算时间. 模型计算表明,在水平层状介质模型上,三维计算结果与解析解或二维数值解十分吻合,计算精度满足误差要求. 在二维山脊上的二极剖面或三维山谷上的中间梯度剖面上,其三维计算结果与相应模型的土槽实验结果或边界元法计算结果也非常接近.     

大地电磁场三维地形影响的数值模拟

用边界单元法解决了大地电磁(MT)三维地形影响的数值模拟问题.采用特殊的网格剖分,数值模拟方法,可以在微机上运行.在二维地形上,其结果与二维模拟结果基本吻合,对三维地形的模拟结果表明,三维模拟与二维模拟有较大差别.     

针对复杂地形的三种地震波走时算法及对比

复杂地形条件下地震波走时算法对于研究复杂地形地区的成像问题有着重要的意义.为了得到精度高且适应于复杂地形的走时算法,首先提出阶梯网格迎风差分法.然后将该方法与不等距网格有限差分法和混合网格线性插值法进行对比研究,得出如下结论:混合网格线性插值法的计算精度最高,但其计算效率最低;阶梯网格迎风差分法的计算精度最低,但其计算效率最高;不等距网格有限差分法的计算精度和计算效率均居中;而究竟选取哪种算法作为给定复杂地形模型的地震波走时算法,应该综合考虑地形的特点、所研究问题对计算精度及计算效率的要求等因素.最后通过一个计算实例验证了三种算法在面对复杂地形、近地表及地下复杂介质等复杂地质条件时均有很好的适应性和稳定性.     

Geochemical evidence for strong ongoing methanogenesis in Sanhu region of Qaidam Basin

Sanhu depression of Qaidam Basin is the largest biogenic gas production region in China.Headspace samples were collected from two wells in this region,and hydrogen and propylene compounds were detected in these samples with a certain concentration.The stable hydrogen isotope ratio of H 2 is relatively light (-700‰--820‰).The stable carbon isotope ratio of propylene ranges from -27‰ to -40‰,which coincides with the rule of change of the stable carbon isotope of kerogen at the corresponding horizon.The charac...     

The special sealing mechanism of caprock for Quaternary biogenetic gas in Sanhu area,Qaidam Basin,China

The Quaternary biogenetic gas reservoirs in the east of Qaidam Basin have many characteristics such as late forming time, shallow burial depth, low diagenetic grade, high porosity and high permeability and so on. It cannot be considered as caprock according to the traditional evaluation criterion. However, the large scale and high efficient biogenetic gas reservoirs of the Qaidam Basin are really formed under these kinds of caprocks, so it does have some specialty in its sealing mechanism. Aiming at the special sealing mechanism, some simulating experiments have been done. The research results show that the sealing ability of biogenetic gas caprock is related with water saturation, the caprock that is saturated with salt water can effectively block seepage and diffusion. Furthermore, the multiple reservoir-caprock groups have accumulated sealing effect, causing the formation of big gas fields. The evaluation method with traditional caprock parameters cannot be adopted in evaluating the study area.     

The special sealing mechanism of caprock for Quaternary biogenetic gas in Sanhu area,Qaidam Basin,China

The Quaternary biogenetic gas reservoirs in the east of Qaidam Basin have many characteristics such as late forming time,shallow burial depth,low diagenetic grade,high porosity and high permeability and so on.It cannot be considered as caprock according to the traditional evaluation criterion.However, the large scale and high efficient biogenetic gas reservoirs of the Qaidam Basin are really formed under these kinds of caprocks,so it does have some specialty in its sealing mechanism.Aiming at the special sealing mechanism,some simulating experiments have been done.The research results show that the sealing ability of biogenetic gas caprock is related with water saturation,the caprock that is saturated with salt water can effectively block seepage and diffusion.Furthermore,the multiple reservoir-caprock groups have accumulated sealing effect,causing the formation of big gas fields.The evaluation method with traditional caprock parameters cannot be adopted in evaluating the study area.     

The special sealing mechanism of caprock for Quaternary biogenetic gas in Sanhu area,Qaidam Basin,China

The Quaternary biogenetic gas reservoirs in the east of Qaidam Basin have many characteristics such as late forming time, shallow burial depth, low diagenetic grade, high porosity and high permeability and so on. It cannot be considered as caprock according to the traditional evaluation criterion. However, the large scale and high efficient biogenetic gas reservoirs of the Qaidam Basin are really formed under these kinds of caprocks, so it does have some specialty in its sealing mechanism. Aiming at the special sealing mechanism, some simulating experiments have been done. The research results show that the sealing ability of biogenetic gas caprock is related with water saturation, the caprock that is saturated with salt water can effectively block seepage and diffusion. Furthermore, the multiple reservoir-caprock groups have accumulated sealing effect, causing the formation of big gas fields. The evaluation method with traditional caprock parameters cannot be adopted in evaluating the study area.

     

Geochemical fractionation of rare earth elements in lacustrine deposits from Qaidam Basin

Based on the concentration and distribution pattern analysis of the rare earth elements (REEs) at the Shell Bar section from Qaidam Basin, we studied the geochemical fractionations of REE in paleolake deposits and their paleo-environmental significance. Our results show that the REE concentration in AS (acid soluble) and AR (acid residual) fractions are 20.9 μg/g and 95.4 μg/g (except element Y) individually, showing a strong REE differentiation between AS and AR fractions. However, the two types of fractions (AS & AR) have similar REE distribution patterns, which are slightly rich in light earth rare elements (LREEs), with slightly right-tilting and negative Eu anomaly. The LREE of AR is richer than that of AS. There were no significant correlations between the REE in AS, AR and other proxies. It indicated that the lacustrine deposition had different material sources and experienced varying geochemical procedures. Correlation analysis between the REE and the content of fine grain-size (<4 μm) of the sediments, Rb/Sr ratio and Mn concentration showed high correlation coefficients. Our results demonstrated that the REE in acid soluble and residual fractions bear different environmental significances, which are strongly dependent on local environment. The REE of acid soluble fraction is closely related to the paleoclimatic changes in the lake catchment and the evolutionary processes of lake itself. Our results suggest that δCe and (La/Yb)_n could be used as proxies of the reduction-oxidation conditions and furthermore the temperature change and the lake level fluctuations. Using these proxies, we reconstructed the paleoclimate and water level fluctuation history during the high lake level period lasting between 43.5 and 22.4 cal. ka BP.     

柴达木盆地东部中新生代两期逆冲断层作用的FT定年

柴达木盆地东部都兰一带,由于一系列北东倾的、向南西推覆的逆冲断层作用和温泉断裂的右行走滑作用,将柴达木盆地的变质基底和前中生界沉积基底翘起,从而使得柴达木盆地与共和盆地分隔开来. 逆冲断层中磷灰石裂变径迹（FT）测年结果反映了柴达木盆地基底断层作用的规律性. 从FT测年结果可以看出,FT年龄分为2组,也就是2个活动时期：第1期为108 Ma至61 Ma;第2期为26.6 Ma至17.8 Ma. 第1期反映的断裂活动具有明显的规律性：从柴达木盆地南缘的东昆仑开始,向柴北缘方向,逆冲推覆的断层作用时间逐渐变年轻,从东昆仑的108.0±9.6 Ma（柴达木南缘断裂）变为63.7±4.4 Ma（柴北缘断裂）,之后可能有小的跳动. 第2期,在原有的一些逆冲断层上形成了新的活动,或形成了一些新的逆冲断层,总体上具有无序或跳跃式变动的特点.     

可可西里岩石圈向北俯冲到柴达木地幔的地震学证据

利用中美德INDEPTH IV合作项目2007—2009年间布置于青藏高原中、北部140个宽频地震台站记录到的天然地震数据,经过接收函数成像处理,获得了3条穿过昆仑—阿尼玛卿缝合带清晰的壳幔结构图像.结果显示柴达木南缘莫霍面位于约50 km深度,羌塘地块、可可西里地块、东昆仑造山带莫霍面位于约65 km深度,昆仑—阿尼玛卿缝合带以北约50 km存在莫霍面深度突变.在可可西里和柴达木岩石圈地幔之间观测到北倾界面,这可能是可可西里岩石圈向北俯冲到柴达木地幔之下的证据.可可西里地块地壳内宽缓的负转换震相带是低速带的反映,其向北挤入到东昆仑山下发生挤压增厚,可能是东昆仑山隆升的原因;由于刚性柴达木岩石圈的阻挡,物质向东改向,则可能是该地区向东旋转的构造应力场产生的原因.本文研究结果不支持亚洲岩石圈地幔在东昆仑—柴达木交界处向南俯冲,据此,我们提出了新的东昆仑造山模式.

     

Geochemical evidence for non-marine depositional environment of foraminiferal fossils from the Nihewan Basin, China

From the first finding in 1970s, the findings of foraminiferal fossil assemblages in inland basins have been reported from time to time, especially in recent years. The debates on the depositional environment of foraminiferal fossils have become the hot spot of researches again in China. Based on the researches of trace element geochemistry and electron scanning microscope of shells of Quaternary foraminiferal fossils from the Xiaodukou section in the inland Nihewan basin, the original ⁸⁷Sr/⁸⁶Sr and other geochemical information of shells were believed to be preserved well and could be used to rebuild the geochemistry of contemporary waters where foraminifera deposited, although there existed some effects of burial diagenesis on the geochemistry of shells to a certain extent. The ⁸⁷Sr/⁸⁶Sr ratios of well-preserved Xiaodukou foraminiferal shells were measured, giving a range of 0.711190±25–0.712018±14, apparently higher than the value of contemporary seawater (0.709087–0.709147) and similar to that of the Sanggan River, proving that it represented the value of the ancient lacustrine water. The hyperbolic mixing models of ⁸⁷Sr/⁸⁶Sr-palaeosalinity and ⁸⁷Sr/⁸⁶Sr-Sr/Ca indicated that the contemporary waters where Xiaodukou foraminifera inhabited was an inland lake and there was no seawater input to the depositional environment.     

Geochemical evidence for sundering of the West Mariana arc in miocene ash from the Parece Vela Basin

Glass and mineral fragments from discrete volcanic ash layers were sampled from DSDP/IPOD Site 450 in the Parece Vela Basin, Philippine Sea and analyzed by electron microprobe. The ashes are interpreted as eruptive products of the adjacent West Mariana arc system between 25 and 14 Ma B.P., and have compositions between basaltic andesite and rhyolite, and rarely, boninite. ‘Continuous’ chemical trends appear to reflect mixing of mafic and silicic magmas. ‘Discontinuous’ trends between these end-members are relatively few, and are consistent with ‘liquid lines’ produced by fractional crystallization. Andesitic tephra become progressively richer in MgO and CaO through the middle Miocene, while boninite appears towards the end of the sequence, between 14 and 15 Ma B.P. Coeval rhyolitic glasses become richer in K₂O and Na₂O, with maximum concentrations at about 15 Ma B.P. Chronologic changes in fractionation type and composition of parent magmas are interpreted to reflect the subaerial volcanic evolution of the West Mariana arc. The appearance of boninite is believed to signal early stages of arc sundering, and corresponds temporally with regional uplift of the sea floor above the carbonate compensation depth, precursor to a new pulse of back-arc spreading.     

横波岩石物理地震随机模拟预测天然气藏分布——以柴达木盆地三湖地区北斜坡为例

柴达木盆地三湖地区北斜坡第四系生物气区是中国最大的生物气区,上第三系为天然气勘探潜力区.第四系地层为弱成岩的疏松砂泥岩薄互层.在以2D地震为主的三湖地区北斜坡预测岩性气藏,利用叠后波阻抗反演及叠前AVO反演预测含气分布难度很大.而综合地质、测井、地震、生产测试等多学科资料,通过对地震道的岩石物理参数进行随机模拟,获得对含气异常比较敏感的地下岩石物理属性体,并分析最敏感的岩石物理属性产气、产水的阈值区间以建立气水识别模式,再利用气水识别模式对岩石物理属性三维模型进行流体置换,最终得到可能的关于生物气藏的三维空间分布特征.利用本文的方法成功预测和落实了上第三系新领域的风险探井.钻探实效表明,该方法适用性强,可以在很大程度上降低天然气藏预测的多解性和不确定性,为进一步优选有利钻探区带、提出可供钻探的目标与井位、缓解该区储量不足的矛盾提供有力的技术支撑.