Question

我正在使用MODFLOW-2000运行地面沉降模型。但是，沉降文件的输出为二进制数据。因为我正在为模型做数百种场景，所以有什么方法可以使用python脚本将其转换为文本。

Answer 1

SUB包的二进制输出与MODFLOW二进制头文件具有相同的格式。您需要知道写入二进制文件的输出文本字符串的名称。有关SUB软件包的信息，请参见MODFLOW-2005 online documentation中的表1，以确定给定的SUB软件包二进制文件的文本字符串。

下面显示了如何使用Z DISPLACEMENT和flopy将二进制沉降文件中的numpy数据转换为ascii文件：

import numpy as np
import flopy

# open the binary file
sobj = flopy.utils.HeadFile('model.zdisplacement.bin', 
                            text='Z DISPLACEMENT')

# get all of the available times in the file
times = sobj.get_times()

# extract the data for the last time in the file
zd = sobj.get_data(totim=times[-1])

# save the z-displacement for the first layer (layer 0) to an ascii file
# zd is a 3D numpy array with a shape of (nlay, nrow, ncol)
np.savetxt('layer0.zdisplacement.txt', zd[0])

如果一层以上，则需要保存每一层的数据。

您可以使用以下命令输出文件中的所有数据：

for t in sobj.get_times():
    zd = sobj.get_data(totim=t)
    for k in range(nlay):
        fpth = 'layer{}_{}.zdisplacement.txt'.format(k, t)
        np.savetxt(fpth, zd[k])

使用FloPy从MODFLOW-2000二进制输出中提取沉降数据

1 个答案: