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

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

基于精确Zoeppritz方程的储层含油气性预测方法

流体指示因子和泊松比作为重要的指标参数在储层含油气性预测中发挥着至关重要的作用,大量学者开展了这两个参数的直接反演研究.然而,现有反演方法主要是以精确Zoeppritz方程的近似公式为正演方程,近似公式诸多的假设条件及较低的计算精度极大地限制了这类方法在复杂储层的应用效果.因此,为了提升储层含油气性预测精度,文中提出了一种新的基于精确Zoeppritz方程的流体因子和泊松比反演方法.首先,借助敏感流体因子定量分析法对现有流体指示因子进行优选,并将传统形式的精确Zoeppritz方程改写为包含该优选流体指示因子和泊松比的新形式.然后,基于新方程构建贝叶斯理论框架下的非线性反演目标函数.同时,为了进一步提升流体因子和泊松比对储层的刻画精度,在假设背景先验模型服从高斯分布的同时引入服从微分拉普拉斯分布的块约束项.最后,借助泰勒级数展开对上述非线性目标函数进行求解.合成数据和油田数据验证结果表明新方法能够稳定合理地估计流体指示因子和泊松比,且精度远高于基于近似公式的传统方法.此外,实际数据测试表明新方法能够有效提升储层含油气性预测精度,降低不确定性.

     

用切贝谢夫配点法消除地球自振常微分方程组在地心处的奇异性

用切贝谢夫配点法求解地球自振常微分方程组,无需进一步改化即可消除这组方程在地心处的奇异性,并能获得高精度的结果。     

不同方案求解非线性化学动力学方程组的比较

在一个简化的气相化学模式的基础上,比较了Hybrid、QSSA、Sklarew三种计算方案求解非线性化学动力学方程组的差异,并引入误差分配和线性组合两种质量守恒技术,分析它们对模拟结果的影响。研究结果表明：不同方案对有机烃类浓度的计算几乎没有影响,但对其它气态物浓度有一定的影响。用QSSA和Slarew方案预测的结果更为相近,而采用Hybrid方案得到的气态物浓度的平衡时间略为提前。从计算效率来看,QSSA方案优于其它两种方案,两种质量守恒技术均能改善模拟结果。     

二维全球气候模式的发展——由对流层向平流层的扩展

采用Arakawa的σ-p混合坐标,将李维亮等人发展的一个二维全球原始方程模式由对流层扩展到平流层中层(10 hPa)。垂直分层采用等lnp方案。发展后的模式不仅在计算上具有很好的稳定性,而且能较好地模拟平均温度场、平均纬向风场和平均经向风场及其季节性变化。同时,在模式对流层之上加上平流层后改善了对流层上层环流的模拟。     

地震波数值模拟的非规则网格PML吸收边界

以格子法为基础,以声波方程为例研究非规则网格PML(Perfectly Matched Layer)方法.本方法的核心是建立局部坐标系下的分裂方程和基于积分近似的微分方程弱形式.该非规则网格模拟方法允许在计算域内设置任意形状的人工边界.对于二维半空间问题,与采用矩形人工边界相比,采用半圆形人工边界可减少计算量20%以上.采用光滑的曲边界,不仅可减少计算区域,还可避免常规的PML吸收边界在吸收带角点区域的特殊处理.本方法事先计算和存储边界单元的局部几何参数,在计算的每一时间步查表调用这些参数,与常规的直边界PML方法相比,不增加任何计算量.     

Na-K-Mg三角图修正与Na-K温标选取

Na-K-Mg三角图作为水-岩平衡状态初步判定的方法,未考虑具体温标公式和矿物成分对平衡状态的影响,因此通过基于不同Na-K温标公式修正后的Na-K-Mg三角图判断水岩平衡状态,并利用适合的Na-K温标公式计算地下热储温度,对划分地热系统成因类型和开发利用地热资源具有重要意义。研究发现,修正后的Na-K-Mg三角图中任意一点,对某种温标公式确定的曲线处于平衡状态,对另一温标公式确定的曲线则有可能处于部分平衡状态;且温度越低,这种现象越显著。最贫Mg上边界和最富Mg下边界是由不同的温标公式确定的,随温度具有明显的分段性。利用各Na-K温标公式作出T-lg（Na/K）曲线图,绘制曲线交点对应的基于各\     

SPOT5影像严格几何纠正方法研究

根据SPOT卫星的成像原理和特点,讨论了对SPOT5卫星影像按共线方程法进行有控制点的几何纠正的方法,经对关中、陕北和陕南三种地形的影像纠正结果进行对比检验,证明了该方法精度达到了较高水平。     

徐州市侧柏毒蛾危害的气象条件及预测

利用现有实验结果和统计、调查资料，分析了侧柏(学名Playcladus Orientalis(L．)Franco)林主要虫害与气象因子的关系，建立起预报方程，并对林业部门开展了专业气象服务，为害虫防治工作争取了有利时机，拓展了气象服务领域，具有明显的社会经济效益。     

Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.