版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
《電法勘探及地電學(xué)》基礎(chǔ)理論實(shí)驗(yàn)報告學(xué)號:姓名:實(shí)習(xí)項(xiàng)目一中梯法的電阻率正演模擬一、實(shí)習(xí)目的:掌握中梯法球體正演場特征二、要求:設(shè)計球體的空間位置、物性參數(shù)和計算空間根據(jù)中梯法球體的正演公式及設(shè)計空間編寫程序;程序如下:#include<stdio.h>#include<math.h>voidmain(){ FILE*fp; fp=fopen("中間梯度電位.txt","w"); floatL,h,I,p,a,u1,u2,r0; floatx[101][101]={0},y[101][101]={0},z[101][101]={0},z1[101][101]={0}; inti,j; printf("請輸入中間梯度法參數(shù),r0,u2:");scanf("%f%f",&r0,&u2); h=2*r0; for(i=0;i<101;i++) for(j=0;j<101;j++) { x[i][j]=-15.0+j*0.3; y[i][j]=-15.0+i*0.3; z1[i][j]=1+2*(u2-1)/(2*u2+1)*r0*r0*r0*(h*h+y[i][j]*y[i][j]-2*x[i][j]*x[i][j])/pow((h*h+y[i][j]*y[i][j]+x[i][j]*x[i][j]),2.5); fprintf(fp,"%f%f%f\n",x[i][j],y[i][j],z1[i][j]); } fclose(fp); }球的空間位置(0,0,10),球體半徑r0=3,情況一:u2=ps/p1=10得下面電位圖:實(shí)習(xí)項(xiàng)目二二層電測深電阻率正演模擬一、實(shí)習(xí)目的:掌握電測深電阻率正演方法及二層電測深曲線特征二、要求:根據(jù)二層電測深正演公式及設(shè)計參數(shù)編寫程序;設(shè)計參數(shù)表:第一層電阻率506060100250403020505第二層電阻率403020505506060100250第一層厚度246810246810r:{1.5,2,3,4,5,6,9,12,15,20,30,45,60,90}編寫程序如下:#include<stdio.h>#include<math.h>main(){ FILE*fp; fp=fopen("測深數(shù)據(jù).txt","w"); floatr[14]={1.5,2,3,4,5,6,9,12,15,20,30,45,60,90}; floatp1[10]={50,60,60,100,250,40,30,20,50,5}; floatp2[10]={40,30,20,50,5,50,60,60,100,250}; floath[10]={2,4,6,8,10,2,4,6,8,10}; floatps[14]={0};inti,j,a,n=30; floats,sun=0,k; for(i=0;i<10;i++) { k=(p2[i]-p1[i])/(p2[i]+p1[i]); for(a=0;a<14;a++) { for(j=1;j<=30;j++) { s=pow(k,j)*r[a]/pow((r[a]*r[a]+2*j*h[i]*2*j*h[i]),1.5); sun=sun+s; } ps[a]=p1[i]*(1+2*r[a]*r[a]*sun); sun=0; printf("%f%f\n",log10(r[a]/h[i]),log10(ps[a]/p1[i])); fprintf(fp,"%f%f\n",log10(r[a]/h[i]),log10(ps[a]/p1[i])); } printf("\n"); fprintf(fp,"\n"); } fclose(fp); }最終在EXCEL成圖:橫坐標(biāo)為log10(r/h),縱坐標(biāo)為log10(ps/p1)(r=AB/2,h為深度,ps為視電阻率,p1為第一層點(diǎn))實(shí)習(xí)項(xiàng)目3電偶極子電位正演一、實(shí)習(xí)目的:掌握電偶極子正演場特征二、要求:設(shè)計電偶極子的空間位置、物性參數(shù)和計算空間根據(jù)電偶極子正演公式及設(shè)計空間編寫程序;編寫程序如下:#include<stdio.h>#include<math.h>voidmain(){ FILE*fp; fp=fopen("偶極子電位.txt","w"); floatL,h,I,p,a; floatx[11][11]={0},y[11][11]={0},z[11][11]={0}; inti,j; printf("此程序用于求偶極子電流在地面的電位和電位梯度,設(shè)偶極子中心為原點(diǎn),長度為L,中心到地面距離為h,電流大小為I,電阻率為p,偶極子傾斜角度a\n "); printf("請分別給出參數(shù),長度,距離,電流大小,電阻率,角度:\n"); scanf("%f%f%f%f%f",&L,&h,&I,&p,&a); for(i=0;i<11;i++) for(j=0;j<11;j++) { x[i][j]=-10.0+j*2; y[i][j]=-10.0+i*2; z[i][j]=(I*p/(2*3.1415))*(1/pow(((x[i][j]+L*cos(a))*(x[i][j]+L*cos(a))+y[i][j]*y[i][j]+(h+L*sin(a))*(h+L*sin(a))),0.5)-1/pow(((x[i][j]-L*cos(a))*(x[i][j]-L*cos(a))+y[i][j]*
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 航運(yùn)物流企業(yè)市場拓展策略考核試卷
- 金屬標(biāo)牌的快速模具設(shè)計與制造考核試卷
- 眼鏡生產(chǎn)過程中的環(huán)境保護(hù)與節(jié)能減排考核試卷
- 服裝批發(fā)商品牌差異化競爭策略考核試卷
- 物流設(shè)備管理與實(shí)踐考核試卷
- 露天游樂場所設(shè)備抗風(fēng)性與穩(wěn)定性分析考核試卷
- 新車銷售展示廳布局與設(shè)計考核試卷
- 2024年度湖南省安全員之C證(專職安全員)能力提升試卷B卷附答案
- 2024年度甘肅省安全員之A證(企業(yè)負(fù)責(zé)人)全真模擬考試試卷A卷含答案
- 法人授權(quán)委托合同(2篇)
- 第4課《公民的基本權(quán)利和義務(wù)》(課件)-部編版道德與法治六年級上冊
- 二十四節(jié)氣(教學(xué)設(shè)計)-2023-2024學(xué)年六年級下冊綜合實(shí)踐活動長春版
- 備用金管理制度(規(guī)章制度)
- 精裝修施工圖的深化設(shè)計管理辦法
- 農(nóng)業(yè)資源生態(tài)保護(hù)2022 年耕地輪作項(xiàng)目績效評價報告
- 國開(甘肅)2024年春《地域文化(專)》形考任務(wù)1-4終考答案
- 《中國慢性阻塞性肺疾病基層診療與管理指南(2024年)》解讀
- 應(yīng)用PDCA提高入院宣教的知曉率
- 加油站服務(wù)意識培訓(xùn)PPT課件.ppt
- 汽輪機(jī)各輔助設(shè)備的作用
- 工程交接證書
評論
0/150
提交評論