Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

A simple hypoplastic model for normally consolidated clay

The paper presents a simple constitutive model for normally consolidated clay. A mathematical formulation, using a single tensor-valued function to define the incrementally nonlinear stress–strain relation, is proposed based on the basic concept of hypoplasticity. The structure of the tensor-valued function is determined in the light of the response envelope. Particular attention is paid towards incorporating the critical state and to the capability for capturing undrained behaviour of clayey soils. With five material parameters that can be determined easily from isotropic consolidation and triaxial compression tests, the model is shown to provide good predictions for the response of normally consolidated clay along various stress paths, including drained true triaxial tests and undrained shear tests.     

库车坳陷侏罗系煤成气动力学模拟研究

应用黄金管—高压釜封闭体系热模拟实验与GC、GC-IRMS分析技术,结合KINETICS专用软件,对库车坳陷侏罗系煤成气进行了动力学模拟研究。库车坳陷侏罗系煤具有高的产气性,在高演化阶段主要产甲烷气;侏罗系煤热解气甲烷碳同位素为-36‰～-25‰,乙烷碳同位素为-28‰～-16‰;甲烷、C2-C5气态烃的生成活化能分别为(47～64k)calm/ol、(55～72k)calm/ol,频率因子各为5.265×1013s-1、5.388×1018s-1。在此基础上,进一步探讨了克拉2气田天然气的成因。研究认为,克拉2气田天然气属阶段捕获的煤成气,主要聚集了5～1Ma时期的天然气,其成熟度Ro分布范围为1.3%～2.5%。     

张家口地下热水资源成因探讨

张家口地热资源丰富，通过对区内地热成矿机制，成矿规律的分析认为，以祁吕系为主的多构造体系复合部位控制了本区地下热水的形成与分布，本区地下热水是大气降水经深循环而成，经计算深循环的深度达为2600m，张家口南部地区有很多个热水异常点，是具有远景意义的地热区。     

邻苯二酚螯合树脂的合成及水样中铝的测定

合成并表征了新螯合树脂———邻苯二酚螯合树脂(XAD-2-Catechol),研究了XAD-2-Catechol吸附铝的特性和茜素红-铝的显色反应,在pH4.5的HAc-NaAc缓冲介质中,茜素红和铝(Ⅲ)反应生成红色络合物,λm ax=500 nm,铝的含量在0-50μg/25 mL内符合比耳定律。建立了邻苯二酚螯合树脂分离/富集-茜素红分光光度法测定天然水样中铝的新方法,对水样中铝形态进行测定,结果满意。     

酶法催化餐饮废油制备生物柴油的研究

以废物资源转化为宗旨,利用固定化脂肪酶催化餐饮废油与乙醇反应制备生物柴油,通过实验获得最佳酯化反应条件:反应温度47℃、有机溶剂为正己烷、醇油比3∶1,5次投加乙醇,酶用量为0.3g,反应时间32h时,生物柴油产率可达81%。     

REE/trace element characteristics of sandstone-type uranium deposits in the Ordos Basin

1 Introduction The Ordos Basin is the second largest sedimentary basin in China. During the last 10 years, a great progress has been achieved in the aspects of tectonic evolution, dynamics process, inner and outer geological processes during Mesozoic-Cen…     

Detrital Zircon of 4100 Ma in Quartzite in Burang, Tibet

A detrital zircon aged 4.1 Ga is discovered by the SHRIMP U-Pb method in a quartzite in Burang County, western Tibet. This is presently the oldest single-grain detrital zircon in China. The Th-U ratios of the two testing points of the 〉4.0 Ga zircon are between 0.76 and 0.86, indicating their magmatic origin. This discovery has offered an important age for investigating the geological evolution of the Qinghai-Tibet Plateau.     

重力斜导数方法在时间域中的理论模型与研究

Miller和Singh在1994年首次提出了斜导数(Tilt梯度)的定义,J.Derek Fairhead和Chris M.Green又于2004年在斜导数的基础上提出了Tilt梯度水平导数的概念。国内学者已对斜导数及斜导数水平导数展开研究,设定了多体模型、无限延伸倾斜脉模型、有限延伸倾斜台阶模型及埋深相同无限延伸垂直厚板及薄板模型,得出一些结论。针对一些地质体产状的判断,提出非绝对值斜导数水平导数的新概念,并应用于理论模型中验证,得出正确的结论。在上述理论模型基础上补充了断层模型和无限延伸倾斜台阶模型,应用非绝对值水平导数的理论解释地质体的产状,并对上述模型提出新的观点和看法。     

VSP资料约束地层反滤波研究

浅地表的低速带和地层滤波作用导致地震信号的能量衰减,地震子波的高频快速衰减,速度频散的子波变形,对地震底层分析带来了严重的影响。地层反滤波算法就是为了消除地层对子波的影响,利用VSP中的直达下行波推算出每个地层滤波算子,计算出各个地层滤波算子的反算子,然后利用这些反算子恢复出地表记录中损失掉的高频成分,减弱地层滤波对波形的影响,增加地表记录的频率成分和波形的准确度。分析实际地震资料表明,经过地层反滤波处理后的记录频率成分明显提升,地质构造更加准确、清晰。