地球物理层析成像四边激发四边接收研究思路及程序源代码.docxVIP

层析成像实验报告实验课程 层析成像学院名称 地球物理学院 专业名称 勘查技术与工程 学生姓名 廖征 学生学号 201205060105 指导教师 曹俊兴教授 二〇一五年 九月程序说明：本次试验程序设计是根据模型模拟地震层析成像，采用单边放炮和三边放炮两种模型，程序设计思路是根据每条炮射线都会与每条竖线相交，然后如果和第n个接收点相交，就会生成n个点，然后还会和中间一些横轴相交，产生m个交点，然后把这n+m个点组成数组的横坐标排序，然后根据两点之间距离公式求出相邻元素两点间的距离，从而得出正演结果。如上图中，对黄色线分析，此线段会与表格有13个点，通过两条直线的交点可以求出每个点的具体坐标，同时对每个点标号，然后对相应的点的很坐标排序:,,,…,然后根据新的序列，求出每相邻两个点的距离并求出没小段距离传播的时间，在对线段所对应的方块进行赋值，将每个方块的旅行时间得出，由此每条射线都会有12*12个方块的旅行时间，由此做成正演。程序流程图：四边激发，四边接收正演源程序：/************************************************************************************************************//*******************************层析成像（四边激发，四边接收）***********************************************//****************************************正演************************************************************//************************************************************************************************************/#include stdio.h#include stdlib.h#include math.h#define PI 3.1415926#define E 2.71828main(){ int i,j,m,n,p,q,a11,b11; //设置计数参数 double v[12][9]={0.0}; //设计速度模型 double a1[4]={0.0}; //设置激发点数组 double b1[4]={0.0}; double c1[4]={0.0}; double d1[4]={0.0}; double a2[9]={0.0}; //设置接收点数组 double b2[9]={0.0}; double c2[12]={0.0}; double d2[12]={0.0}; double z=0; /***********************对接收点激发点赋初值*************************/ for(i=0;i4;i++) { a1[i]=7.5+10*i; } for(i=0;i4;i++) { b1[i]=7.5+10*i; } for(i=0;i4;i++) { c1[i]=7.5+6*i; } for(i=0;i4;i++) { d1[i]=7.5+6*i; } for(i=0;i9;i++) { a2[i]=2.5+5*i; } for(i=0;i9;i++) { b2[i]=2.5+5*i; //printf(%lf\n,b2[i]); } for(i=0;i12;i++) { c2[i]=1.5+3*i; } for(i=0;i12;i++) { d2[i]=1.5+3*i; } /****************************设置交点坐标数组*********************************/ double abx[4][9][20]={0.0}; double aby[4][9][20]={0.0}; double acx[4][9][20]={0.0}; double acy[4][9][20]={0.0}; double adx[4][9][20]={0.0}; double ady[4][9][20]={0.0}; double bax[4][9][20]={0.0}; double bay[4][9][20]={0.0}