我认为这是可能的4台是的,而不是三个。
(大约)均匀的假设速度我们得到以下问题。假设假设v_p美元和v_s美元地震和你所有的接收器之间是常数,产量收益率以下系统:\{对齐}开始d & = \√6{\δx ^ 2 + \δy ^ 2 + \δz ^ 2} ={3已知3未知数}\文本,文本\四\{为}\左\{{数组}{rl} \ \开始δx & =间{地震}\文本间\文本{接收机}\ \ \δy & = y_ \{地震}- y_ \文本{接收机}\ \ \δz & = z_ \文本{地震}- z_ \文本{接收机}\{数组}{分裂}\ \开始结束结束{分裂}\。,\ \ \δt & = t_s - t_p ={测量}\文本,\ \ p & = \离开(\压裂{1}{v_s} - \压裂{1}{v_p} \右)= \文本未知{1}。结束\{对齐}重申你的问题然后方程给出了提出问题:$ $ \δt = p d = p \√6{\δx ^ 2 + \δy ^ 2 + \δz ^ 2}。$ $$ $ \δt = p d。$ $例如,对于v_p = 2000美元米/秒,v_s = 1200美元m / s,我们发现p = 1/3000美元和可以计算如d = 600美元那$ \δt = 600/3000 = 0.2美元年代。
现在,我们必须设置=逆问题。你通常至少需要尽可能多的方程未知数。在这种情况下,我们有4个未知数(速度和(x, y, z)美元地震位置)。这意味着我们需要4(独特的)方程。我认为最简单的方法获得4独特的方程结合测量:$ $ p {pmatrix} d_1开始\ \ \ d_2结束\ \ d_3 \ \ d_4 \ {pmatrix} = {pmatrix} \ \开始δδt_2 t_1 \ \ \ \ \ \δδt_4 t_3 \ \ \ \ {pmatrix}, $ $表示与美元间的{1,2,3,4}$和美元y_ {1, 2, 3, 4} $接收器的纬度和经度,假设测量表面,我们得到:$ $ p {pmatrix} \ \开始√6{(间\文本{地震}-x_1) ^ 2 + (y_ \文本{地震}-y_1) ^ 2 + z_文本{地震}^ 2}\ \ \ \√6{(间\文本{地震}-x_2) ^ 2 + (y_ \文本{地震}-y_2) ^ 2 + z_文本{地震}^ 2}\ \ \ \√6{(间\文本{地震}-x_3) ^ 2 + (y_ \文本{地震}-y_3) ^ 2 + z_文本{地震}^ 2}\ \ \ \√6{(间\文本{地震}-x_4) ^ 2 + (y_ \文本{地震}-y_4) ^ 2 + z_文本{地震}^ 2}\ \结束{pmatrix} = {pmatrix} \ \开始δδt_2 t_1 \ \ \ \ \ \δδt_4 t_3 \ \ \ \ {pmatrix}结束。$ $
最后,我们必须解决逆问题。不幸的是,该系统是非线性(未知的“隐藏”在广场和不能被分离的根源)。你必须使用一些非线性规划求解。例如,在MATLAB可以做到这一点
%真实地震位置x_earthquake = 0;y_earthquake = 0;z_earthquake = 1000;%真实站位置x_1 = 1000;y_1 = 1000;x_2 = 500;y_2 = -300;x_3 = -400;y_3 = -100;x_4 = -10; y_4 = 800; X = [x_1;x_2;x_3;x_4]; Y = [y_1;y_2;y_3;y_4]; % True velocity structure p0 = 1/3000; % Forward equation (F(1)=x_eartquake, F(2)=y_earthquake, F(3)=z_earthquake, F(4)=p. t =@(F) F(4) * sqrt( (F(1)-X).^2 + (F(2)-Y).^2 + F(3).^2); % True recordings measured_times = t([x_earthquake,y_earthquake,z_earthquake,p0]); % Misfit function (=0 at optimum) G =@(F) t(F) - measured_times; % Invert选项%的选项= optimoptions (“fsolve”、“FiniteDifferenceType”,“中央”);%使用这些将在MATLAB给出更好的结果F_inv% F_inv= fsolve (G(1000、1000、1000、0.1)选项)F_inv = fsolve (G, [1000、1000、1000、0.1])x_earthquake_inv = F_inv (1) y_earthquake_inv = F_inv (2) z_earthquake_inv = F_inv (3) p_inv = F_inv (4)
结果如果你把这个复制到https://octave-online.net/例如,它将找到一个地震在(0,0,1000)和p = 1/3000(与,一些小的错误经过12位数)由于非线性求解器进行求解。