基于新的差分结构的时-空域高阶有限差分波动方程数值模拟方法

传统的高阶有限差分波动方程数值模拟方法采用高阶差分算子近似空间偏导数,能有效抑制空间频散.然而,传统的有限差分法仅采用二阶差分算子近似时间偏导数,这使得地震波场沿时间外推的精度较低.当采用较大的时间采样间隔,传统的有限差分法模拟波场会出现明显的时间频散,甚至不稳定.本文基于新的差分结构和中心网格剖分,发展了一种空间任意偶数阶精度、时间四阶和六阶精度的时空域有限差分方法.基于对离散后的频散关系进行泰勒展开,本文推导了时空域高阶有限差分算子的差分系数.相速度分析表明时间四阶、六阶精度的差分方法能显著地减小传统时间二阶精度差分方法的时间频散.在相同的精度下与传统差分法比较,本文发展的时间四阶、六阶有限差分方法的计算效率比传统方法高.均匀和非匀均介质中的波场数值模拟实验进一步证实本文研究的时空高阶有限差分方法的优越性.     

贵州岩屋坪汞矿床有机流体成矿的包裹体证据

岩屋坪汞矿床是铜(仁)凤(凰)汞矿带内重要的汞矿床之一。对该矿床中的流体包裹体进行了研究,结果表明:矿床内存在多种类型包裹体,包括H2O包裹体、烃H2O包裹体、烃类包裹体和固体辰砂包裹体。气液H2O包裹体均一温度变化于82.4~176.5 ℃,流体盐度为13.40％~21.75％NaCl,密度主要集中在1.05~1.10 g/cm3,成矿流体属于一种高盐度、高密度的油田卤水。矿石中沥青和辰砂晶体常分布在同一微裂隙中,辰砂晶体中常捕获有沥青包裹体,沥青包裹辰砂颗粒,烃类包裹体和固体辰砂包裹体密切共生,说明有机质与汞矿化的关系极为紧密,为研究有机质在汞的成矿过程中的作用提供了重要依据。     

甘肃龙首山碱交代型铀矿床绿泥石特征及意义

绿泥石化是龙首山铀矿床重要的蚀变类型之一。通过对龙首山碱交代型铀矿床的绿泥石等蚀变矿物进行的岩相学和电子探针成分分析研究,确定了龙首山地区绿泥石的化学类型主要为铁镁绿泥石,少数为蠕绿泥石。依据绿泥石成因或与共生矿物的关系,绿泥石可被划分为黑云母蚀变型、长石蚀变型、沥青铀矿共生型和副矿物共生型等4种类型。泥质岩是本区绿泥石的主要原岩类型,是多期次地质作用形成的产物。研究认为,龙首山地区碱交代型铀矿床的成矿过程可表述为矿前期在相对较高温度的热液流体作用下,黑云母发生绿泥石化蚀变,随后热液继续交代长石,形成长石蚀变型绿泥石,进而在成矿期热液温度相对较低的条件下形成与沥青铀矿紧密共生的绿泥石。绿泥石在铀成矿过程中不但活化了花岗岩里的铀,而且还给铀矿化供应了相对良好的积淀环境。     

塔北隆起X型走滑断裂成因机制的新解释

塔里木盆地塔北隆起发育两组呈小角度相交(40°)的透入性X型走滑断裂,分别沿着NNE走向和NNW延伸。在对塔北哈拉哈塘地区三维地震资料解释的基础上,对走滑断层的几何展布特征以及断层的剖面变形特征进行研究;同时重点解析了RP6断裂和HA13断裂,分析比较NNW向与NNE向断层的变形及发育特征差异;结合盆地重磁资料以及周缘造山带的活动特征,对塔北隆起小角度的X型走滑断层的发育机制以及演化进行了分析。研究表明,塔北隆起走滑断层在垂向上具有明显的分层变形特征,分为三个构造层:震旦系-中寒武统下构造层(TH₃界面以下)、上寒武统-中奥陶统中构造层(TH₃-TO₃t界面)和上奥陶统-石炭系上构造层(TO₃t-TP界面)。断层在下构造层和中构造层中整体处于压扭环境,多发育正花状构造;上构造层中断层主要发育负花状构造或正断层,整体处于张扭环境。两组断裂比较,NNW向断裂活动性强,在各构造层中均有显著的断裂特征发育,垂向连通性强,发育先存基底断裂,而NE向断层主要发育在中构造层,在下构造层和上构造层中断层发育不明显。活动性分析表明,断层的形成与演化具有多期性,走滑断层的形成经历了三期主要活动:中寒武世末、中晚奥陶世和志留纪-石炭纪。塔北隆起X型走滑断裂的形成受到了NNW向基底断裂和薄弱带的控制,NNW向先存基底断裂带或薄弱带优先发育走滑断裂,基底断裂与主挤压应力方向的夹角小于45°-Φ/2,NNE断层的发育受NNW向先存断裂限制,最终形成小角度相交的X型断裂。     

湿斜压大气中暴雨中尺度系统发展的一种可能机制

