《反射地震记录道》由会员分享,可在线阅读,更多相关《反射地震记录道(7页珍藏版)》请在金锄头文库上搜索。
1、反射地震记录道实验报告学年学期:20172018学年第1学期院 别:土木工程学院专 业:勘查技术与工程专业方向:岩土工程方向班 级:勘查1501学生:余庭垚学号:3150604120指导教师:江凡实验报告内容一、实验要求1、利用 VB 进行褶积过程的编程。2、设计两个地层,而地层的反射系数序列和地震子波的褶积结果就是反射记 录道。二、实验目的1、增强利用 VB 解决现实问题的能力,更加熟练的使用 VB。2、在使用 VB 进行褶积编程的过程中,进一步熟悉了解地层反射系数序列和 地震子波。3、在编写代码的过程中对简单的褶积过程更加了解,认识到该模型可制作正 演理论地震记录。三、实验过程Dim l1
2、(1000) As DoubleDim l2(1000) As DoubleDim n As IntegerDim v1 As DoubleDim v2 As DoubleDim v3 As DoubleDim r1 As DoubleDim r2 As DoubleDim r3 As DoubleDim fs(1000) As DoubleDim r(100) As DoubleDim f As DoubleDim dt As DoubleDim ltw As IntegerDim gt As DoublePrivate Sub Command1_Click()f = Text1.Textd
3、t = 0.01ltw = 50For t = 0 To ltwr(t)二 20 * Exp(3.14 * f * dt * (t 一 Itw / 2)八 2) * (1 一 2 * (3.14 * f * dt * (t - Itw / 2)八 2)Next tPicture3.ScaIe (一20, 0)一(20, 50)For t = 0 To Itw 一 1Pic ture3.Li ne (r( t), t)-(r(t + 1), t + 1), vbRedNext tPicture1.AutoRedraw = TruePicture1.ScaIe (0, 0)一(99, 99)For
4、 i = 0 To 99I1(i) = 0.2 * i + 20I2(i) = 80Next iFor i = 0 To 99Pic ture1.Li ne (i, l1(i)-(i + 1, l1(i + 1)Pic ture1.Li ne (i, l2(i)-(i + 1, l2(i + 1)Next iv1 = 1000v2 = 2000v3 = 3000r1 = 1.1r2 = 1.2r3 = 1.3gt1 = (v2 * r2 一 v1 * r1) / (v2 * rou2 + v1 * r1)gt2 = (v3 * r3 一 v2 * r2) / (v3 * rou3 + v2 *
5、 r2)n = 150Call Convolution(r, ltw, fs, n)Picture2.AutoRedraw = TruePicture2.Scale (一15, 0)一(1000, 120)For i = 0 To 99 Step 2For j = 0 To nfs(j) = 0Next jfs(Int(l1(i) / v1 * 1000) + 1) = gt1fs(Int(l2(i) / v2 * 2000) + 1) = gt2Call Convolution(r, ltw, fs, n)For j = 0 To n - 1Picture2.Line (fs(j) + i
6、* 15, j)-(fs(j + 1) + i * 15, j + 1), vbBlue Next jNext iEnd SubPrivate Sub Command2_Click()Text1.Text = Picture1.ClsPicture2.ClsPicture3.ClsEnd SubFunctionConvolution(x()As Double,lenX As Integer,y() As Double,lenYAs Integer)Dim x2(5000) As DoubleDim y2(5000) As DoubleDim z(5000) As DoubleFor i = 0
7、 To lenY - 1y2(i) = y(lenY - i - 1)Next iFor i = 0 To lenX - 1x2(i + lenY) = x(i)Next iFor i = 0 To lenX + lenY - 2For j = 0 To lenY - 1z(i) = z(i) + x2(i + j + 1) * y2(j)Next jNext iDim k As Integerk = Int(lenX / 2)For i = 0 To lenY - 1y(i) = z(i + k)Next iEnd Function四、实验分析Forml清楚的分辨出两个地层,这也是我们进行褶积的原因,可以很清晰的分辨出两个地层。五、实验结论由图像可得,固定一个dt值,改变频率f,地震子波随之改变,波长变短,其褶积最后图像也变得越来越窄,变得越来越尖锐,可作为地震记录。而频率越高的,图像中所显示的地层界面越精确,越能精确地反映出地层信息。