帮我写一个随机生成三维多孔介质结构的代码

创建一个随机生成三维多孔介质结构的代码，我们可以使用Python的NumPy库来生成。这里提供一个简单的示例，使用正则网格和随机分布的方式来创建一个模拟的多孔介质：

import numpy as np

def generate_porous_media(size=(100, 100, 100), porosity=0.2):
    # 定义介质大小和孔隙率
    volume_size = size
    voxel_size = 1.0  # 单位假设为米
    total_volume = np.prod(volume_size)

    # 创建均匀的立方体网格
    grid = np.indices(volume_size).astype(np.float32)
    grid -= np.array(grid.shape) / 2  # 将中心移动到原点

    # 计算每个像素是否为空洞
    mask = np.random.uniform(size=volume_size) < porosity

    # 把空洞区域填充为0，实心部分填充为1
    porous_grid = 1 - mask * (1 - porosity / total_volume)

    return porous_grid

# 示例：生成一个100x100x100的多孔介质，孔隙率为0.2
porous_structure = generate_porous_media()

这个代码生成了一个指定尺寸的三维数组，其中随机分布的孔隙被标记为0，实心部分为1。你可以通过调整`size`参数改变体积大小，`porosity`参数改变孔隙率以模拟不同的材料特性。