魏建 孫祥娥

摘要：基于泰勒展開的旅行時公式是計算P波旅行時的一種常用方法。泰勒展開的結(jié)果僅在計算點附近準(zhǔn)確，而地層中具有垂直對稱軸的橫向各向同性（transverse isotropy with vertical symmetry axis， VTI）介質(zhì)進一步放大了旅行時計算的不準(zhǔn)確性。本文提出了一種改進的P波旅行時計算方法。該方法以基于泰勒展開的偏移距八階旅行時方程為基礎(chǔ)，利用系數(shù)匹配方法分別處理偏移距八階項與六階項的系數(shù)、偏移距六階項與四階項的系數(shù)，用整合處理優(yōu)化得到的參數(shù)結(jié)構(gòu)，形成了基于系數(shù)匹配與參數(shù)整合的偏移距四階P波旅行時計算方法。水平層狀VTI介質(zhì)模型測試結(jié)果表明，該方法能夠得到更小的計算誤差，保證其在遠偏移距處具有較好的表現(xiàn)。三維垂直地震剖面數(shù)據(jù)實驗結(jié)果表明，該方法經(jīng)參數(shù)調(diào)整后能夠進行數(shù)據(jù)分析，在中偏移距處的表現(xiàn)相對穩(wěn)定，在遠偏移距處具有一定優(yōu)勢。

關(guān)鍵詞：旅行時計算；系數(shù)匹配；參數(shù)整合；VTI介質(zhì)

doi：10.13278/j.cnki.jjuese.20220309

中圖分類號：P631

文獻標(biāo)志碼：A

收稿日期：2022-11-10

作者簡介：魏建（1991-），男，博士研究生，主要從事油氣信息探測與儀器裝備方面的研究，E-mail： 201873033@yangtzeu.edu.cn

通信作者：孫祥娥（1970-），女，教授，博士生導(dǎo)師，主要從事地震信號處理方面的研究，E-mail： sxeyangtzeu@163.com

基金項目：國家自然科學(xué)基金項目（51978078）

Supported by the National Natural Science Foundation of China （51978078）

P-Wave Traveltime Calculation Equation Based on Coefficient

Matching and Integration ProcessingWei Jian， Sun Xiange

College of Electronic Information， Yangtze University， Jingzhou 434023， Hubei， China

Abstract： The traveltime formula based on Taylor expansion is a common method to calculate P-wave traveltime. The Taylor expansion is accurate only near the calculation point， and the inaccuracy of the traveltime approximation is further magnified in the medium of?? transverse isotropy? with vertical? symmetry axis?? （VTI）. An improved method for calculating P-wave traveltime is presented in this paper. This method utilizes the 8th-order offset-traveltime equation based on the Taylor expansion， and the coefficient matching method is used to determine the coefficients of the series expansion. Then the integration processing is used to optimize the parameters. Finally， the 4th-order P-wave traveltime calculation method is formed? based on? coefficient matching and parameter integration. The test results for horizontal layered VTI medium models show that the proposed method has relatively small numerical errors and relatively good performance at the far offset. The experimental results for 3D VSP data show that after parameter adjustment， the method effectively analyzes the data， exhibiting stable performance at middle offsets and a certain advantage at far offsets.

Key words： traveltime calculation; coefficient matching; parameter integration; VTI medium

0 引言

旅行時計算公式是地球物理學(xué)的重要研究課題之一，其計算精度影響著速度分析、時差校正、成像和正向建模等關(guān)鍵處理步驟^[1]。為了在具有垂直對稱軸的橫向各向同性（transverse isotropy with vertical symmetry axis， VTI）介質(zhì)中獲得更準(zhǔn)確的P波旅行時，P波旅行時方程已由雙曲線形式逐漸轉(zhuǎn)變?yōu)榉请p曲線形式。

