海洋CSEM法探测二维海底天然气水合物响应特征

针对海底天然气水合物开展了可控源电磁探测法(Controlled-source Electromagnetic Method,CSEM)响应研究,根据实际地质资料设计二维地电模型并进行正演计算,讨论了频率变化对CSEM响应结果的影响以及海底地形起伏时的CSEM响应特征。综合天然气水合物二维正演模拟结果,论证了海洋CSEM探测海底天然气水合物的可行性,研究成果可为今后的海底天然气水合物勘探提供技术支持。     

一种层状大地瞬变电磁响应正演计算的改进方法

频谱法是瞬变电磁法正演计算中常用的计算方法之一,其关键在于内层含贝塞尔函数的积分计算。对于重叠回线装置的瞬变电磁法正演的内层积分,以往采用在积分区间内寻找贝塞尔函数的零点分布,依次在相邻零点之间采用一般的数值积分法在求得各自的积分值后再叠加的方法进行计算。这种方法精度较高,但效率低。利用贝塞尔函数的大宗量渐近特性,提出一种新的计算方法。计算结果表明,该方法计算效率高、方法简单、精度可靠。     

均匀半空间回线源非中心点频率域电磁场计算

在回线源频率域非中心点磁场计算过程中,会遇到双贝塞尔函数积分问题。这里给出了一种基于静磁场毕奥-萨伐定律计算双贝塞尔函数积分的数值方法,将频率域电磁场写成由回线源直接激发的一次场和由地下导电介质感应出的二次场的和的形式,根据一次场和二次场被积函数的形态区别,对两部分的积分采用不同的计算方法。为了验证该方法的准确性,将回线源剖分成若干个电偶极子,用电偶极子叠加的电磁场与本文提出方法计算的结果对比。对比结果表明,二者的相对误差很小,这有助于研究回线源电磁场的分布规律。     

水力输送法开采海底浅层天然气水合物技术研究

海洋天然气水合物占全球水合物总储量的99%,海洋天然气水合物大部分沉积在海底浅层没有良好的覆盖层,无法采用开采传统油气资源的方法进行开采.为了开采海底浅层天然气水合物,分析了海洋天然气水合物的存在类型及地质特征,探讨了目前开采海洋天然气水合物的理论方法及局限性.根据海底天然气水合物成矿的地质特点及水合物的物理特性,提出了以水力输送技术为主要理论的海洋浅层天然气水合物开采方法.该方法克服了其他方法开采海底天然气水合物面临的技术困难,为海洋天然气水合物的开采提供了技术基础及理论支撑.     

结合有限项多指数级数近似的点源格林函数闭式算法

针对水平层状介质点电流源格林函数,提出了一种快速计算技术。该方法利用点电流源层状介质格林函数谱为单调实函数的特点,采用多指数级数近似格林函数谱,进而借助Somerfeld-Lipschtz积分实现格林函数的闭式计算。计算过程中,格林函数谱以等对数间隔抽样,近似多指数级数的指数常数保持固定并被等对数间隔抽样,近似过程只是指数级数系数谱的计算过程。计算结果表明,采用如下参数能实现点电流源格林函数高精度闭式计算:(1)λmin=10-6,λmax=10,采样点N=128,λ对数抽样;(2)T2min=10-7,T2max=103,T2对数抽样,抽样个数M=64;(3)最优正则化参数α2=10-5,采用最简单模型对解进行约束。     

冲绳海槽天然气水合物稳定带厚度的计算

水合物稳定带（HSZ）的研究对天然气水合物的成矿与分布规律以及资源评价研究都具有重要的指导意义，文章根据天然气水合物的相平衡条件和相应的压力-温度方程，参照有关的身体算法，讨论了海底天然气水合物稳定带厚度的理论计算方法，并根据稳定带厚度的理论计算方法，在Windows中文平台上编制了简单实用的程序，利用该程序对冲绳海槽水合物稳定带的厚度进行了计算，并讨论了冲绳海槽水合物稳定带厚度的分布特征。     

双重贝寒尔函数积分的数值计算

