盖革只由到时确定地震震中的概率方法.docVIP

只由到时确定地震震中的概率方法 盖革（L. C. Geiger），1910 [译者说明：译文根据Geiger,L.，1912, Probability method for the determination of earthquake epicenters from the arrival time only，1912， Bulletin St. Louis University 8, 60-71. (由F. W. L. Peebles 和 A. H. Corey根据Geiger, L.， 1910，Herdbestimmung bei Erdbeben aus den Ankunftszeiten. Nachrichten von der K?niglichen Gesellschaft der Wissenschaften zu G?ttingen, Mathematisch-Physikalische Klasse, 331-349. 从德文翻译成英文)。 在上述英文译文中没有给出德文原文所附的一些表格,故本译文中在个别文字上作了修改。译文中的公式标号按原英文本的样式，文章中间有10个公式的标号在公式前面。而在之前还有10个公式的标号在公式的后面。但连续阅读不会造成混淆。读者如发现本译文中的错误，请发邮件给 zhzhao@seis.ac.cn 赵仲和] 符号：λ=经度，φ=纬度，t=发震时刻，以时、分、秒计。T=波的走时，以分、秒计。*表示相伴的值是假定的，是要被校正的。△是在地球表面上的角距离。?只使用一次，表示观测值，当附在时间上时，表示到达一台站的观测到时。数字下标和通用下标n表示台站数据，或者是假定值或者是绝对值。下标0表示震中的函数。 为使用这一方法，我们必须首先假定一个近似的震中位置，这一假定值越接近最终位置，则校正值越小，定位结果会越准确，但最初假定值有几度的误差将可以计算出相当好的位置。这一过程在很大程度上依赖于第一个初至波的走时，当然，它不是距离的线性函数。 如果我们以任意方便的方法假定了、 和，并假定每个值中有误差，我们用、 和表示，于是有式（1）： （1） 在这些等式中 ，左手边的量是要寻找的量，右手边的第一项是没有改变的假定量，它们受第二项的校正。于是问题本身变成确定δ量。为确定这些δ量，我们进行下述过程： 首先利用球面余弦定律计算从假定的震中到每个台站的距离： (2) 在本讨论最后的制表工作中，这是用取对数的方法计算的，在余弦定律右边的两项用I和II 表示，这两项要相加。在这样得到 之后，由自然余弦表得到相应的以角度测量的。 现在我们有了从台站n到一个假定震中位置的近似角距离或者说弧距离，于是我们从一个走时表得到近似的走时。根据我们的符号，这一走时称作。近似的走时加上近似的发震时刻便给出近似的到时，即 (3) 如果在我们的假定位置中没有误差，在记录的时间中没有误差，而且观测中也没有误差的话，那么，*tn将是我们的观测到时?tn。但实际上存在误差，所以在二者之间将有差值Fn，即 Fn=*tn-?tn (4) 由此可清楚地看到，如果没有、 和这些误差，而且?tn是准确的，那么就不会存在Fn。 我们的下一步是建立Fn和假定由于Fn的存在而造成的误差之间的关系。 在此时我们假定，除了经度，关于震中位置和时间的其他数据都是正确的，并假定由于经度差造成了Fn，于是 就是说，Fn将是一个时间校正值，等于每单位经度的时间校正值乘以经度校正的单位数。在当前考虑经度变化的情况下，如果我们检查每度经度的走时变化率，并利用这一比值，以一个等于这一比值乘以实际经度校正值（以相同的单位）的量来校正我们的时间，那么将得到满意的Fn值。现在，将刚才关于经度所说的话以同等效力独立地应用于纬度，于是有。另一方面，如果我们假定发震时刻*t0有一个改变，则到时*tn也必定有一个类似的改变，因些到时随发震时刻的变化率总是为单位1。然而，为了便于数字处理，将这一比值按其他比值那样对待，因此有第三个校正值，尽管该值只不过是。由这一考虑得知，Fn实质上是以差分形式存在的，因为它可以等于前述的差分；因此，当δ项变得无限小时，Fn将趋于以0作为极限。现在，由泰勒定理可以看到，如果Fn分别依赖于、和，则Fn可等于它们之和。在构成这一等式时，用偏微分符号?代替比值中的δ，以指出对每一个量的对待犹如其他量均为常数，并在最后使用总和。在我们下文列出的工作中将会看到，对于纬度和经度，我们严格遵循这一过程，对于时间我们也同样这样做，但事实上时间比值总是等于1，因此只