根据暴雨中尺度系统发生时的大气运动特征,在绝热、无摩擦条件下,从湿斜压原始方程出发,考虑大气运动的涡度及散度演化不仅受动力场的制约,而且还受到热力场的约束,采用不同于传统研究涡度及散度方程的分析方法,导出了新形式的涡度及散度方程。在此基础上,分析了湿斜压大气中对流层中低层气流旋转与辐合持续增长的动力特征,初步揭示了湿斜压大气中暴雨中尺度系统发展的动力机制。     

岩爆及采矿诱发岩体失稳破坏过程数值模拟新方法

岩石介质的宏观非线性主要是由非均质性和各向异性造成的,应用新的数值计算软件RFPA(2D)——基于岩石宏观非线性行为可能是由具弹脆性特征的细观微元体不断破裂造成的,造成微元体不断破裂(并非某一时刻同时破裂)的原因是微元体材料性质(微元体强度、弹性模量和泊松比等参数)的非均匀性,是一个数学上相对简单但能充分考虑岩石介质复杂性的方法,对采动引起岩体失稳破坏的全过程进行了数值模拟研究,表明应用RFPA(2D)是切实可行的。      

Changjiang River sediment delivering into the sea in response to water storage of Sanxia Reservoir in 2003

1 IntroductionHuman activitieshavebeen stronglyaffecting theworld river sedimentsupply to the oceans since the19thor20thcentury(StanleyandW arne,1993;Fanos,1995;Sm ithandW inkley,1996;Stanley,1996;Yang,Shi,etal.,2004;Yang,Zhao,etal.,2002).Thesedi-m enttra…     

星座美洲海鞘次级代谢产物的研究

从采自青岛近海的星座美洲海鞘(Amaroucium constellatum)70%乙醇提取物中分离得到9个次级代谢产物.经波谱(MS、NMR)分析和文献比对,确定了它们的化学结构分别为尿嘧啶(1)、胆甾醇(2)、十九碳甘油醚(3)、胸腺嘧啶(4)、反-4-对羟基苯基-3-丁烯-2-酮(5)、反-4-苯基-3-丁烯-2-酮(6)、对胺基苯乙酸(7)、胸腺嘧啶脱氧核苷(8)和N-甲基尿嘧啶核苷(9).其中,化合物2～3,5～7,9均为首次从该海鞘中分离得到.     

苏州城区场地等效剪切波速计算深度取值探讨

鉴于场地等效剪切波速的计算深度是否需要由地表以下20 m增至30 m的争论尚未有定论,依据苏州城区场地143组实测剪切波速资料的统计分析,对GB50011-2010《建筑抗震设计规范》场地类别分类标准与欧美抗震设计规范场地分类标准进行了对比,比较了苏州城区场地等效剪切波速计算深度取为20 m和30 m的差异性,结果表明:苏州城区场地剪切波速在地表以下20 ～ 60m范围的变异性较大,采用分段函数描述剪切波速随土层深度的变化关系是合适的;Ⅲ类场地等效剪切波速大小的分布偏向Ⅲ类场地类别的下界限值,当计算深度由20 m增至30 m时,场地等效剪切波速的均值增大16％,如果直接沿用欧美规范的场地分类界限值,将会整体提高苏州城区场地类别的划分标准,苏州城区Ⅲ类和Ⅳ类场地的等效剪切波速分界值取为170 m/s是适宜的.     

2008年汶川8.0级地震孕震机理研究

用三维流变非连续变形与有限元相结合（DDA+FEM）的方法,在青藏高原及邻近地区三维构造块体相互制约的大背景中,考虑了龙门山断裂带东西两侧地势、地壳厚度和分层的明显变化,及断裂带东侧四川盆地及鄂尔多斯块体坚硬地壳阻挡的影响,通过用GPS资料做位移速率边界约束和震源机制约束,计算得到研究区的速度场和应力场与该地区GPS测量结果和震源机制分布结果基本一致.在此基础上,模拟计算现今构造块体边界断层上表征剪应力及法向应力等综合影响的危险度分布.结果表明,上、中地壳层危险度分布中危险度较高的地段多数与近几十年来发生的七级以上大震区域基本一致.包括2008年汶川8.0级等大震的发震断层.通过分别对龙门山断裂带东西两侧的两种不同构造格局进行试算表明,龙门山断裂带东西两侧地势、地壳厚度、分层与物性明显变化对汶川大震的孕育发生均起了关键性作用.计算得到的应变率强度分布图可见,高原东部整个边缘地带均接近应变率强度的陡变带.其中以龙门山断裂带上的陡变最为明显,西侧应变率强度是东侧的近4倍,而且断裂带东侧应变率强度等值线衰减比西侧快.反映了汶川大震逆冲型发震断层地区独特的特征.此外,由计算得到的应变能密度分布图可见,龙门山断裂带在上、中地壳层中均位于宽度相同、其走向与龙门山断裂带走向一致的高应变能密度带中,在上地壳层这个带的东西两侧则是应变能密度较低的地区,而在中地壳层,其强度在断裂带东侧逐渐向东衰减,西侧应变能密度高,而东侧应变能密度较低.表明在印度板块强烈推挤作用和高原各构造块体相互制约及龙门山断裂带东西两侧特殊构造环境中,高原地壳物质向东水平运动,受到龙门山断裂带东侧介质刚性强度较大的四川盆地阻挡,使得汶川大震发震断层在大震前已积累了相当水平的应变能,并同时处于力学上的不稳定状态.