在非雙曲線形式的旅行時公式中，基于泰勒級數(shù)展開的P波近似公式是最基本的一種計算方法，如：Castle^[2]通過改進高階旅行時公式提出的平移雙曲線近似法，Alkhalifah等^[3]利用聲學(xué)近似理論提出的三參數(shù)偏移距四階旅行時計算方法，Ursin等^[4]使用一次系數(shù)匹配方法提出的偏移距四階旅行時計算方法，以及Ursin等^[5]隨后提出的改進的幾何擴展近似方法。近期，Xu等^[6]提出了基于偏移距泰勒級數(shù)展開與聲學(xué)近似的計算方法。基本的計算方法受限于偏移距高階項的使用，影響了其在介質(zhì)中計算的準(zhǔn)確性。除了基于泰勒級數(shù)展開的基本旅行時計算方法之外，其他能夠提高旅行時計算精度的方法相繼被提出。Fomel^[7]通過研究相速度與群速度的關(guān)系，提出了基于橢圓近似的旅行時計算方法。為了提高擾動方法的準(zhǔn)確性，隨后提出了優(yōu)化后的擾動方法^[8-9]。隨著研究的深入，F(xiàn)omel等^[10]提出了基本廣義時差近似法，在一定條件下與基于橢圓近似方法具有相同形式。針對基本廣義時差方程存在計算參數(shù)不對稱及準(zhǔn)確性不高的情況，提出了通過增加參與計算參數(shù)以及射線數(shù)目的改進方法，有：修正的廣義時差近似法^[11]、擴展的廣義時差近似法^[1]、基于三射線的廣義時差近似法^[12]以及改進的廣義時差近似法^[13]。之后，Ravve等^[14]提出了基于五參數(shù)的旅行時計算方法。Abedi^[15]基于有理近似方法提出了多階旅行時計算方程，隨后與Pardo通過計算無限偏移距附近的漸近級數(shù)構(gòu)建了新的旅行時計算方法^[16]?；诼眯袝r計算相關(guān)理論的快速發(fā)展，進一步改善旅行時計算準(zhǔn)確性的方法將會被提出并廣泛應(yīng)用。

在VTI介質(zhì)中，一般泰勒級數(shù)展開的傳統(tǒng)旅行時方程的計算會變得不準(zhǔn)確，如何改善這種情況一直是研究的重點。基于此，本文以一般泰勒級數(shù)P波旅行時方程為基礎(chǔ)，將兩次系數(shù)匹配方法與參數(shù)整合處理相結(jié)合，提出了適用于VTI介質(zhì)的旅行時計算新方法，并通過模型和實際地震數(shù)據(jù)的測試與分析證明該方法的優(yōu)勢，為旅行時計算研究提供新思路。

1 理論方法

1.1 偏移距四階旅行時計算公式

基于泰勒級數(shù)展開的旅行時與偏移距之間的一般P波近似形式可表示^{[4， 16-17]}為

其中，

式中：t₀為雙程旅行時；x為偏移距；v_n為正常時差（normal moveout， NMO）速度；S₂為非均質(zhì)性參數(shù)，與各向異性參數(shù)η^[7]的關(guān)系為S₂=1+8η。

Ursin等^[4-5]為了提高P波旅行時計算的準(zhǔn)確性，根據(jù)Alkhalifah等^[3]提出的計算公式，利用方程（1）的偏移距六階項提出一種偏移距四階旅行時計算公式^[4]：

其中，

式中，S₃為非均勻性參數(shù)。

由于S₃不能直接從地震數(shù)據(jù)中估計，所以參數(shù)S₃與S₂的近似關(guān)系可以通過垂直慢度平方公式與聲學(xué)近似理論^[5]得出：

1.2 基于系數(shù)匹配與參數(shù)整合處理的旅行時計算方法

在方程（3）的推導(dǎo)過程中，Ursin等^[4-5]利用一次系數(shù)匹配法來處理偏移距的六階項，將方程的最高階近似為一般的偏移距四次方項。按照這個思路，可以做進一步處理。為了避免推導(dǎo)更為復(fù)雜的偏移距系數(shù)和非均勻性參數(shù)的近似關(guān)系，同時不局限于泰勒級數(shù)旅行時方程的高階近似形式，我們以式（3）引入六次方項的處理方式為基礎(chǔ)，目前僅將式（1）中偏移距的階數(shù)提高到八次方來嘗試提高最終旅行時的計算準(zhǔn)確性。若保持偏移距的最高階數(shù)與方程（1）相同，那么系數(shù)c₂、c₃和c₄可以利用兩次系數(shù)匹配方法進行處理。

首先，將方程（1）寫為偏移距最高階為八次方項的截斷形式：

其中：