双重贝塞尔函数积分则由于贝塞尔函数互乘项的强振荡与慢衰减特性而以以应用通常的数值积分算法。本文将被积区间[0,∞]划分为[0，λ0]、[λ0，∞]两部分，应用贝塞尔函数的克尔函数表述式及后者的大宗量渐近特性，区间[λ0，∞]的双重贝寒尔函数积分可被转化为Fourier正（余）弦变换，并可利用各种快速算法对其进行数值计算；区间[0，λ0]上的双重贝塞尔函数积分的计算可直接应用一般的数值积分算法并能获得较高的计算精度；当需大量计算有共同参量的双重贝塞 尔函数积分时，其计算效率仍显不足。此时，可应用贝塞尔函数的导数关系式对[0,λ0]内的双重贝塞尔函数积分进行恒等变换，再用差商近似导数，将其转化为对贝塞尔函数本身的积分，而该积分又仅需计算一次。故本算法对双重贝塞尔函数积分的计算效率有明显提高。     

海底天然气渗漏系统水合物成藏动力学及其资源评价方法

根据海底天然气渗漏系统沉淀水合物的动力学模型,计算了海底渗漏系统中渗漏天然气沉淀为水合物的比例。如果BushHill渗漏系统海底水合物或喷溢天然气的化学组成是渗漏系统气源天然气通过沉淀水合物转变而成,则需要约3. 3% ~21. 7% (平均12. 9% )的渗漏天然气在海底沉淀为水合物。结合渗漏系统的活动时间(1万年)和天然气流量(800t/a),沉淀在墨西哥湾GC185区BushHill渗漏系统中的水合物天然气资源为0. 37×109 ~2. 43×109 m3,平均为1. 45×109 m3,与体积和含量评价方法获得的结果基本一致。     

海底天然气水合物赋藏区域地热基本模型探讨

在分析海底天然气水合物赋藏区域地热场特征的基础上,借鉴盆地地热模型理论和方法,建立了海底天然气水合物赋藏区域的地热基本模型,并详细介绍了边界条件和导热系数的确定方法.由于天然气水合物的形成、分布与实际地质和地热条件密切相关,因此加强天然气水合物赋藏海域的地质与地热基础调查是构建准确可靠的地热模型的关键.     

任意形状工作面沉陷预测的概率积分法及其应用

概率积分法是煤矿开采沉陷预测的重要方法,可适用于缓倾斜、倾斜煤层开采地表移动分布曲线符合正态分布规律的地表移动和变形计算。针对原方法对任意形状工作面存在的不足,在原有走向和倾向概率积分函数的基础上,以走向和倾向为基准划分积分函数区域,在该区域内剖分非结构化的三角形单元作为二重积分的基本单位;通过坐标旋转变换的二重积分换元法,转化为以直线为上、下限的二重积分,计算旋转后的影响半径,对新的二重积分采用复合辛普森数值积分公式,积分得出地表任意点在基本三角形开采单元下的积分下沉值;最后,通过叠加计算完成任意点、任意形状工作面概率积分法的沉陷预测计算。算法在GIS系统下实现,实例应用效果良好,可为“三下”采煤提供地表移动预测计算的支持。     

基于不同半径分段组合的微动探测方法研究

空间自相关法是一种传统的微动勘查方法。笔者在双重圆形观测台阵中不同半径组合方式下,利用空间自相关法分别求出自相关系数以及频散曲线,对比不同半径的频散曲线,结合前人研究成果,总结出大半径组合方式在低频部分的可靠性更高,小半径组合方式在高频部分的可靠性更高,进而提出了分频段拟合贝塞尔函数、分半径平均频散曲线的方法,一定程度上扩宽了可用频率范围,提高了探测精度。通过场地实验证实了方法的有效性。     

Optimization of the rate of return of the mine production process

Summary The basic foundations of an integrational method of modelling and estimating the economic effectiveness of a deep coal mine production process is presented. The method is founded on the concept of integration of smaller structural elements into larger ones in order to obtain a representational model of the whole mine production process. It involves analysis of the spatial, engineering and time structures of the mining process. The basic structural unit of the production process is the elementary time interval,T, together with a linear function of the approximate intensity of financial outlays or effects. The repeated integration of these elementary time intervals within the spatial, engineering and time structures of the mine production process can be extended to model the whole mine or discrete spatial structure elements.The model may be used to investigate the trends of financial outlays and effects, and their current values updated for arbitrarily assumed updating moments. The model can also be used as an investigation tool for analysing the economical effectiveness of the mine production process.     

