资源描述：
《复杂形体重磁异常正演》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、直立长方体模型如图2-1建立三维直角坐标系，x轴指向北方向，y轴指向东方向，z轴铅垂向下。在xy平面下存在一规模为2a*2b*(h2-h1)直立长方体，其在各坐标轴上的投影起始坐标为：ξ1、ξ2--x坐标的起点和终点η1、η2--x坐标的起点和终点ζ1、ζ2--x坐标的起点和终点，即h1、h2如下图所示：图2-1直立长方体直立长方体的重力异常正演公式直立长方体引起的重力异常表达式为：Δgr=Gσ{ξ•ln[η+r]+η•ln[ξ+r]-ζ•arctan[ξηζr]}∣x1x2∣y1y2∣z1z2在上式中，采用CGSM制，重力
2、异常Δgr单位为Gal，长度单位为cm，剩余密度σ单位为gcm3，引力常量G为6.672×10-8cm3(g.s2)。直立长方体的磁力异常正演公式直立长方体引起的重力异常表达式为：ΔT(x,y,z）=μ04πM{k1·ln[r+(ξ-x)]+k2·ln[r+(η-y)]+k3·ln[r+(ζ-z)]+k4·arctanξ-x（η-y）（ξ-x）2+rζ-z+（ζ-z）2+k5·arctanξ-x（η-y）（η-y）2+rζ-z+（ζ-z）2+k6·arctanξ-x（η-y）rζ-z}∣x1x2∣y1y2∣z1z2其中，L0
3、、M0、N0及α、β、γ分别为地磁场及总磁化强度的方向余弦；则有：k1=M0·γ+N0·β，k2=L0·γ+N0·α，k3=L0·β+M0·α，k4=L0·α，k5=M0·β，k6=-N0·γ。同样，在上式中，采用CGSM制，导磁系数μ0=1。输入文件格式1.1输入参数文件格式设计输入参数文件用于记录场的分量标识、导数的阶数、观测面的形态标识、场源参数文件名、测点文件名、输出文件名等信息。本次实验使用的一个参数文件样例cmd.dat为：请选择需要计算的场分量：1重力或磁力异常2重力或磁力异常的导数1请输入需要计算的导数阶次：
4、1请选择观测面的形态：1平面规则网2平面非规则网3曲面规则网4曲面非规则网1场源参数文件名：source.dat测点位置文件名（规则网时为grd文件，非规则网为xyz文件）：BXYZL.GRD重力异常输出文件：anomaly_gra.GRD磁力异常输出文件（倾斜磁化异常和化极磁异常）：anomaly_mag.GRDanomaly_pol.GRD1.2场源参数文件格式设计本次实验的场源参数保存在“source.dat”中。第一列为剩余密度(g/cm3)；第二列为磁化强度(10-6CGSM)；第三列为磁化方向倾角(DEG)；第四
5、列为磁化方向与x轴的夹角(DEG)；第五列为磁化方向与y轴的夹角(DEG)；第六列～第七列为x坐标的起点和终点(km)；第八列～第九列为y坐标的起点和终点(km)；第十列～第十一列为z坐标的起点和终点(km，向下为正)。如下所示：0.7,20000,50,85,5,-9,-3,4,8,2,70.8,34000,50,85,5,0,7,0,5,3,70.9,17000,50,85,5,-4,8,-8,-3,2,71.3计算点输入坐标格式设计程序能对如下四种测网进行正演计算：1平面规则网2平面非规则网3曲面规则网4曲面非规则网其
6、中，规则网计算点坐标保存在grd文件中文件后缀为grd；非规则网计算点坐标保存在xyz文件中，文件后缀为txt。对于平面测网，采用只读取第一点高程数据的办法来节省读取外部文件的时间。注：测点坐标文件的后缀必须严格按照如上规定，否则程序将报错并要求重新修改参数文件。输出文件格式计算完成后程序将自动输出重力异常、磁力异常以及化极磁力异常到外部文件，规则网数据输出到grd文件中，非规则网数据输出到xyz(txt)文件中。PROGRAMComplex_body_forwardIMPLICITNONECHARACTER*80cmdfi
7、le,sourcefile,stationfile,output_grafile,&output_magfile,output_polfileINTEGERcomponent,order,plane!场的分量标识，导数的阶数，观测面的形态标识INTEGERn_source,n_station,i,j,mpoint,line,cs,cxyzREALxmin,xmax,ymin,ymaxREAL,ALLOCATABLE::SOURCE(:,:),XYZ(:,:)REAL,ALLOCATABLE::ANOM_GRA(:),ANOM
8、_MAG(:),ANOM_POL(:)cmdfile='cmd.dat'cs=11cxyz=3CALLINPUT_CMDFILE(component,order,plane,cmdfile,sourcefile,&stationfile,output_grafile,output_magf
