高阻球體上不同拉線(xiàn)方向的電測(cè)深曲線(xiàn)與視電阻率斷面圖(FORTRAN代碼)_第1頁(yè)
高阻球體上不同拉線(xiàn)方向的電測(cè)深曲線(xiàn)與視電阻率斷面圖(FORTRAN代碼)_第2頁(yè)
高阻球體上不同拉線(xiàn)方向的電測(cè)深曲線(xiàn)與視電阻率斷面圖(FORTRAN代碼)_第3頁(yè)
高阻球體上不同拉線(xiàn)方向的電測(cè)深曲線(xiàn)與視電阻率斷面圖(FORTRAN代碼)_第4頁(yè)
高阻球體上不同拉線(xiàn)方向的電測(cè)深曲線(xiàn)與視電阻率斷面圖(FORTRAN代碼)_第5頁(yè)
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說(shuō)明：本文檔由用戶(hù)提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、!低阻球體上沿X軸拉線(xiàn)的電測(cè)深曲線(xiàn)PROGRAMPS_CSIMPLICITNONEEXTERNALPS4REALX0,Y0,P1,P2,H0,r0,MN,AB,PS4!測(cè)點(diǎn)xy坐標(biāo)，圍巖電阻率，球體電阻率，PARAMETER(P1=1,P2=10,H0=2,r0=1)!以下為測(cè)點(diǎn)在x軸上不同位置時(shí)的電測(cè)深數(shù)據(jù)；*X0=0Y0=0OPEN(1,FILE=P-AB-x001.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/

2、P1,AB/2ENDDOCLOSE(1)*X0=0.5Y0=0OPEN(1,FILE=P-AB-x002.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)*X0=1Y0=0OPEN(1,FILE=P-AB-x003.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/

3、2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)*X0=1.5Y0=0OPEN(1,FILE=P-AB-x004.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)*X0=2Y0=0OPEN(1,FILE=P-AB-x005.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS

4、4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)*X0=3Y0=0OPEN(1,FILE=P-AB-x006.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)!以下為測(cè)點(diǎn)在X軸上時(shí)的電阻率斷面圖數(shù)據(jù)*OPEN(1,FILE=X-AB-001.DAT,ACTION=WRI

5、TE)Y0=0DOX0=-3,3,.05DOAB=0.1,15,0.5MN=AB/100WRITE(1,*)X0,AB/6,&PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.,P1,P2,r0,H0)/P1ENDDOENDDOCLOSE(1)!以下為測(cè)點(diǎn)在Y軸上不同位置時(shí)的電測(cè)深數(shù)據(jù)I*X0=0Y0=0OPEN(1,FILE=P-AB-y001.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,

6、P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)*X0=0Y0=0.5OPEN(1,FILE=P-AB-y002.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)*X0=0Y0=1OPEN(1,FILE=P-AB-y003.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS4(X0-AB/2,0.,X0-MN/2

7、,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)*X0=0Y0=1.5OPEN(1,FILE=P-AB-y004.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)*X0=0Y0=2OPEN(1,FILE=P-AB-y005.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WR

8、ITE(1,*)PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)*X0=0Y0=3OPEN(1,FILE=P-AB-y006.DAT,ACTION=WRITE)DOAB=0.1,12,0.5MN=AB/100WRITE(1,*)PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.&,P1,P2,r0,H0)/P1,AB/2ENDDOCLOSE(1)!以下為測(cè)點(diǎn)在Y軸上時(shí)的電阻率斷面圖數(shù)據(jù)*OPEN(1,FILE=Y-AB-001.DAT,

9、ACTION=WRITE)X0=0DOY0=-3,3,.5DOAB=0.1,15,0.05MN=AB/100WRITE(1,*)Y0,AB/6,&PS4(X0-AB/2,0.,X0-MN/2,0.,X0+MN/2,0.,X0+AB/2,0.,P1,P2,r0,H0)/P1ENDDOENDDOCLOSE(1)ENDPROGRAMFUNCTIONPS2(XA,YA,XM,YM,P1,P2,r0,H0)implicitnoneREALPS2,XA,YA,XM,YM,P1,P2,r0,H0,RR,K,d,r,T,COS0RR=SQRT(XM-XA)*2+(YM-YA)*2)K=(P2-P1)/(P1+

10、2*P2)d=SQRT(XA*2+YA*2+H0*2)r=SQRT(XM*2+YM*2+H0*2)T=r0*3/(d*2*r*2)COS0=(d*2+r*2-RR*2)/(2*d*r)PS2=P1*(1+2*RR*k*T*COS0)ENDFUNCTIONFUNCTIONPS3(XA,YA,XM,YM,XN,YN,P1,P2,r0,H0)implicitnoneREALPS3,XA,YA,XM,YM,XN,YN,P1,P2,r0,H0,AN,AM,MN,PS2AN=SQRT(XN-XA)*2+(YN-YA)*2)AM=SQRT(XM-XA)*2+(YM-YA)*2)MN=SQRT(XN-XM)*2+(YN-YM)*2)PS3=ABS(AN*PS2(XA,YA,XM,YM,P1,P2,r0,H0)-&AM*PS2(XA,YA,XN,YN,P1,P2,r0,H0)/MNENDFUNCTIONFUNCTIONPS4(XA,YA,XM,YM,XN,YN,X

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽，若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間，僅對(duì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請(qǐng)與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論