根據(jù)方程（3）的形式，在方程（6）的基礎(chǔ)上利用一次系數(shù)匹配方法來處理系數(shù)c₃和c₄，并省略式（6）中偏移距的八階項，可以得到：

其中，

式中，S₄為非均勻性參數(shù)。

在得到式（8）后，按照同樣的處理思路能夠?qū)懗鲎罱K的表達式：

式（8）中引入了參數(shù)b₁，在一定程度上并不適合直接近似處理式（8）中的偏移距六次方項。基于此，可以繼續(xù)使用系數(shù)匹配方法處理系數(shù)c₂和c₃，得到間接參數(shù)b₂：

為了優(yōu)化參數(shù)b₁、b₂在旅行時方程應(yīng)用中的結(jié)構(gòu)，可在式（11）的基礎(chǔ)上通過參數(shù)整合得到參數(shù)B的新形式，即

B=b₁·b^-1₂。??? （12）

對于式（12），根據(jù)垂直慢度平方公式與聲學(xué)近似理論^[5]能夠得到非均勻性參數(shù)S₄與S₂的近似關(guān)系：

然后，將式（5）（13）代入到式（9）（11）中，再利用S²₃近似替換S³₂來區(qū)別直接簡化的形式，即可得到化簡后的參數(shù)b₁、b₂：

最后，將式（12）（14）帶入式（10）中，就能夠得到基于兩次系數(shù)匹配與參數(shù)整合（twice coefficient matching and integration， TCMI）的P波旅行時平方公式：

2 實驗結(jié)果

2.1 模型結(jié)果對比

文中使用深度為4 km、偏移距為10 km的水平層狀模型（圖1），具體參數(shù)如表1所示。

誤差實驗的對比方法為：基于兩次系數(shù)匹配與參數(shù)整合的方法、擾動理論近似方法^[8-9]、擴展的廣義時差近似方法^[1]以及三階有理近似方法^[15]，分別記為TCMI、PT（perturbation theory）、EGMA（extended generalized moveout approximation）、RA（rational approximation）。驗證4種方法的評價指標(biāo)為參數(shù)方程^[1]計算出的P波準(zhǔn)確時間。

其中：

式中：x（p）為參數(shù)方程計算出的偏移距；t（p）為參數(shù)方程計算出的P波準(zhǔn)確時間；p為射線參數(shù)；t_0，i、v_n，i、η_i分別為每層的旅行時、NMO速度與各向異性參數(shù)^[16]。計算誤差^[16]為

式中，t_x為P波旅行時計算值。

2.1.1 均勻?qū)訝頥TI介質(zhì)模型

對于均勻?qū)訝頥TI介質(zhì)模型，通過式（18）計算得出的誤差曲線如圖2所示。在整體偏移距上，PT方法與RA方法的計算誤差較大，與PT、RA方法相比，TCMI方法具更小的計算誤差。當(dāng)偏移距為0～8.6 km時，TCMI方法的誤差小于EGMA方法；當(dāng)偏移距為8.6 km時，TCMI方法與EGMA方法的誤差曲線趨向重合；當(dāng)偏移距大于8.7 km時，EGMA方法的誤差為4種方法中最小。在遠偏移距10 km處：PT、RA、TCMI和EGMA方法的誤差分別為2.34%、1.63%、1.33%、0.95%；與PT、RA方法相比，TCMI方法的計算精度分別提高了1.01%、0.30%，相比于EGMA方法，計算精度降低了0.38%。綜上可知：在均勻?qū)訝罡飨虍愋越橘|(zhì)中，TCMI方法在偏移距0～8.6 km范圍內(nèi)具有較小的誤差，在偏移距8.6～10 km范圍內(nèi)的誤差僅小于EGMA方法。

2.1.2 非均勻?qū)訝頥TI介質(zhì)模型

非均勻?qū)訝頥TI介質(zhì)模型的誤差曲線如圖3所示。在圖3中：TCMI方法與其他3種方法的計算誤差相比仍然較??；由于RA方法與PT方法計算誤差的差值較為接近，會出現(xiàn)視覺上的重合，實際上PT方法的誤差更大。在遠偏移距10 km處：PT、RA、EGMA和TCMI方法的誤差分別為5.74%、5.73%、5.03%、1.83%；與PT、RA和EGMA方法相比，TCMI方法的準(zhǔn)確性分別提高了3.91%、3.90%、3.20%。綜上可知：在非均勻VTI介質(zhì)中，TCMI方法在整體以及遠偏移距處均具有更小的計算誤差。

2.2 實際地震數(shù)據(jù)結(jié)果對比

VTI介質(zhì)是一種能夠近似表示成具有水平層狀的周期性薄互層介質(zhì)，在實際地震數(shù)據(jù)中較為常見。實驗對比使用的某研究區(qū)三維垂直地震剖面（3D vertical seismic profile， 3D VSP）數(shù)據(jù)同樣具有類似性質(zhì)，其檢波器共160級，觀測深度為275.84～2 700.00 m。

在實際地震數(shù)據(jù)驗證中，使用孫祥娥等^[18]提出的時差對比方法，即將3D VSP的初至?xí)r間作為衡量標(biāo)準(zhǔn)，通過時差判斷旅行時公式的表現(xiàn)。時差分析方法如下：

t_d=t_x-t_fb。??? （19）

式中：t_d為時差；t_fb為初至?xí)r間。

2.2.1 方法適用性

針對實際數(shù)據(jù)的研究重點，以深度1 784.76 m的實驗結(jié)果為例，初至?xí)r間如圖4所示。TCMI方法通過式（19）計算出的時差結(jié)果如圖5所示。由圖5可知：偏移距1 500 m之后的時差并未保持相對水平，而是出現(xiàn)向下彎曲的情況。

基于彎曲情況以及分析TCMI方法的組成，發(fā)現(xiàn)：該情況是由整合處理得到的參數(shù)B（式（12））造成的。式（12）的提出主要基于兩點考慮：一是Ursin等^[4-5]提出的分式結(jié)構(gòu)的思想；二是以模型結(jié)果為導(dǎo)向。為了使TCMI方法更好地應(yīng)用于實際數(shù)據(jù)，基于以上兩點考慮，此時應(yīng)動態(tài)調(diào)整參數(shù)B的組合形式：

B=b₂·b^-1₁。??? （20）

此時，TCMI方法計算出的時差如圖6所示。由圖6可看出，調(diào)整參數(shù)后的時差結(jié)果相對水平，與圖5差異明顯;說明調(diào)整參數(shù)后的TCMI方法能夠應(yīng)用于實際地震數(shù)據(jù)中。

基于模型與實際數(shù)據(jù)的結(jié)果，此時可將參數(shù)B稱為動態(tài)參數(shù)，這為后續(xù)將TCMI方法優(yōu)化為自適應(yīng)的形式奠定了基礎(chǔ)。

2.2.2 方法對比

EGMA、RA、PT方法和調(diào)整參數(shù)后的TCMI方法通過式（19）計算出的時差對比結(jié)果如圖7所示。由圖7可看出：TCMI方法與經(jīng)典方法，即EGMA、RA、PT方法，計算出的時差均在一定范圍內(nèi)，說明調(diào)整參數(shù)后TCMI方法有效。

為了區(qū)分TCMI方法與其他3種方法的差別，我們計算了EGMA、RA、PT方法時差與TCMI方法時差的差值，結(jié)果如圖8所示。分析圖8可知：TCMI方法與EGMA、RA、PT方法的時差之差均較小，但在中偏移距500～2 500 m范圍內(nèi)，因TCMI方法具有較大的時差計算值，所以EGMA、RA、PT方法與TCMI方法時差差值曲線會出現(xiàn)小于0的情況；而在遠偏移距，即大于2 500 m時，TCMI方法的計算值較小，可知該方法在該范圍具有優(yōu)勢。對于其他3種方法，在中偏移距范圍內(nèi)，EGMA方法與RA方法的差值較小，均與PT方法的差值較大；在遠偏移距范圍內(nèi)，EGMA方法的優(yōu)勢大于RA和PT方法，PT方法的表現(xiàn)一般。

3 討論

