地球物理课程设计

1、课程设计程序及成果图件报告姓名：学号：专业：地球物理学指导老师： 2014.1.7前言：重磁异常正演是在给定地质体的形状大小、空间位置及物性参数的条件下，求它在它外部空间任意点上的物理场值。电子计算机的应用使正演方法进入了一个新的阶段，它可以对任意形状物体进行计算；对物性“非均匀”的地质体，也可以通过将其分解成若干个“均匀”物性的物体的方法和分区块均匀的数值方法来计算，使正演方法可以达到相当高的精度。以下是我对截面为多边形的水平柱体模型的重磁异常的正演程序。C语言代码： #include "stdio.h"#include "math.h"#define

2、 pi 3.14159 #define n 101 /*测点数*/#define N 360 /*拟合多边形边数*/double fun(double a,double b) /*关于atan(b/a)的讨论*/ doublethet; if(fabs(a)>1.0e-15) if(b/a)>0.0) if(a>0.0&&b>0.0)thet=atan(b/a);else if(a<0.0&&b<0.0)thet=atan(b/a)-pi; else if(b/a)<0.0)thet=atan(b/a)+pi; else

3、 if(fabs(a)<=1.0e-15) if(b<0.0)thet=-pi/2;else if(b>=0.0)thet=pi/2; return(thet); void main() FILE *fp;fp=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

4、,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;i<N;i+) x1i=x0+r*sin(alpha*i);z1i=z0-r*cos(alpha*i); for(k=0;k<n;k+) xkk=20.0*k;for(i=0;i<N-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

5、+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+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(

6、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)/(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);deltag

7、k=-2*G*sigma*gk; for(i=0;i<n;i+) printf("%lf %lf %lfn",deltagi,deltaxi,deltazi); fprintf(fp,"%lf %lf %lfn",deltagi,deltaxi,deltazi); fclose(fp);数据结果：G X Z0.011686 7.067300 7.0673500.011924 7.062859 7.3543370.012167 7.049233 7.6497860.012415 7.025691 7.9536400.012667 6.991471 8

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.382523 10.6639110.014851 6.231952 11.0321550.015142 6.062661 11.4050340.015437 5.873697 11.7817270.015736 5.664118 12.1613120

9、.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 14.7983110.018192 3.144038 15.1542070.018501 2.713665 15.4995370.018808 2.256373 15.8322420.019113 1.772162 16.1501690.01

10、9416 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 17.7683580.021381 -3.004993 17.8846130.021629 -3.697400 17.9653790.021867 -4.404783 18.0086530.022094 -5.124601 18.01260

11、10.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 -10.255658 16.8342670.023384 -10.960629 16.4878160.023469 -11.647843 16.0978720.023536 -12.313736 15.6656380.023584 -12.95

12、4885 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 11.6308850.023377 -16.504704 10.9437840.023272 -16.851025 10.2389090.023150 -17.152950 9.5198990.023010 -17.409984 8.7903

13、630.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 -17.980408 3.6823880.021362 -17.899392 2.9902320.021103 -17.782879 2.3153580.020835 -17.633012 1.6598150.020559 -17.45203

14、5 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.2686850.018478 -15.511596 -2.7257010.018168 -15.165991 -3.1558000.017858 -14.809824 -3.5591480.017548 -14.445052 -3.9360540

15、.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