2024华数杯国际数学建模竞赛A题：放射性废水

最新推荐文章于 2024-06-19 11:46:05 发布

蟹堡王首席大厨

最新推荐文章于 2024-06-19 11:46:05 发布

阅读量1k

点赞数 20

文章标签：数学建模

本文链接：https://blog.csdn.net/qq_41489047/article/details/135663275

版权

2024华数杯国际数学建模竞赛A题（代码+思路）持续更新中

问题重述

2011年3⽉，⽇本东海岸发⽣了地震，引发了福岛第⼀核电站事故，导致三个核反应堆熔毁，并在⼀场巨⼤海啸中冲毁了电站的冷却系统，核燃料产⽣融化的碎⽚。为了冷却熔化的核燃料，海⽔不断注⼊反应堆，导致⼤量放射性核素污染的冷却⽔。尽管全球各国⼈⺠反对，⽇本政府于2023年8⽉24⽇开始强制排放经过处理的福岛放射性废⽔到太平洋。受核素污染的放射性废⽔总量超过100万吨。整个项⽬预计将⾄少持续30年。附录是⽇本政府公布的四轮排放计划。
在这里插入图片描述
这些核废⽔含有氚，⼀种可以在环境中存在很⻓时间的放射性同位素。放射性元素的污染程度是指环境中存在的放射性元素的数量以及对⼈类和⽣态系统的潜在危害程度。通常通过测量放射性元素的浓度、辐射⽔平、半衰期等参数来评估。氚的⾼放射性使其具有通过辐射损害细胞和组织的潜⼒。在海洋环境中，氚将被⽣物吸收并进⼊⻝物链，导致对⽣态系统中的物种造成辐射损害，影响海洋⽣物的繁殖和⽣态平衡。放射性废⽔在海⽔中的扩散路径受到许多因素的影响，包括⽔流、海床地形、⽔深、潮汐和季节变化以及环境条件。了解放射性废⽔在环境中的传输和扩散可以帮助我们评估对周围海洋⽣态系统和⼈类健康的巨⼤影响。通过建⽴数学模型和分析，我们可以预测核废⽔的扩散范围和路径，制定环境保护措施和应急计划。

1. 建⽴扩散数学模型，描述海⽔中放射性废⽔扩散的速率和⽅向，考虑⽔流、环境条件和其他影响因素。已知截⾄2023年8⽉27⽇上午12:00，从⽇本排放的放射性废⽔量为1095吨。如果不再排放放射性废⽔，请预测到2023年9⽉27⽇时在⽇本附近海域的放射性废⽔污染范围和程度。

问题⼀要求建⽴⼀个数学模型来描述海⽔中放射性废⽔的扩散情况。这个问题不能简单的套⽤机器学习模型，应该使⽤机理分析的⽅法来进⾏建模。我们可以考虑使⽤扩散⽅程（Diffusion Equation）来描述放射性物质在⽔体中的传播，扩散⽅程可以采⽤⼀维或⼆维的形式，取决于具体情况。
在这里插入图片描述

import numpy as np
import matplotlib.pyplot as plt

def simulate_diffusion(L, T, D, delta_x, delta_t):
    # 模型参数
    num_points = int(L / delta_x) + 1
    num_steps = int(T / delta_t) + 1

    # 网格和初始条件
    x = np.linspace(0, L, num_points)
    C = np.zeros((num_points, num_steps))

    # 设置初始条件（瞬时释放）
    C[:, 0] = np.where((x >= L/2 - 5) & (x <= L/2 + 5), 1, 0)

    # 数值模拟
    for n in range(0, num_steps - 1):
        for i in range(1, num_points - 1):
            C[i, n+1] = C[i, n] + D * delta_t / delta_x**2 * (C[i+1, n] - 2 * C[i, n] + C[i-1, n])

    return x, C

# 模拟参数
L_simulation = 200  # 海域长度
T_simulation = 50  # 模拟总时间
D_simulation = 0.1  # 扩散系数
delta_x_simulation = 2  # 空间步长
delta_t_simulation = 0.5  # 时间步长

# 运行模拟
x_result, C_result = simulate_diffusion(L_simulation, T_simulation, D_simulation, delta_x_simulation, delta_t_simulation)

# 绘制结果
plt.figure(figsize=(12, 6))
for n in range(0, C_result.shape[1], int(C_result.shape[1] / 10)):
    plt.plot(x_result, C_result[:, n], label=f'Time: {n * delta_t_simulation:.1f} hours')

plt.title('Radioactive Wastewater Diffusion in Seawater')
plt.xlabel('Distance (km)')
plt.ylabel('Concentration')
plt.legend()
plt.grid(True)
plt.savefig("./第一题.png")
plt.show()

在这里插入图片描述

2. 在2023年，⽇本政府已经三次排放了放射性废⽔。如果将来不再排放，请建⽴数学模型研究三次排放后放射性废⽔的扩散路径。考虑海洋环流模式、⽔流动⼒学、海床地形、⽔深变化、潮汐影响和季节波动等因素。估计污染中国领海需要多⻓时间。

在这里插入图片描述

import numpy as np
import matplotlib.pyplot as plt
# 模型参数
D = 0.01 # 扩散系数
L = 100 # 空间⻓度
T = 30 # 模拟的总时间
Nx = Ny = Nz = 100 # 空间⽹格数
Nt = 300 # 时间步数
# 空间和时间步⻓
dx = dy = dz = L / Nx
dt = T / Nt
# 初始化浓度场
C = np.zeros((Nx, Ny, Nz, Nt+1))
# 设置初始条件
C[:, :, :, 0] = 0.0 # 初始浓度为零
# 数值求解
for t in range(Nt):
 for x in range(1, Nx-1):
 for y in range(1, Ny-1):
 for z in range(1, Nz-1):
 C[x, y, z, t+1] = C[x, y, z, t] + D * dt / dx**2 * (C[x+1,
y, z, t] - 2*C[x, y, z, t] + C[x-1, y, z, t]) + \
 D * dt / dy**2 * (C[x, y+1, z, t] - 2*C
[x, y, z, t] + C[x, y-1, z, t]) + \
 D * dt / dz**2 * (C[x, y, z+1, t] - 2*C
[x, y, z, t] + C[x, y, z-1, t])
# 绘制结果（路径追踪）
fig = plt.figure()
ax = fig.add_subplot(projection='3d')
ax.scatter(xs=path_x, ys=path_y, zs=path_z, c='r', marker='o', label='Pat
h')
ax.set_xlabel('空间位置 (x)')
ax.set_ylabel('空间位置 (y)')
ax.set_zlabel('空间位置 (z)')
ax.set_title('放射性废⽔三维扩散模拟及路径追踪')
plt.legend()
plt.show()