本文在利用含有八次方項的旅行時公式推導(dǎo)出一種新的旅行時計算方法。在此基礎(chǔ)上，還可以將偏移距的階數(shù)提高多少是需要思考的^[17]，是否仍能使用類似的處理方法也是需要著重研究的。在計算P波旅行時的時候，還可以通過計算模型空間中每一點的到時來分析每一點到時的誤差分布。在實際數(shù)據(jù)的實驗中，我們提出了與模型試驗中形式不同的參數(shù)B。這兩種形式使得TCMI方法更加豐富且完善，為后續(xù)將其優(yōu)化為自適應(yīng)的計算方法奠定了基礎(chǔ)。另外，與射線追蹤方法相比，旅行時公式仍存在一定不足，所以，研究旅行時公式應(yīng)集中在它的應(yīng)用上，如相關(guān)參數(shù)的反演^[19]或成像^[20-22]。使用哪種公式進行分析需要重點分析，對于TCMI方法，還可以嘗試聯(lián)合其他旅行時公式同時對數(shù)據(jù)進行分析，這是一個可以深入研究的方向。

4 結(jié)論

該研究基于一般的偏移距八階旅行時計算公式，利用系數(shù)匹配方法得到了處理高階項后的參數(shù)b₁、b₂，然后利用參數(shù)整合、垂直慢度與聲學(xué)近似理論，結(jié)合實際地震數(shù)據(jù)實驗得到了參數(shù)B的兩種形式，形成了基于兩次系數(shù)匹配與參數(shù)整合的P波旅行時計算方法，得到以下結(jié)論：

1）對于各向異性參數(shù)相同的層狀VTI介質(zhì)，與擴展的廣義時差近似方法、三階有理近似方法相比，基于兩次系數(shù)匹配與參數(shù)整合的P波旅行時計算方法在遠偏移距10 km處的計算精度分別提高了1.01%、0.30%；相比于擴展的廣義時差近似方法，計算精度降低了0.38%。對于各向異性參數(shù)不同的VTI介質(zhì)，與擴展的廣義時差近似方法、三階有理近似方法和擴展的廣義時差近似方法相比，準(zhǔn)確性分別提高了3.91%、3.90%、3.20%。

2）對于3D VSP數(shù)據(jù)，通過調(diào)整參數(shù)B的形式，使基于兩次系數(shù)匹配與參數(shù)整合的P波旅行時計算方法適用于實際數(shù)據(jù)分析。通過實驗對比，基于兩次系數(shù)匹配與參數(shù)整合的P波旅行時計算方法在遠偏移距范圍內(nèi)具有一定優(yōu)勢。

3）基于兩次系數(shù)匹配與參數(shù)整合的P波旅行時計算方法在模型與實際數(shù)據(jù)實驗中的表現(xiàn)說明，系數(shù)匹配與參數(shù)整合處理可以在提高旅行時計算準(zhǔn)確性上起到重要作用，但在一定程度上也會影響方法的計算穩(wěn)定性。

4）結(jié)合模型與實際數(shù)據(jù)實驗中的表現(xiàn)可知：基于兩次系數(shù)匹配與參數(shù)整合的P波旅行時計算方法可以成為分析旅行時的一種選擇。如何提高該方法的穩(wěn)定性，以及將參數(shù)B的兩種形式更好地與方法結(jié)合中是下一步的研究重點。

參考文獻（References）：

[1]Abedi M M， Stovas A. Extended Generalized Nonhyperbolic Moveout Approximation [J]. Geophysical Journal International，2019，26（2）：1428-1440.

[2]Castle R J. A Theory of Normal Moveout[J]. Geophysics， 1994， 34（6）： 859-881.

[3]Alkhalifah T， Tsvankin I. Velocity Analysis for Transversely Isotropic Media [J]. Geophysics，1995，60（5）：1550-1566.

[4]Ursin B， Stovas A. Traveltime Approximations for a Layered Transversely Isotropic Medium[J]. Geophysics， 2006， 71（2）： D23-D33.

[5]Ursin B， Stovas A. Improved Geometrical Spreading Approximation in Layered Transversely Isotropic Media [J]. Geophysics， 2009， 74（5）： D85-D95.

[6]Xu S， Stovas A， Mikada H. Traveltime and Relative Geometric Spreading Approximation in Elastic Orthorhombic Medium [J]. Geophysical Journal International， 2020， 85（5）： C153-C162.

[7]Fomel S. On Anelliptic Approximations for qP Velocities in VTI Media [J]. Geophysics， 2004， 52（3）： 247-259.

[8]Xu S， Stovas A， Hao Q. Perturbation-Based Moveout Approximations in Anisotropic Media[J]. Geophysical Prospecting， 2017， 65（5）： 1218-1230.

