《反演与层析成像》实验 - 图文

{ }

b++;

for(c=0;c

for(k=0;k

if(d[k].x>d[k+1].x) { }

temp1=d[k].x; d[k].x=d[k+1].x; d[k+1].x=temp1; temp2=d[k].y; d[k].y=d[k+1].y; d[k+1].y=temp2;

b=0;

while(d[b].x!=100) {

ray[12*i+j].point[b].x=d[b].x;

ray[12*i+j].point[b].y=d[b].y;

fprintf(fp4,\

\\n\

b++;

}

for(a=1;a

mm=sqrt(pow((ray[12*i+j].point[a-1].x-ray[12*i+j].point[a].x),2)+pow((ray[1

2*i+j].point[a-1].y-ray[12*i+j].point[a].y),2)) ;

if(ray[12*i+j].n>=0)

ray[12*i+j].transform[(int)(ray[12*i+j].point[a-1].y/3)][(int)(ray[12*i+j].point[a-1].x/5)]=mm;

if(ray[12*i+j].n<0)

ray[12*i+j].transform[(int)(ray[12*i+j].point[a].y/3)][(int)(ray[12*i+j].point[a-1].x/5)]=mm;

} }

}

for(i=0;i<144;i++) { }

fp3=fopen(\矩阵A.txt\ for(i=0;i<144;i++)

{ }

for(j=0;j<12;j++) {

for(y=0;y<9;y++)

for(j=0;j<12;j++) { }

for(y=0;y<9;y++) { }

ray[i].time+=ray[i].transform[j][y]*xc[j][y];

fprintf(fp3,\

fprintf(fp3,\}

fprintf(fp3,\

fp5=fopen(\矩阵A2.txt\ for(i=0;i<144;i++) { }

fp6=fopen(\走时.txt\ for(i=0;i<144;i++) {

for(j=0;j<12;j++) { }

fprintf(fp5,\

for(y=0;y<9;y++)

fprintf(fp5,\

} }

fprintf(fp6,\

三、 结果图

实验六 射线追踪-向后处理

一、 实验原理

在向后进行反向追踪时，也要考虑来自各个方向的射线。找出来自把该接收点围起来的各个边界射线中，到达该接收点旅行时最小的射线，确定此射线来自的边界以及边界上的插值点，把该插值点作为下一个接收点，以此类推，至到接收点为炮点，就完成了一个接收点到炮点的射线追踪。在此过程中，要记录每条射线经过每个格子的长度。

实验七~八 反演方法

一、 实验原理

计算理论走时

利用每条射线经过每个格子的长度以及经过格子的慢度，便可求得每条线的理论走时。

上述的改进虽然可以解决原LTI算得的节点旅行时不是最小的问题，却降低了计算效率。在一些复杂的模型里，它的效率会是原LTI的1/5甚至更低。所以为了兼顾精度和效率的要求，可以设定循环次数。另外，在比较粗糙的初始模型初始迭代可以用原LTI来进行射线追踪，到后面反演得到的模型精度比较高之后，再用改进LTI。

此外，在实际工程当中，每炮和它的接收点对应的射线只会覆盖整个区域很小的一部分，所以为了提高效率，在对每一炮进行射线追踪时，可以只计算这些射线覆盖的区域，而不必计算整个区域。这个区域可以由炮点最左边的接收点和最右边的接收点来确定，考虑到有回波的可能，此区域可适当地扩大。

下面是计算循环次数不相同得到的射线追踪图。原模型是在坐标为（1500m ，300m）到(4000m,300m),以及从(0,700m)到（1500m,700m）有一水平的高速反射界面。在（1500m,300m）到（1500m,700m）有一坚直的相同速度的高速反射界面，水平反射界面以上是从地表初速度800m/s以每增深10m,速度增大30m/s。炮点在原点，图上所标的S的位置。每个接点相隔100m,第一个接收点在炮点右侧200m处。

(a) 循环0次