曹书锦
- 作品数：33 被引量：57H指数：5
- 供职机构：湖南科技大学更多>>
- 发文基金：国家自然科学基金湖南省自然科学基金湖南省教育厅科研基金更多>>
- 相关领域：天文地球自动化与计算机技术文化科学建筑科学更多>>
- 基于粒子群算法的重力张量反演研究被引量：3
- 2013年
- 重力梯度张量的定义是对重力位求二阶导数,相比于传统的布格重力异常,它在反映地下密度异常分布上有着更高的灵敏度,且能够进一步准确直接地反映目标体的边界。但对单个张量分量作反演时,可能会丢失一部份有用的信息,从而造成反演结果的误差。而全张量反演是将重力的五个梯度张量联合起来进行反演,这样做可以综合更丰富的场源信息。相比于传统做法中的布格重力异常反演和单重力张量分量反演,得到的反演结果不仅有了更高的分辨率,在识别目标体特征上也有更好的效果。粒子群算法是一种基于群体智能的优化迭代算法,这里利用粒子群算法对重力张量单分量、布格重力异常和全张量分别进行反演,并对结果进行简要分析。
- 孙中科朱自强鲁光银曹书锦
- 关键词：粒子群优化算法反演
- 基于傅立叶变换的B样条函数欧拉解概率密度成像方法
- 本发明提供一种基于傅立叶变换的B样条函数欧拉解概率密度成像方法，确定待测区域范围，测量重力矢量数据及重力梯度张量数据；构建三维重力梯度张量欧拉反褶积方程，确定滑动窗口大小，构建线性方程组，并获取欧拉解，输出欧拉解集；将欧...
- 曹书锦毛雅静鲁光银杨博郭晓旺陈新跃朱自强陶川溢邓意怀
- 文献传递
- 基于单元细分H-自适应有限元全张量重力梯度正演被引量：7
- 2010年
- 在重力梯度各阶张量正演中,难于求取复杂地质体的解析式,转而利用已知解析解的简单规则几何形体近似剖分地质体,其会引起较大的几何拟合误差.本文首先从全张量重力梯度解析公式出发结合等参变换基于有限元分析给出了满足地下复杂地质体正演的型函数,引入体积分的误差指示器用以估计全张量重力梯度误差,提出了适用于全张量重力梯度正演的自适应计算的策略及迭代算法,数值算例证明了本文方法的正确性和有效性.
- 曹书锦朱自强鲁光银曾思红郭文斌
- 关键词：非结构化网格自适应有限元
- 一种基于多层级方法重磁快速反演方法
- 本发明公开了一种基于多层级方法重磁快速反演方法，首先确定观测高度和地球物理模型，离散地球物理模型，提出虚拟线线观测网格/虚拟面面观测网格，获得含多个BTTB矩阵重磁正演核矩阵；构造基于多尺度Haar小波的插值算子和限制算...
- 曹书锦郭晓旺杨博毛雅静陈新跃朱自强鲁光银邓意怀马致远
- 重力场欧拉反褶积最优解提取被引量：1
- 2017年
- 在传统重力场欧拉反褶积中,单一构造指数难于表征多个异弱常源;而枚举构造指数易于导致欧拉解过度发散、欧拉解解集庞大且其中谬解极多。为此,在采用非预测欧拉反褶积方法,避免枚举构造指数的基础上,采用截断奇异值分解法解欧拉齐次方程,以评价欧拉解的不稳定性,在此基础上,构建估计函数以保留位场总场效应相对较弱的异常源。最后通过数值例子进一步表明了方法的有效性和可靠性。
- 周勇曹书锦侯萍萍杨子龙
- 关键词：欧拉反褶积最优解奇异值分解
- 一种基于多变量核密度欧拉解概率密度成像方法
- 本发明提供一种基于多变量核密度欧拉解概率密度成像方法，首先通过测量或位场变换获得重力矢量及重力梯度张量数据，构建三维重力/重力梯度张量欧拉反褶积方程，逐一移动滑动窗口，获得一系列含异常源空间位置[x，y，z，N]欧拉解集...
- 朱自强杨博鲁光银曹书锦毛雅静张亮张东兴郭晓旺陈新跃
- 文献传递
- 基于OpenMP的重力张量并行正演被引量：1
- 2013年
- 重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏。对于复杂的重力密度模型来说,张量的解析公式是很难推导的。为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量。采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度。有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题。这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路。
- 余涛朱自强鲁光银曹书锦
- 关键词：正演OPENMP
- 一种重磁快速正演方法
- 本发明公开了一种重磁快速正演方法，首先确定观测高度和地球物理模型，离散地球物理模型，提出虚拟线线观测网格/虚拟面面观测网格，获得含多个BTTB矩阵核矩阵，使用多尺度Haar小波构建一系列对应于不同层的核矩阵，从而实现基于...
- 曹书锦毛雅静郭晓旺杨博朱自强鲁光银马致远陈新跃邓意怀
- 重力场欧拉反褶积最优化滑动窗口选择
- 2021年
- 针对欧拉反褶积法中使用单一构造指数难于表征多个异常源,而使用不同尺寸的滑动窗口或通过枚举构造指数易于导致欧拉解集庞大,增加不必要的后续处理解释工作。为此,利用聚类分析划分欧拉反褶积解集,采用截断奇异值分解法解欧拉齐次方程,以评价因欧拉反褶积所构成的大型线性方程组条件数过大而造成欧拉解的非唯一性和不稳定性问题。在此基础上,进一步对比分析相同尺寸异常源在不同尺寸滑动窗口下所产生的欧拉解集中优解的比例,借此分析滑动窗口尺寸和异常源尺寸间的比例关系,并以此为依据选择最优化窗口。模型试验结果表明:当测量点距为异常源尺寸的1/10时,最优化窗口应为异常源尺寸的3倍且小于勘探区尺寸的1/2。
- 杨博曹书锦毛雅静周勇周勇朱自强
- 关键词：欧拉反褶积
- 重力梯度张量种植反演被引量：1
- 2017年
- 针对在重力梯度张量正演中计算耗时过长和核矩阵内存消耗过大等制约反演实施的瓶颈问题,在L1范数的基础上,引入种植反演,用累加求和分析替换迭代求解,避免计算或存储反演核矩阵,以减少内存占用和加快反演迭代;针对种植反演容易导致相邻异常源相互侵入的问题,引入一个基于位场水平衰减特性加权函数来限制密度吸引子的作用范围,以期使密度吸引子忽略较远的异常源,抑制相邻异常源相互干扰。反演结果及分析表明重力及重力梯度张量种植反演所需计算机内存小和水平衰减特性加权函数能有效的抑制相邻异常源的侵入。
- 曹书锦朱自强鲁光银
- 关键词：重力梯度张量L1范数