[9]Abesdi M M， Riahi M A， Stovas A. Three-Parameter NMO Correction in Layered Anisotropic Media： A Stretch-Free Approach [J]. Geophysics， 2019， 84（3）： C129-C142.

[10]Fomel S， Stovas S A. The Generalized Non-Hyperbolic Moveout Approximation [J]. Geophysics， 2010， 75（2）： U9-U18.

[11]Stovas A， Fomel S. The Modified Generalized Moveout Approximation： A New Parameter Selection [J]. Geophysics， 2017， 65（3）： 687-695.

[12]Abedi M M， Riahi M A， Stovas A. A New Parameterization for Generalized Moveout Approximation Based on Three Rays [J]. Geophysical Prospecting， 2019， 67（5）： 1243-1255.

[13]Jin S， Stovas A， Hu X. Improved Generalized Moveout Approximation with a Novel Parameterization [J]. Geophysics， 2022， 87（3）： C49-C61.

[14]Ravve I， Koren Z. Traveltime Approximation in Vertical Transversely Isotropic Layered Media [J]. Geophysical Prospecting， 2017， 65（6）：1559-1581.

[15]Abedi M M. Rational Approximation of P-Wave Kinematics： Part 1： Transversely Isotropic Media [J]. Geophysics， 2020， 85（5）： C175-C185.

[16]Abedi M M， Pardo D. Large-Offset P-Wave Traveltime in Layered Transversely Isotropic Media [J]. Geophysics， 2021， 86（3）： C65-C74.

[17]魏建，孫祥娥. VTI介質(zhì)中基于降階補償處理的旅行時計算方法[J]. 山東科技大學(xué)學(xué)報（自然科學(xué)版），2022，41（1）：86-91.

Wei Jian， Sun Xiange. A Method of Traveltime Calculation Based on Descending Order Compensation Processing in VTI Medium[J]. Journal of Shandong University of Science and Technology （Natural Science）， 2022， 41（1）： 86-91.

[18]孫祥娥，凌云，高軍，等. 各向異性假設(shè)條件下三維 VSP 與地面地震數(shù)據(jù)旅行時差分析及參數(shù)求取[J]. 石油地球物理勘探，2009，44（6）：720-725.

Sun Xiange， Ling Yun， Gao Jun， et al. Travel Time-Difference Analysis and Parameter Extraction for 3D VSP Data and Surface Seismic Data in Anisotropic Media [J]. Oil Geophysical Prospecting， 2009， 44（6）： 720-725.

[19]張利振，孫成禹，王志農(nóng)，等. 面波信息約束的初至波走時層析反演方法[J]. 物探與化探，2023，47（5）：1198-1205.

Zhang Lizhen， Sun Chengyu， Wang Zhinong， et al. First Arrival Wave Travel Time based Tomography Inversion with Surface Wave Information as Constraints[J]. Geophysical and Geochemical Exploration， 2023， 47（5）：1198-1205.

[20]張志鵬，孫章慶，韓復(fù)興，等. 影響高斯束偏移成像質(zhì)量和效率的主要因素分析[J]. 世界地質(zhì)，2022，41（2）：349-364.

Zhang Zhipeng， Sun Zhangqiang， Han Fuxing， et al. Analysis of Main Factors Affecting Imaging Quality and Efficiency of Gaussian Beam Migration[J]. World Geology， 2022， 41（2）： 349-364.

[21]李洪麗，劉財，田有，等. 中國東北地殼結(jié)構(gòu)的地震層析成像[J]. 吉林大學(xué)學(xué)報（地球科學(xué)版），2022，52（1）：270-280.

Li Hongli， Liu Cai， Tian You， et al. Seismic Tomography of the Crustal Structure in Northeast China[J]. Journal of Jilin University （Earth Science Edition）， 2022， 52（1）： 270-280.

[22]萬曉杰，鞏向博，成橋，等. 基于DeCNN的逆時偏移低頻噪聲壓制方法[J].吉林大學(xué)學(xué)報（地球科學(xué)版），2023，53（5）：1593-1601.

Wan Xiaojie， Gong Xiangbo， Cheng Qiao， et al. Eliminating Low-Frequency Noise in Reverse-Time Migration Based on DACNN[J]. Journal of Jilin University （Earth Science Edition）， 2023， 53（5）： 1593-1601.