华罗庚数学思维训练导引三年级[1] 2_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
喜欢此文档的还喜欢
华罗庚数学思维训练导引三年级[1] 2
阅读已结束，如果下载本文需要使用
想免费下载本文？
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
你可能喜欢Execl宏在井斜数据计算中的应用50
上亿文档资料，等你来发现
Execl宏在井斜数据计算中的应用50
第30卷第6期；物探化探计算技术；2008年11月；文章编号:(；Execl宏在井斜数据计算中的应用；王高成；1,2；,林承焰,郭玲,许磊,戴鹤伟；1341；(1.中国石油大学,山东东营.中；摘要:为改进在大斜度井及井斜方位角变化大的情况下；进行测斜计算所造成的计算结果误差较大,以及地下构；常用的直线法,在
第30卷　第6期物探化探计算技术　　　　　　2008年11月文章编号:(94―06Execl宏在井斜数据计算中的应用王高成1,2,林承焰,郭　玲,许　磊,戴鹤伟1341(1.中国石油大学,山东东营　.中石油浙江油田分公司,浙江杭州　.四川石油管理局物探公司,四川成都　.中国石油勘探开发研究院,北京　100083)摘　要:为改进在大斜度井及井斜方位角变化大的情况下,采用目前地质研究中常用的直线法进行测斜计算所造成的计算结果误差较大,以及地下构造形态严重歪曲的情况,利用MicrosoftExcel内嵌的宏语言编写了Devprocess宏程序为原理,,实现了自动计算多井的井斜数据。经与Petrel,度高、数据输入输出方便考价值关键词:;直线法;井斜校正;Excel宏中图分类号:TE122.35　　文献标识码:A常用的直线法,在大斜度井及井斜方位角变化大的+0　前言在油藏描述工作中常常遇到斜井,而斜井将导致井底发生位移。此外,斜井测深都大于实际的铅直深度,尤其是在大角度造斜情况下,如不适当处理斜井,势必会在编制剖面图、构造等值线图乃至[1]等厚度图时,造成地下构造形态的严重歪曲。在地质研究中,需要通过对测斜数据计算,得出以井口位置为坐标原点的各测量点的正北、正东和垂直位置,以及水平位移、位移方位等,从而为进一步井斜校正作准备。常用的测斜数据计算方法有二[2]大类:①将相邻二井斜测点视为一直线,该算法较简单,精度差;②将相邻二井斜测点视为平面曲线,算法稍复杂,精度较高,如最小曲率法、圆柱螺线法和校正平均法等。从方法原理分析,平面曲线类的算法精度高于直线类的算法。作者在本文中采用我国钻井行业标准规定使用的校正平均法,通[3,4]过ExcalVBA编程技术,实现快速计算多井测斜数据。通过对比,该方法能弥补目前地质研究中收稿日期:　　　　改回日期:情况下,所造成的计算结果误差较大以及地下构造形态严重歪曲的情况,同时该程序使用起来也更加方便。1　原理在实际工作中,为了方便计算测斜数据,我们常采用直线段的方法。直线法把斜井段视为直线,斜井段在垂直及水平方向上投影也为直线段(见下页图1及图2)。主要的计算公式为:斜井垂深:D=D0+n[5]6ΔDi=1minmicos&i东西向偏移:dEW=南北向偏移:dSN=6ΔDi=1nαsin&i?sinαsin&i?cos6ΔDi=1mi斜井水平偏移距:d=ΔDmisin&i?sinαi)+(ΔDmisin&i?cosαi)i=1i=16n26n26期王高成等:Execl宏在井斜数据计算中的应用　　495位置,将大大提高计算精度。校正平均法假设井斜段形状为一条圆柱螺线,其二端与上、下二测点处井眼方向相切,在水平投影图上是圆弧,在垂直剖面图上也是圆弧(如图3所示)。郑基英教授在这个假设的基础上推导出“圆柱螺线法”,经数学处理便提出了“校正平均法”,测斜段计算公式如[6～7]下:ΔD=fD?ΔDm?cosαcΔLP=fD?ΔDm?sinαcΔN=fH?ΔDm?sinαc?cos&cΔE=fH?ΔDm?sinαc?sin&c　图1　斜井直线法水平投影Fig.1　Horizontalprojectionofdeviationwellofstraight2linemethod其中2αfD=1-Δ/24Δα+Δ&)/24fH=1-(αc=(αi-1+αi)/2c=(i-1+&&i&i-1,只是在平均角法的公式基础上乘上校正系数fH和fD,因而称之为校正平均角法。对各测斜段计算得ΔN、到的正北、正东和垂直投影方向的长度ΔD、ΔE进行相加,就可以得到以井口位置为坐标原点的各测量点的正北、正东和垂直位置了。按照这个原理,利用Excel宏语句编写了可以实现多井测斜数据计算的宏程序―Devprocess,精度较直线法有了极大提高。在井深处于2000m左右,井斜角大于35°的情况下,利用该宏程序计算出来的结果,与Schlumberger公司的Petrel等商业化软件计算结果之间的误差能控制在0.5m之内。而该宏程序以MicrosoftExcel电子表格为基础,因而其使用界面简洁友好,数据输入和结果输出十分方便,运行方法简单,兼容性和可移植性强,使用起来比22图2　斜井直线法垂向投影Fig.2　Verticalprojectionofdeviationwellofstraight2linemethod然而,实际的斜井井眼轨迹是展布在三维空间中的一条曲线,是连续弯曲的,这在井斜角较大时,利用直线段方法计算的结果,与实际情况之间偏离较远,差别较大。据统计,在井深处于2000m左右,斜井的井斜角达到35°时,利用直线法计算的垂深和实际垂深之间的差别,能达到10m以上。这将严重扭曲地下实际地质形态,给后续研究工作带来错误的认识。用校正平均法计算井筒的空间图3　校正平均角法示意图[6]Fig.3　Sketchmapofcorrectedaverageprinciple496物探化探计算技术30卷Petrel等商业化软件方便得多,更适合地质研究人打开宏的列表,或者按住“ALT+F8”,进入VBE环境,找到宏名称-“Devprocess”,点击“运行”。2.3　数据输出界面Devprocess宏程序根据输入的数据列,能自动员用于井斜校正工作。2　Devprocess的使用方法及步骤2.1　数据输入界面计算出垂深列、横偏移距列和纵偏移距列(如图4所示),将其复制到一个新建的文件中保存,以免影响该程序的后续使用。同时,对计算结果稍加整理,就可实现测点数据的自动计算功能。需要输入的数据包括井名列、测点序号、斜深、井斜角和方位角(图4)。2.2　运行宏Excel宏语言是由微软开发出来,能在其桌面3　Devprocess程序的源代码校正平均法原理简单,程序实现起来也较容易,但是在Devprocess程序的编写中还需考虑以下三个问题:(1)应用程序中执行通用自动化(OLE)任务的编程语言[8～9],编写简单,省去了对文件的输入、输出操作,而且结果显示在Excel表中,易于控制。在使用时需将安全设置为低,这样宏将被自动启用。从MicrosoftExcel菜单中顺序选择:工具|宏|宏就可图4　Devprocess宏数据输入(A、B、C、E、F列)及运行结果(D、G、H列)界面Fig.4　Interfaceofdataloading(A,B,C,E,F)andcalculatedresult(D,G,H)ofDevprocess6期王高成等:Execl宏在井斜数据计算中的应用　　497关数据,从而实现多井自动计算的功能。(2)二相邻测点间的方位角的变化应限于±180°之间,任何大于这个范围的变化意味着在井轨迹反方向的变化,在这种情况下计算&c和Δ&i时,需采用以下方法:当　&i-&i-1&180°时Δ&i=&i-&i-1-360°　　　　&c=(&i+&i-1)/2-180°当　&i-&i-1&-180°时Δ&i=&i-&i-1+360°　　　　&c=(&i+&i-1)/2+180°(3)针对有些井有第0测点,而有些没有第0测点的情况,需要对没有第0测点的的井添加第0测点。第0测点是第1测段的上部测点,该测点不是实测的,而是人为规定的。在编程时,应先判断出没有第0测点的井,然后若第1测点的井深大于25m,则第0测点的井深应比第1测点的井深小于25m,且井斜角规定为零;1[6]Selection.InsertShift:=xlDownCells(i+1,1)=Cells(i+2,1).ValueCells(i+1,7)=0Cells(i+1,8)=0IfCells(i+2,3).Value&=25Then第’1测点井深小于或等于25mCells(i+1,3).Value=0Cells(i+1,5).Value=0EndIfIfCells(i+2,3).Value&25Then′第1测点井深大于25mCells(i+1,3).Value=Cells(i+2,3).Value-25Cells(i+1,5).Value=0EndIfi+1,+2,6).Value+1,4i3).ValueIf等于25m时,零。’3333333SubDevprocess()Dimii,i,j,k,m,wellrow,maxrowAsIntegerDimi1AsIntegerDimx,y,azim,incl,md,tvd,f1,f2,azimz,in2clz,inclc,azimc,piAsDoublepi=3.14159Cells(1,4).Value=&垂深&Cells(1,7).Value=&横偏移距&Cells(1,8).Value=&纵偏移距&m=Range(&A65536&).End(xlUp).RowFori=1Tom-1　IfCells(i,1)&&Cells(i+1,1)Then′根据IfCells(i+1,Cells(i+1,Cells(i+1,Cells(i+1,Cells(i+1,EndIfIfCells(i+1,5).Value&&0Then′没有第0测点,填加第0测点Range(&A&&i+1&&:&&&O&&i+1).Select′对第0测点预处理完毕　　k=1x=0:y=0:azim=0:incl=0:md=0:f1=0:f2=0:azimz=0:inclz=0:inclc=0:azimc=0α,inclz-Δ&f1-fD,f2-fH,′md-ΔDm,azimz-Δinclc-&c,azimc-αctvd=Cells(i+1,4).ValueDoWhile(Cells(i+k,1)=Cells(i+k+1,1)Andi+k+1&=mAndCells(i+k+1,1).Value&&&&)′判断出同属一口井的测斜数据,并利用校正平均法对该口井进行计算i1=i+k+1md=Val(Cells(i1,3).Value)-Val(Cells(i1-1,3).Value)IfVal(Cells(i1,6).Value)-Val(Cells(i1-1,6).Value)&180Then′对同一测段内井斜方位井号,判断出同一口井5).Value=0Then′对第0测点预处理,有第0测点的情况6)=Cells(i+2,6).Value4)=Cells(i+1,3).Value7)=08)=0角变化的绝对者超过180°的做预处理　azimz=Val(Cells(i1,6).Value)-Val(Cells(i1-1,6).Value)-360　azimz=azimz/1803pi　azimc=(Val(Cells(i1,6).Value)+Val(Cells(i1-1,6).Value))/2-180　azimc=azimc/1803pi　　ElseIfVal(Cells(i1,6).Value)-Val(Cells(i1-1,6).Value)&-180Then　azimz=Val(Cells(i1,6).Value)-Val(Cells(i1-1,6).Value)+360498物探化探计算技术f2=1-(azimz3azimz+inclz3inclz)/24tvd=f13md3Cos(inclc)+tvdx=f23md3Sin(inclc)3Sin(azimc)+xy=f23md3Sin(inclc)3Cos(azimc)+yCells(i1,7).Value=xCells(i1,8).Value=yCells(i1,4).Value=tvdk=k+130卷　azimz=azimz/1803pi　azimc=(Val(Cells(i1,6).Value)+Val(Cells(i1-1,6).Value))/2+180　azimc=azimc/1803pi　　Else　azimz=Val(Cells(i1,6).Value)-Val(Cells(i1-1,6).Value)　azimz=azimz/1803pi　azimc=(Val(Cells(i1,6).Value)+Val(Cells(i1-1,6).Value))/2　azimc=azimc/1803pi　　EndIfinclz=Val(Cells(i1,5).Value)-Val(Cells(i1-1,5).Value)inclz=inclz/1803piinclc=(Val(Cells(i1,5).Value)+Val(Cells(i1-1,5).Value))/2inclc=inclc/1803pif1=1-)　Loop　i=i+k-1　EndIfNextiEndSub’3333程序分析结束3334(最大的井)进行计算(见表1),并将计算结果46°表1　CHEN3～60井测斜数据及计算结果分析表Tab.　CalculatedresultofsurveyscomputationofwellCHEN3260序号012345测深(m)0.36.方位角Devprocess宏程序计算结果(°)0.000.220.450.730.580.38(°)229...90291.10直线法计算结果Petrel计算结果垂深/m△E(m)△N(m)垂深(m)△E(m)△N(m)垂深(m)0...36394.190.00-0.40-0.52-0.76-1.08-1.320.00-0.35-0.47-0.64-0.69-0.640...36394.200.000.00-0.08-0.21-0.56-0.850.000.00-0.07-0.25-0.36-0.360...36394.19△E(m)0.00-0.40-0.51-0.76-1.08-1.32△N(m)0.00-0.35-0.48-0.62-0.68-0.64(略去6～33号数据)242...85.685.673.553.585.0599.2098.80.0818.24.-4.-4.-5.-5.-5.-3.17.23.-4.01-4.14-4.54-5.11-5.60-4.55.297.355.0318.23.-4.51-4.64-5.10-5.63-5.36-3.62101.05.4.7019.10.01(略去40～68号数据)7555.-13..60-13..945..24-15..67-14..846..32-17..77-16..846..41-18..86-18..846..51-20..94-19..145..92-22..34-21..045..72-23..14-22.45567.47.116.149.-14.54-16.10-17.77-19.52-21.42-22.75-23.95468....60包含各类专业文献、行业资料、各类资格考试、生活休闲娱乐、高等教育、应用写作文书、专业论文、幼儿教育、小学教育、Execl宏在井斜数据计算中的应用50等内容。　
　EXECL 测量超高计算 2页 2下载券 Execl宏在井斜数据计算中... 7页 免费喜欢...EXECL 计算使用方法 1、查找重复内容公式:=IF(COUNTIF(A:A,A2)&1,&重复&,...&5,104.6191.42
&17,851.240173.67
&26,871.940184.300
新赛股份代码：600540
16.400.20&1.23%
昨收盘:16.20
今开盘:16.10
最高价:17.05
最低价:16.01
均价:16.49
成交:143909手
流通:3.03亿
市值:59.41亿元
换手:4.75%
委比:-89.0
委差:-3662