论文-空洞探测模型的建立与求解.docVIP

空洞探测模型的建立与求解 摘要 　本文论述在探测过程中, 对山体、坝体、隧洞等某些内部定位问题。通过构造出一个测定空洞位置的代数模型，并以矩阵形式来确定存在的空洞路径，然后由路径的交点来确定空洞的位置。 关键词 代数模型 矩阵 空洞路径 一、问题的重述 山体，隧洞，坝体等的某些内部结构可以用弹性波测量来确定。一个简化问题可描述为一块均匀成的矩形平板内有一些充满空气的空洞。在平板的两个邻边分别等距地设置若干波源，在它们的对边对等地安放同样多的接收器，记录弹性波由每个波源到达对边上每个接收器的时间。根据弹性波在介质中和在空气中不同的传播速度，来确定板内空洞的位置。现考察如下具体问题： 一块240m×240m的平板 ,在AB边等距地设置7个波源Ｐi （ｉ＝１,２，…，７）边上对等地安放７个接收器Ｑj（ｊ＝１，２，…，７），记录由Ｐi发出的弹性波到达Ｑj 的时间t ij 秒 在AD边等距地设置7个波源Ri i 1,…,7 ,BC边对等地安放7个接收器S j j 1,…,7 ,记录由Ri 发出的弹性波到达S j的时间r ij 秒 。 已知弹性波在介质和空气中的传播速度分别为2880（米/秒）和320（米/秒），且弹性波活板边缘的传播速度与在介质中的传播速度相同。 表1 　波源Pi 发出的弹性波到达接收器Qj 的时间t ij s tij Q1 Q2 Q3 Q4 Q5 Q6 Q7 P1 0.0611 0.0895 0.1996 0.2032 0.4181 0.4923 0.5646 P2 0.0989 0.0592 0.4413 0.4318 0.4770 0.5242 0.3805 P3 0.3052 0.4131 0.0598 0.4153 0.4156 0.3563 0.1919 P4 0.3221 0.4453 0.4040 0.0738 0.1789 0.0740 0.2122 P5 0.3490 0.4529 0.2263 0.1917 0.0839 0.1768 0.1810 P6 0.3807 0.3177 0.2364 0.3064 0.2217 0.0939 0.1031 P7 0.4311 0.3397 0.3566 0.1954 0.0760 0.0688 0.1042 表2 　波源Ri 发出的弹性波到达接收器S j 的时间rij s rij S1 S2 S3 S4 S5 S6 S7 R1 0.0645 0.0602 0.0813 0.3516 0.3867 0.4314 0.5721 R2 0.0753 0.0700 0.2852 0.4341 0.3491 0.4800 0.4980 R3 0.3456 0.3205 0.0974 0.4093 0.4240 0.4540 0.3112 R4 0.3655 0.3289 0.4247 0.1007 0.3249 0.2134 0.1017 R5 0.3165 0.2409 0.3214 0.3256 0.0904 0.1874 0.2130 R6 0.2749 0.3891 0.5895 0.3016 0.2058 0.0841 0.0706 R7 0.4434 0.4919 0.3904 0.0786 0.0709 0.0914 0.0583 据上提出两个问题： 问题一: 确定该平板内空洞的大概位置。 问题二：只根据由Pi发出的弹性波到达Qj的时间tij i,j 1,…,7 ，能确定空洞的位置吗；讨论在同样能够确定空洞位置的前提下，减少波源和接受器的方法。 二、问题的分析 对于空洞探测问题，图中R和P表示波源，Q和S表示接收器，而四边是边长为240m的正方形平板。需要探测在平板中是否存在空洞，可根据以下原理： 由于弹性波在介质和在空洞中的传播速度是不同的，所以可以根据在介质中测的传播时间来判断在某个波源到接收器这段距离中是否存在空洞。然后可以把问题所测时间列为一个时间矩阵。因此，我们首先通过波在边界上传播所用时间的计算值与其最小观测值的差，求出所允许误差的最大值。具体所求出的所用时间为0.0833秒，而数据中的最小值为0.0583，所以其差值为0.025.为其所允许的最大误差。 由已知条件可以算出PI到QJ的距离,得出一个距离矩阵.根距离矩阵,算出一个在没有空洞的前提下波由PI到QJ的时间距阵 由于该模板为正形,由对称可以得出横波也有这样的时间矩阵.则可以分别由测定的横波和纵波传播的实际时间矩阵减去没有空洞时算出的时间矩阵,得出两个时间差距阵.由时间差距阵内元素的正负可判定从这条波上是否存在空洞.从横纵两个方向上的有空洞波的路径交点可得出空洞的大体位置。以及讨论在能够确定空洞位置的前提下,