地球物理課程設(shè)計(jì)

1、課程設(shè)計(jì)程序及成果圖件報(bào)告姓名：學(xué)號(hào)：專業(yè)：地球物理學(xué)指導(dǎo)老師：前言：重磁異常正演是在給定地質(zhì)體的形狀大小、空間位置及物性參數(shù)的條件下，求它在它外部空間任意點(diǎn)上的物理場(chǎng)值。電子計(jì)算機(jī)的應(yīng)用使正演方法進(jìn)入了一個(gè)新的階段，它可以對(duì)任意形狀物體進(jìn)行計(jì)算；對(duì)物性“非均勻”的地質(zhì)體，也可以通過(guò)將其分解成若干個(gè)“均勻”物性的物體的方法和分區(qū)塊均勻的數(shù)值方法來(lái)計(jì)算，使正演方法可以達(dá)到相當(dāng)高的精度。以下是我對(duì)截面為多邊形的水平柱體模型的重磁異常的正演程序。C語(yǔ)言代碼： #include stdio.h#include math.h#define pi 3.14159 #define n 101 /*測(cè)點(diǎn)數(shù)*/

2、#define N 360 /*擬合多邊形邊數(shù)*/double fun(double a,double b) /*關(guān)于atan(b/a)的討論*/ doublethet; if(fabs(a)1.0e-15) if(b/a)0.0) if(a0.0&b0.0)thet=atan(b/a);else if(a0.0&b0.0)thet=atan(b/a)-pi; else if(b/a)0.0)thet=atan(b/a)+pi; else if(fabs(a)=1.0e-15) if(b=0.0)thet=pi/2; return(thet); void main() FILE *fp;fp=

3、fopen(輸出.txt,w+); double xN=0.0,zN=0.0,Pn=0.0,Qn=0.0; deltagn=0.0,deltaxn=0.0,deltazn=0.0; double M=2000.0,G=6.67e-5,sigma=2.67,r=200.0,x0=1000.0,z0=1000.0; inti,k; doublexk101,x1N,z1N; double gn=0.0,alpha=(2*pi)/N; doubleMx=M*sin(pi/4),Mz=M*cos(pi/4);for(i=0;iN;i+) x1i=x0+r*sin(alpha*i);z1i=z0-r*co

