辛普森公式的推广形式及应用.pdfVIP

第 28卷 第 1期 山东 理 工 大 学 学 报 (自然 科 学 版) VoI．28No．1 2014年 1月 JournalofShandongUniversityofTechnology(NaturalScienceEdition) Jan．2014 文章编号 ：1672—6197(2014)01—0029—03 辛普森公式的推广形式及应用 刘智颖 ，王向公 ，任威龙。，龙耀萍。 (1．长江大学 油气资源与勘探技术教育部重点实验室，湖北 武汉 430100； 2．长江大学 地球物理与石油资源学院，湖北 武汉430100； 3．北京吉奥特能源项 目责任有限公司，北京 100195) 摘 要 ：鉴于当前常用数值积分方法的不足，以及直角坐标系下辛普森公式的局限性，研究并提出 了辛普森公式在三维柱坐标系和球坐标系下的推广形式，并将其应用于地球物理测井 中自然伽马 射线强度 函数数值积分运算． 关键词 ：数值积分；辛普森公式；三维柱坐标；三维球坐标 中图分类号：O241．4 文献标志码：A ThededucingformsofSimpson formulaanditsapplication LIU Zhi—ying ，W ANG Xiang—gong ，REN W ei—long。，LONG Yao—ping。 (1．KeyLaboratoryofExplorationTechnologiesforOilandGasResources，MinistryofEducation， YangtzeUniversity，Wuhan430100，China； 2．CollegeofGeophysicsandOilResources，YangtzeUniversity，Wuhan430100，China； 3．BeijingGeoTechEnergyProjectLimitedCompany，Beijing100195，China) Abstract：ConsideringthelimithnessofnumericalintegratorandtheapplicationoftheSimpson formulainrectangularcoordinatesystem ，weproposedthededucingformsoftheSimpsonformu- la，whichcouldbeusedin3D cylindricalcoordinatesand3D sphericalcooridinates，andapplied them totherelevantalgorithm intheintegralofGR intensityfunction． Keywords：numericalintegrator；Simpsonformula；3D cylindricalcoordinate；3D sphericalcoor— dinates 当前，简单而常见的数值积分方法是无限细分法(亦称 “黎曼和”)．在应用过程 中，需将区间划分成尽可 能多的份数 以保证计算精度．然而 ，较细 的划分会导致计算时间增加 ，而且 当被积函数变化幅度较大时，其误 差增大，故辛普森公式应运而生．常见的辛普森公式大都只适用于直角坐标系．事实上很多积分需要在柱坐 标和球坐标下进行，如果做坐标变换后再求解则十分繁琐而且会造成误差累积．因此，本文分析、研究并推导 了柱坐标系和球坐标系下的辛普森公式 ，且成功应用于地球物理测井 中 自然伽马射线强度 函数的数值积分 运算． 1 三维柱坐标系和球坐标系下辛普森公式的推广 1．1 辛普森公式的常见形式及其适用范围 常见的辛普森公式 (公式 (1))只能应用于直角坐标系下一元函数的数值积分 收稿 日期 ：2013一O7—26 作者简介 ：刘智颖，男，lzykkww100@163．corn；通信作者：王向公，男，wxgll@yangtzeu．edu．an 30