地震反射道.docVIP

工程物探VB课程报告 学年学期： 2016 ~ 2017学年第1学期 院 别： 土木工程学院 专 业： 勘查技术与工程 专业方向： 岩土方向 班 级： 勘查1401 学 生： 任霄瑶 学 号： 3140604120 指导教师： 江凡 课设目的 认识学习地震反射记录道 二 、课设内容 解决思路： 当地下介质近似为层状介质时，并存在大量的薄层，这表明每个层界面都可以产生反射并成为一个反射面。由于在层界面的两侧速度和密度存在一定的差异，地震波到达层界面之后会出现反射、折射等现象，因此可以计算得到不为零的反射系数，由此可见，一个实际地震记录道就是由这些无数多个子波组成的复合震动。显然，振动的幅值大小与界面的反射系数成正比。利用地震子波与反射系数做卷积，就得到了单道地震记录道。在地面上接收到N个反射波总和的一个地震道记录X(t)可表示为一个反射记录道是地层反射系数序列和地震子波的褶积（卷积）结果。于是反射道可表示为： 2.求每个反射界面的反射系数，即求得反射系数序列 波阻抗z[i]=p[i]*v[i] 反射系数R[i]=(z[i+1]-z[i])/( z[i+1]+z[i] 选取适当的子波. 4.程序说明： Dim layer(100) As Double Dim nt As Integer Dim v1 As Double Dim v2 As Double Dim rou1 As Double Dim rou2 As Double Dim flag As Integer Dim tralen(2000) As Double Dim w(2000) As Double Dim f As Double Dim dt As Double Dim ltw As Integer Dim ref As Double Const PI = 3.1415926535 Private Sub Form_Load() flag = 0 Picture1.AutoRedraw = True Picture1.Scale (0, 0)-(99, 99) Dim oridep As Integer oridep = 50 For i = 0 To 99 layer(i) = oridep Next i For i = 0 To 99 Picture1.Line (i, layer(i))-(i + 1, layer(i + 1)) Next i v1 = 1000# v2 = 2000# rou1 = 1.1 rou2 = 1.2 ref = (v2 * rou2 - v1 * rou1) / (v2 * rou2 + v1 * rou1) f = 100 dt = 0.001 ltw = Int((4# / f / dt)) + 1 Call ricker(w(), f, dt, ltw) ricker Picture3.AutoRedraw = True Picture3.Scale (-20, 0)-(20, 40) For i = 0 To ltw - 1 Picture3.Line (w(i), i)-(w(i + 1), i + 1), vbRed Next i tralen(Int(oridep / v1 * 1000) + 1) = ref nt = 120 Call Convolution(w, ltw, tralen, nt) Picture2.AutoRedraw = True Picture2.Scale (-15, 0)-(1000, 120) For i = 0 To 99 For j = 0 To nt - 1 Picture2.Line (tralen(j) + i * 10, j)-(tralen(j + 1) + i * 10, j + 1), vbBlue Next j Next i End Sub Private Sub Picture1_MouseDown(Button As Integer, Shift As Integer, x As Single, y As Single) flag = 1 Dim pointX As Integer Dim pointY As Integer pointX