基于交错网格有限差分算法的页岩声波各向异性校正方法研究

页岩具有层理发育的特征, 这会引起强烈的声波各向异性, 导致直井与水平井声波测井数据之间差异明显, 因此在水平井储层参数计算中无法直接应用基于直井的岩石物理解释模型.为了解决这一问题, 本文以页岩波速各向异性实验数据为基础, 引入交错网格有限差分算法, 首先在直井井孔模型中(VTI介质)模拟了声波的发射和接收, 随后通过弹性系数矩阵的Bond变换, 模拟了在井斜角不为0的情况下(TTI介质)井孔中的声场传播, 以任意井斜角与井斜角为0情况下纵波慢度差值相对值为纵坐标, 以相对应的井斜角为横坐标, 建立起了纵波各向异性校正公式.模拟与应用结果表明: 井孔中波形曲线与实轴积分法(RAI)得到的波形曲线一致, 同时利用慢度相似相关算法(STC)得到的地层纵波慢度与给定的实验测量值吻合很好, 在28种地层弹性参数的情况下, 平均相对误差为2.3%;纵波慢度差值相对值与井斜角关系曲线显示, 在井斜角小于30°的条件下, 纵波慢度差值相对值变化较小, 随着井斜角大于30°后, 相对值变化增大, 当井斜角为90°也就是水平井模式下, 纵波慢度差值相对值达到最大.根据纵波各向异性校正公式, 对水平井纵波曲线进行了慢度校正.利用纵波校正前后计算的水平井有效孔隙度与导眼井岩心分析有效孔隙度相对误差分析表明, 纵波校正后计算的有效孔隙度计算精度有了明显的提高, 证明了该方法具有非常好的应用效果, 可用于页岩水平井纵波慢度校正.

     

新疆巴楚地区金伯利质角砾橄榄岩物质组成及含矿性研究

本文讨论出露于新疆巴楚瓦吉里塔格地区的一种角砾状超镁铁岩,其结构、成分复杂,由超镁铁岩包体、斑晶（或捕晶）及基质三部分构成。超镁铁岩包体常见单辉辉石岩、纯橄岩,其次有少量橄榄辉石岩。研究表明均属基性岩浆结晶岩,本次研究未见幔源橄榄岩包体;斑晶主要为橄榄石,次为金云母,基质由微晶（10~40μm）单斜辉石、钙铁榴石、钙钛矿、磁铁矿（或含钛磁铁矿）、蛇纹石、碳酸盐及金属硫化物等组成;捕晶包括单斜辉石、褐色角闪石、磷灰石、含钛磁铁矿等。多种地球化学判别图均指示其属金伯利岩类,但低MgO和低Mg#比值、高TFe2O3和CaO等区别于世界典型金伯利岩。与典型金伯利岩有相似之处,该类岩石均具有向右陡倾的REE配分型式,但（La/Yb）n比值略偏低;微量元素蛛网图也与典型金伯利岩基本一致,仅显示更加富集不相容元素,具有更显著K, Ti负异常,且部分样品出现Rb, Zr, P负异常,指示其源区地幔交代程度偏低。鉴于岩石的产状、结构构造、矿物组合和地球化学性质近似于金伯利岩,但缺少高铬铬铁矿、镁铝榴石、镁钛铁矿等金伯利岩指示矿物,故不属典型的金伯利岩,可称之为金伯利质角砾橄榄岩。就少量研究样品所示信息,该类岩石不具有寻找金刚石的潜在远景,但鉴于巴楚及邻区尚有许多角砾状超镁铁岩岩墙和岩脉出露,该区金刚石成矿条件有待更进一步的研究。     

栅格图形数据叠置分析方法及其在城镇规划中应用

本文介绍用栅格数据进行图形叠置分析的方法,提出矢量数据转换为栅格数据的链拓扑法、减少栅格数据存储量的编码方法、根据模糊综合评判原理建立多因素综合评价数学模型以及叠置评价结果的类型格网合并的边界跟踪算法,最后给出一个实例说明栅格图形数据叠置分析方法在城镇规划中的应用。     

利用二维DLT及光束法平差进行数字摄像机标定

利用共线方程和二维DLT之间的对应关系导出了由二维DLT的8个参数表达的主纵线方程式,讨论了主点初值的求解方法,给出了单张像片摄像机参数分解不惟一性及临界序列的证明。详细推导了利用二维DLT参数分解摄像机外方位元素初值的实用算法,论述了利用光束法平差进行摄像机标定的数学模型。实际图像数据实验取得了很好的结果,验证了本文所提出的摄像机标定算法的可行性。     

超高压变质带金刚石研究进展

金刚石在自然界中非常稀缺,是极其珍贵的矿产资源。其中,变质成因的金刚石陆续在世界范围内多个超高压变质带被发现,更新了人们对超高压变质作用和板块构造运动的认识,推动了超高压变质带动力学的研究,引起了学术界的广泛关注。然而,对于变质成因金刚石的形成机制,还未形成一致的认识。本文在综述前人研究成果的基础上,介绍了金刚石形成的大地构造背景,研究了其形成所需的地质环境和物理化学条件,并着重对超高压变质带金刚石成因机制进行了探讨,全面分析了温压条件、XCO₂、fO₂对变质成因金刚石形成的影响。对比发现,当变质作用的峰期压力达到金刚石稳定域且具有较高的温度时(压力> 3 GPa,温度为600～1 000℃),有利于金刚石的形成,C-O-H流体中的较高CO₂含量和较低氧逸度也是超高压变质带金刚石形成的必需因素。另外,金刚石形成后折返早期的降温作用有利于其保存。     