4、s(alpha*i); for(k=0;kn;k+) xkk=20.0*k;for(i=0;iN-1;i+) xi=x1i-xkk;xi+1=x1i+1-xkk;zi=z1i;zi+1=z1i+1;Qk+=(zi+1-zi)*(xi-xi+1)/(pow(zi+1-zi),2)+pow(xi-xi+1),2)*(fun(xi,zi)-fun(xi+1,zi+1)-(1/2)*pow(zi+1-zi),2)*log(pow(xi+1,2)+pow(zi+1,2)/(pow(xi,2)+pow(zi,2)/(pow(zi+1-zi),2)+pow(xi-xi+1),2); Pk+=(pow(zi+

5、1-zi),2)/(pow(zi+1-zi),2)+pow(xi-xi+1),2)*(fun(xi,zi)-fun(xi+1,zi+1)-(1/2)*(zi+1-zi)*(xi-xi+1)*log(pow(xi+1,2)+pow(zi+1,2)/(pow(xi,2)+pow(zi,2)/(pow(zi+1-zi),2)+pow(xi-xi+1),2);deltaxk=(Mx*Pk+Mz*Qk)/(2*pi);deltazk=(Mx*Qk-Mz*Pk)/(2*pi);gk+=(1/2)*(zi+1-zi)*(xi*zi+1-xi+1*zi)*log(pow(xi+1,2)+pow(zi+1,2)

6、/(pow(xi,2)+pow(zi,2)/(pow(zi+1-zi),2)+pow(xi-xi+1),2)-(xi+1-xi)*(xi*zi+1-xi+1*zi)/(pow(zi+1-zi),2)+pow(xi-xi+1),2)*(fun(xi,zi)-fun(xi+1,zi+1);deltagk=-2*G*sigma*gk; for(i=0;in;i+) printf(%lf %lf %lfn,deltagi,deltaxi,deltazi); fprintf(fp,%lf %lf %lfn,deltagi,deltaxi,deltazi); fclose(fp);數(shù)據(jù)結(jié)果：G X Z0.

7、011686 7.067300 7.0673500.011924 7.062859 7.3543370.012167 7.049233 7.6497860.012415 7.025691 7.9536400.012667 6.991471 8.2658010.012925 6.945778 8.5861210.013186 6.887785 8.9144030.013453 6.816638 9.2503890.013724 6.731454 9.5937600.014000 6.631327 9.9441260.014279 6.515330 10.3010240.014563 6.3825

8、23 10.6639110.014851 6.231952 11.0321550.015142 6.062661 11.4050340.015437 5.873697 11.7817270.015736 5.664118 12.1613120.016037 5.433000 12.5427580.016340 5.179451 12.9249210.016646 4.902621 13.3065440.016954 4.601714 13.6862490.017263 4.276001 14.0625410.017572 3.924832 14.4338060.017882 3.547657

9、14.7983110.018192 3.144038 15.1542070.018501 2.713665 15.4995370.018808 2.256373 15.8322420.019113 1.772162 16.1501690.019416 1.261207 16.4510860.019714 0.723882 16.7326920.020008 0.160769 16.9926400.020297 -0.427326 17.2285520.020579 -1.039362 17.4380450.020855 -1.674053 17.6187530.021122 -2.329861

10、 17.7683580.021381 -3.004993 17.8846130.021629 -3.697400 17.9653790.021867 -4.404783 18.0086530.022094 -5.124601 18.0126010.022308 -5.854081 17.9755920.022508 -6.590239 17.8962250.022694 -7.329901 17.7733630.022865 -8.070195 17.6056940.023020 -8.806829 17.3935070.023159 -9.536516 17.1363240.023280 -

11、10.255658 16.8342670.023384 -10.960629 16.4878160.023469 -11.647843 16.0978720.023536 -12.313736 15.6656380.023584 -12.954885 15.1927260.023612 -13.568017 14.6810740.023621 -14.150072 14.1329370.023611 -14.698194 13.5508990.023581 -15.209814 12.9377990.023532 -15.682675 12.2967040.023464 -16.114835

12、11.6308850.023377 -16.504704 10.9437840.023272 -16.851025 10.2389090.023150 -17.152950 9.5198990.023010 -17.409984 8.7903630.022854 -17.622006 8.0538950.022682 -17.789493 7.3137840.022494 -17.912159 6.5743190.022293 -17.991316 5.8383710.022078 -18.028104 5.1091140.021851 -18.023924 4.3895290.021612

13、-17.980408 3.6823880.021362 -17.899392 2.9902320.021103 -17.782879 2.3153580.020835 -17.633012 1.6598150.020559 -17.452035 1.0253920.020276 -17.242270 0.4136290.019987 -17.006083 -0.1741910.019692 -16.745858 -0.7370270.019393 -16.463975 -1.2740740.019091 -16.162781 -1.7847510.018785 -15.844576 -2.26

14、86850.018478 -15.511596 -2.7257010.018168 -15.165991 -3.1558000.017858 -14.809824 -3.5591480.017548 -14.445052 -3.9360540.017238 -14.073522 -4.2869580.016929 -13.696968 -4.6124110.016622 -13.317006 -4.9130620.016316 -12.935131 -5.1896390.016012 -12.552720 -5.4429400.015711 -12.171032 -5.6738150.015413 -11.791209 -5.8831570.015118 -11.414283 -6.0718890.014826 -11.041176 -6.2409520.014539 -10.672710 -6.3913010.014255 -10.309607 -6.5238920.013975 -9.952497 -6.6396770.013700 -9.601925 -6.7395990.013429 -9.258353 -6.8245820.013162 -8.922171 -6.8955340.012901 -8.593699 -6.9533370.012643 -