地球空间数据集成多尺度问题基础研究

多尺度数据集成是地球空间数据集成中最难处理的问题。将多尺度数据的集成分解为空间和时间多尺度数据集成,在分析应用项目对数据尺度需求的基础上,就两种多尺度数据集成的传统和数据意义上的集成方法进行了详细的探讨。     

“嫦娥一号”第一幅月面遥感影像撞击坑特征

撞击坑是月表的主要地貌形态,按普通碗状型、中心隆起型、平底型及同心圆型4种类型描述"第一幅月面图像"上月表撞击坑的形态特征.应用遥感专业图像处理软件,结合白赛尔大地主题反算法,实现撞击坑分布及直径信息提取.研究发现,研究区域的撞击坑直径主要集中在300～700 m.而且撞击坑空间分布不均匀,直径越小,分布越密集.     

An MHD model for a heliospheric current sheet

A numerical solution of the full set of MHD equations shows the generation of a heliospheric current sheet during the thermal expansion of the corona. Calculations were performed for a compressible plasma taking into account dissipative terms and anisotropy of the thermal conductivity of the magnetized plasma. It is shown that the current sheet is not magnetically neutral. The sheet contains a normal component of the magnetic field, which plays a fundamental role during the formation of the sheet and in the stationary state. The sheet is stable against MHD perturbations, which are apparently carried away by the plasma flow. For the numerical scheme chosen, the minimum sheet thickness is determined by the length of the spatial integration step. The PERESVET code was used for the calculations.     

四川南江竹坝磁铁矿区磁异常数据处理应用效果

文章以四川南江竹坝矿区李子垭磁铁矿床地面高精度磁测的数据处理成果为例,提出利用化极和求取垂向二阶导数的方法,在展示矿体空间赋存状态、缩小探矿范围等方面有一定效果。     

地面-井下双源瞬变电磁立体探测模拟研究

在探测地下、井下或巷道下方的局部异常体时,地面瞬变电磁信噪较低而井下瞬变电磁发射磁距受限,因此探测深度较小。为解决这个问题,本文提出了一种地面大回线仅发射+井下重叠回线发射接收的地面-井下双源瞬变电磁立体探测方法。通过COMSOL软件数值模拟与室内水槽试验,模拟了四种不同的瞬变电磁装置(井下重叠回线,井上大回线发射+井下小回线接收,地面中心回线,地面-井下立体探测)探测。数值模拟结果发现地面-井下双源装置的感应电位曲线幅值最大且灵敏度最高;室内模型试验结果绘制成感应电动势多测道剖面图,分析比较其曲线特征,表明地面-井下瞬变电磁立体探测明显优于其他三种瞬变电磁装置,其效果非常明显、稳定且磁感应强度更高,与数值模拟结果一致。因此,地面-井下双源瞬变电磁立体探测方法可行且值得推广。     

导线环垂直激发场高精度数值模拟

导线环垂直激发场是导线环在地下产生的激发场中垂直于地磁场的分量,在地面核磁共振找水方法中,垂直激发场在地下的分布直接决定地面核磁共振响应,探讨影响垂直激发场地下分布的诸因素进而研究这些因素对地面核磁共振响应的作用,对核磁共振找水正确野外施工具有一定的指导意义,有助于采集更高质量的野外数据.利用经典Chave算法对含有贝塞尔函数的强振荡积分核进行了数值积分试验,将其计算结果和对应的解析解进行了对比,试验结果显示该算法用于强震荡函数无穷积分的计算具有良好的收敛性、较高的精度和较快的计算速度,采用Chave算法对含有贝塞尔函数的强振荡函数进行积分是可行的.在此基础上,进一步计算并数值模拟了导线环垂直激发场及其相位,分别从地磁倾角、磁偏角、地下介质导电性、深度等方面探讨了影响垂直激发场及其相位分布的主要因素.