大洋地幔橄榄岩-铬铁矿中的金刚石和深地幔再循环

全球多地蛇绿岩型地幔橄榄岩和铬铁矿中发现微粒金刚石,并在中国西藏南部和俄罗斯乌拉尔北部的蛇绿岩铬铁矿中发现原位产出的金刚石,认为是地球上金刚石的一种新的产出类型,不同于金伯利岩型金刚石和超高压变质型金刚石。它们与呈斯石英假象的柯石英、高压相的铬铁矿和青松矿等高压矿物以及碳硅石和单质矿物等强还原矿物伴生,指示蛇绿岩中的这些矿物组合形成于深度150~300 km或者更深的地幔。金刚石具有很轻的C同位素组成（δ13C-18‰~-28‰）,并出现多种含Mn矿物和壳源成分包裹体。研究认为它们曾是早期深俯冲的地壳物质,达到>300 km深部地幔或地幔过渡带后,经历了熔融并产生新的流体,后者在上升过程中结晶成新的超高压、强还原矿物组合,通过地幔对流或地幔柱作用被带回到浅部地幔,由此建立了一个俯冲物质深地幔再循环的新模式。蛇绿岩型地幔橄榄岩和铬铁矿中发现金刚石等深部矿物,质疑了蛇绿岩铬铁矿形成于浅部地幔的已有认识,引发了一系列新的科学问题,提出了新的研究方向。     

高密度电法在岩溶探测中的应用

高密度电法是一种适用性广泛的电阻率勘探方法,拥有多种电极排列装置选择。不同装置由于电极分布方式的差异,在针对不同地质目标和测区环境时往往表现出明显不同的探测能力。为探究各排列装置的特点及其适用情况,推导了2.5D电位的变分问题,采用Delaunay三角化算法实现了非结构化网格剖分,进而实现了有限元正演模拟。结合实践,对常见地质情况进行建模,分别使用温纳α、温纳β、施伦贝谢尔、偶极−偶极4种装置开展正反演对比研究。研究发现：(1) 在探测未知区域单个异常体时,温纳β、施伦贝谢尔2种排列的效果更好;(2) 在探测相邻较近的多个异常体时,应优先采用具有更高水平分辨率的施伦贝谢尔和偶极−偶极排列;(3) 对于低阻破碎带的勘探,温纳β和偶极−偶极排列更为适用;(4) 当地层具有较好分层界限时,4种排列均可体现良好的分层效果。因此,高密度电法勘探实践中应综合考虑不同排列装置的特点,根据具体情况选用合适的多种排列方式进行测量,并对采集数据进行综合对比分析,以获得更为准确的解释结果。

     

计算最小走时和射线路径的界面网全局方法

用慢度分块均匀正方形模型将介质参数化,仅在正方形单元的边界上设置计算结点,这些结点构成界面网．根据Ｈｕｖｓｅｎｓ和Ｆｅｒｍａｔ原理,由不断扩张、收缩的波前点扫描代替波前面搜索,在波前点附近点的局部最小走时计算中对波前点之间的走时使用双曲线近似,通过比较确定最小走时和相应的次级源位置,记录在以界面网点位置为指针的３个一维数组中．借助这些数组通过向源搜索可计算任意点（包括界面网以外的点）上的全局最小走时和射线路径．这一方法不受介质慢度差异大小限制,占内存少,计算速度较快,适于走时反演和以Ｍａｓｌｏｖ射线理论为基础的波场计算．     

菱刈矿床矿脉群形成的构造环境

本文主要介绍菱刈矿床矿脉的赋存状态、裂隙体系分析、矿脉群赋存的构造环境。菱刈矿床属脉状浅成热液型金银矿床。矿床成因与裂隙系统的形成关系密切。矿脉群大致位于基底岩石(四万十层群)隆起的中央部位。     

EXPERIMENTS OF A REDUCED GRID IN LASG/IAPWORLD OCEAN GENERAL CIRCULATION MODELS (OGCMs)

Due to the decrease in grid size associated with the convergence of meridians toward the poles in spherical coordinates, the time steps in many global climate models with finite-difference method are restricted to be unpleasantly small. To overcome the problem, a reduced grid is introduced to LASG/IAP world ocean general circulation models. The reduced grid is implemented successfully in the coarser resolutions version model L30T63 at first. Then, it is carried out in the improved version model LICOM with finer resolutions. In the experiment with model L30T63, under time step unchanged though, execution time per single model run is shortened significantly owing to the decrease of grid number and filtering execution in high latitudes. Results from additional experiments with L30T63 show that the time step of integration can be quadrupled at most in reduced grid with refinement ratio 3. In the experiment with model LICOM and with the model’s original time step unchanged, the model covered area is extended to the whole globe from its original case with the grid point of North Pole considered as an isolated island and the results of experiment are shown to be acceptable.