p3

p3/README.md

@@ -0,0 +1,229 @@
+# 任务三：月球殖民地水需求与运输分析
+
+## Task 3: Water Supply Analysis for Moon Colony
+
+本目录包含针对 **MCM 2026 Problem B** 任务三的完整分析：
+
+> "Investigate the water needs for a one-year period once the 100,000-person Moon Colony is fully operational. Use your delivery model to understand the additional cost and timeline needed to ensure the colony has sufficient water for one full year after the Moon Colony is inhabited."
+
+---
+
+## 一、水需求模型
+
+### 1.1 模型假设
+
+| 参数 | 符号 | 取值 | 说明 |
+|------|------|------|------|
+| 殖民地人口 | $N$ | 100,000 人 | 题目给定 |
+| 生存用水 | $w_s$ | 2.5 L/人/天 | 饮用+基本代谢 |
+| 卫生用水 | $w_h$ | 0.4×α L/人/天 | α=舒适度系数 |
+| 医疗应急用水 | $w_m$ | 5.0 L/次 | 每次医疗事件 |
+| 每百人日患病率 | $p$ | 2% | 医疗事件频率 |
+| 水循环回收率 | $\eta$ | 90% | 闭环水系统效率 |
+| 安全缓冲天数 | $T_b$ | 30 天 | 应急储备 |
+
+### 1.2 三种舒适度场景
+
+| 场景 | α值 | 每人日用水 | 含义 |
+|------|-----|-----------|------|
+| **基本生存** | α=1 | 2.9 L | 最低标准 |
+| **舒适生活** | α=50 | 22.5 L | 中等舒适 |
+| **奢侈标准** | α=250 | 102.5 L | 接近地球水平 |
+
+### 1.3 核心公式
+
+**每日循环用水量（水库容量）**：
+$$W_{reservoir} = N \times (w_s + 0.4\alpha)$$
+
+**每日补充量**：
+$$W_{daily} = W_{reservoir} \times (1-\eta) + W_{medical}$$
+
+**首批运输量**：
+$$W_{initial} = W_{reservoir} + W_{daily} \times T_b$$
+
+---
+
+## 二、三场景水需求汇总
+
+| α | 人均日用水 | 水库存量 | 日补充 | 月补充 | 首批运输 | 年补充 |
+|---|-----------|---------|--------|--------|---------|--------|
+| **1** | 2.9 L | 290 吨 | 29.1 吨 | 873 吨 | **1,163 吨** | **10,622 吨** |
+| **50** | 22.5 L | 2,250 吨 | 225.1 吨 | 6,753 吨 | **9,003 吨** | **82,162 吨** |
+| **250** | 102.5 L | 10,250 吨 | 1,025.1 吨 | 30,753 吨 | **41,003 吨** | **374,162 吨** |
+
+---
+
+## 三、四种运输方案
+
+### 3.1 方案定义
+
+| 方案 | 描述 | 日运力 | 比能量 |
+|------|------|--------|--------|
+| **方案1** | 仅空间电梯 (3部) | 1,471 吨/天 | 157.2 GJ/吨 |
+| **方案2** | 仅火箭 (10发射场) | 1,250 吨/天 | 506.1 GJ/吨 |
+| **方案3** | 混合 (电梯+10发射场) | 2,721 吨/天 | 317.5 GJ/吨 |
+| **方案4** | 混合 (电梯+低纬发射场) | 1,971 吨/天 | 243.5 GJ/吨 |
+
+---
+
+## 四、详细结果分析
+
+### 4.1 场景 α=1（基本生存标准）
+
+| 方案 | 首批天数 | 首批能耗 (TJ) | 月补天数 | 月补能耗 (TJ) | 年能耗 (PJ) |
+|------|----------|--------------|----------|--------------|-------------|
+| **方案1: 仅电梯** | **0.79** | **182.8** | 0.593 | **137.2** | **1.67** |
+| 方案2: 仅火箭 | 0.93 | 588.6 | 0.698 | 441.8 | 5.38 |
+| **方案3: 混合全部** | **0.43** | 369.2 | **0.321** | 277.1 | 3.37 |
+| 方案4: 混合低纬 | 0.59 | 283.1 | 0.443 | 212.5 | 2.59 |
+
+**结论**：能耗最优为方案1（仅电梯），时间最优为方案3（混合全部）。
+
+### 4.2 场景 α=50（舒适生活标准）
+
+| 方案 | 首批天数 | 首批能耗 (TJ) | 月补天数 | 月补能耗 (TJ) | 年能耗 (PJ) |
+|------|----------|--------------|----------|--------------|-------------|
+| **方案1: 仅电梯** | 6.12 | **1,415.3** | 4.590 | **1,061.6** | **12.92** |
+| 方案2: 仅火箭 | 7.20 | 4,556.3 | 5.402 | 3,417.6 | 41.58 |
+| **方案3: 混合全部** | **3.31** | 2,858.1 | **2.482** | 2,143.8 | 26.08 |
+| 方案4: 混合低纬 | 4.57 | 2,191.9 | 3.426 | 1,644.1 | 20.00 |
+
+**结论**：水需求增加约8倍，但仍可在1周内完成首批运输。
+
+### 4.3 场景 α=250（奢侈标准）
+
+| 方案 | 首批天数 | 首批能耗 (TJ) | 月补天数 | 月补能耗 (TJ) | 年能耗 (PJ) |
+|------|----------|--------------|----------|--------------|-------------|
+| **方案1: 仅电梯** | 27.87 | **6,445.7** | 20.903 | **4,834.4** | **58.82** |
+| 方案2: 仅火箭 | 32.80 | 20,751.2 | 24.602 | 15,563.8 | 189.36 |
+| **方案3: 混合全部** | **15.07** | 13,016.9 | **11.301** | 9,762.9 | 118.78 |
+| 方案4: 混合低纬 | 20.80 | 9,982.5 | 15.601 | 7,487.0 | 91.09 |
+
+**结论**：高舒适度标准下，水运输成为重要负担，需约1个月完成首批运输。
+
+---
+
+## 五、可视化图表
+
+### 5.1 多场景对比热力图
+
+![multi_scenario_heatmap](multi_scenario_heatmap.png)
+
+### 5.2 首批运输能耗对比
+
+![multi_scenario_initial_energy](multi_scenario_initial_energy.png)
+
+### 5.3 首批运输天数对比
+
+![multi_scenario_initial_days](multi_scenario_initial_days.png)
+
+### 5.4 每月补充能耗对比
+
+![multi_scenario_monthly_energy](multi_scenario_monthly_energy.png)
+
+### 5.5 每月补充天数对比
+
+![multi_scenario_monthly_days](multi_scenario_monthly_days.png)
+
+### 5.6 水需求分解（α=1）
+
+![water_demand_breakdown](water_demand_breakdown.png)
+
+### 5.7 运输方案对比（α=1）
+
+![water_transport_comparison](water_transport_comparison.png)
+
+---
+
+## 六、横向对比汇总
+
+### 6.1 方案1：仅空间电梯（能耗最优）
+
+| α | 首批天数 | 首批能耗(TJ) | 月补天数 | 月补能耗(TJ) | 年能耗(PJ) |
+|---|---------|-------------|---------|-------------|-----------|
+| 1 | 0.79 | 182.8 | 0.593 | 137.2 | 1.67 |
+| 50 | 6.12 | 1,415.3 | 4.590 | 1,061.6 | 12.92 |
+| 250 | 27.87 | 6,445.7 | 20.903 | 4,834.4 | 58.82 |
+
+### 6.2 方案3：混合全部（时间最优）
+
+| α | 首批天数 | 首批能耗(TJ) | 月补天数 | 月补能耗(TJ) | 年能耗(PJ) |
+|---|---------|-------------|---------|-------------|-----------|
+| 1 | 0.43 | 369.2 | 0.321 | 277.1 | 3.37 |
+| 50 | 3.31 | 2,858.1 | 2.482 | 2,143.8 | 26.08 |
+| 250 | 15.07 | 13,016.9 | 11.301 | 9,762.9 | 118.78 |
+
+---
+
+## 七、结论与分析
+
+### 7.1 舒适度系数影响
+
+| 对比 | α=1 → α=50 | α=1 → α=250 |
+|------|------------|-------------|
+| 年补充量增幅 | +674% | **+3,423%** |
+| 占电梯运力 | 1.98% → 15.3% | 1.98% → **69.68%** |
+
+**关键发现**：
+- α=250 时，年水补充量 374,162 吨占电梯运力的 **69.68%**
+- 高舒适度标准将严重影响建设物资运输能力
+
+### 7.2 运输方案选择建议
+
+| 场景 | 推荐方案 | 理由 |
+|------|----------|------|
+| α=1 (基本) | 方案1 (仅电梯) | 能耗低，运力充足 |
+| α=50 (舒适) | 方案4 (混合低纬) | 能耗-时间平衡 |
+| α=250 (奢侈) | 方案3 (混合全部) | 需要最大运力 |
+
+### 7.3 与建设阶段能耗对比
+
+| 场景 | 水供应年能耗 | 建设阶段总能耗 | 占比 |
+|------|-------------|---------------|------|
+| α=1 | 1.67 PJ | 15,720 PJ | **0.01%** |
+| α=50 | 12.92 PJ | 15,720 PJ | **0.08%** |
+| α=250 | 58.82 PJ | 15,720 PJ | **0.37%** |
+
+**结论**：即使在高舒适度标准下，水供应能耗仍不到建设能耗的1%。
+
+### 7.4 Timeline（时间线）总结
+
+| 场景 | 首批运输时间（电梯） | 月补充时间（电梯） |
+|------|---------------------|-------------------|
+| α=1 | **19小时** | 14小时 |
+| α=50 | **6天** | 4.6天 |
+| α=250 | **28天** | 21天 |
+
+---
+
+## 八、文件清单
+
+| 文件 | 描述 |
+|------|------|
+| `water_supply_analysis.py` | 水需求分析核心代码 |
+| `water_transport_results.csv` | 所有场景方案数据 |
+| `water_analysis_report.txt` | 详细分析报告 |
+| `water_demand_breakdown.png` | 水需求分解图 |
+| `water_transport_comparison.png` | 方案对比图（α=1） |
+| `annual_energy_comparison.png` | 年度能耗对比图 |
+| `multi_scenario_heatmap.png` | 多场景热力图 |
+| `multi_scenario_initial_energy.png` | 首批能耗对比 |
+| `multi_scenario_initial_days.png` | 首批天数对比 |
+| `multi_scenario_monthly_energy.png` | 月补能耗对比 |
+| `multi_scenario_monthly_days.png` | 月补天数对比 |
+| `README.md` | 本分析报告 |
+
+---
+
+## 九、如何运行
+
+```bash
+cd /Volumes/Files/code/mm/20260130_b/p3
+python water_supply_analysis.py
+```
+
+可在代码中修改 `comfort_factors = [1, 50, 250]` 分析其他舒适度场景。
+
+---
+
+*报告生成时间：2026-01-31*

p3/water_analysis_report.txt

@@ -0,0 +1,149 @@
+====================================================================================================
+TASK 3: WATER SUPPLY ANALYSIS - MULTI-SCENARIO COMPARISON
+====================================================================================================
+
+## 模型参数
+- 殖民地人口: 100,000 人
+- 生存用水: 2.5 L/人/天
+- 卫生用水: 0.4 × α L/人/天
+- 医疗应急: 5.0 L/次 (每百人每日患病率 2%)
+- 水循环回收率: 90%
+- 安全缓冲: 30 天
+- 舒适度系数场景: α = [1, 50, 250]
+
+====================================================================================================
+## 各场景水需求汇总
+====================================================================================================
+
+α          人均日用水(L)        水库(吨)        日补充(吨)       月补充(吨)       首批(吨)        年补充(吨)      
+----------------------------------------------------------------------------------------------------
+1          2.9             290.0        29.1         873.0        1163.0       10621.5     
+50         22.5            2250.0       225.1        6753.0       9003.0       82161.5     
+250        102.5           10250.0      1025.1       30753.0      41003.0      374161.5    
+
+
+====================================================================================================
+## 场景: 舒适度系数 α = 1
+====================================================================================================
+
+### 水需求参数
+- 每人每日用水: 2.9 L (生存 2.5 + 卫生 0.4)
+- 每日循环水量 (水库): 290.0 吨
+- 每日医疗用水: 0.1 吨 (20 次)
+- 每日循环损耗 (10%): 29.0 吨
+- 每日总补充量: 29.1 吨
+- 每月补充量: 873.0 吨
+- 首批运输量 (水库+30天缓冲): 1163.0 吨
+- 年补充量: 10621.5 吨
+
+### 四种运输方案对比
+
+方案                                       首批天数         首批能耗(TJ)        月补天数         月补能耗(TJ)       
+----------------------------------------------------------------------------------------------------
+方案1: 仅空间电梯                               0.79         182.82          0.593        137.24         
+方案2: 仅火箭（10发射场）                          0.93         588.58          0.698        441.82         
+方案3: 混合（电梯+10发射场）                        0.43         369.21          0.321        277.15         
+方案4: 混合（电梯+低纬发射场）                        0.59         283.14          0.443        212.54         
+
+
+====================================================================================================
+## 场景: 舒适度系数 α = 50
+====================================================================================================
+
+### 水需求参数
+- 每人每日用水: 22.5 L (生存 2.5 + 卫生 20.0)
+- 每日循环水量 (水库): 2250.0 吨
+- 每日医疗用水: 0.1 吨 (20 次)
+- 每日循环损耗 (10%): 225.0 吨
+- 每日总补充量: 225.1 吨
+- 每月补充量: 6753.0 吨
+- 首批运输量 (水库+30天缓冲): 9003.0 吨
+- 年补充量: 82161.5 吨
+
+### 四种运输方案对比
+
+方案                                       首批天数         首批能耗(TJ)        月补天数         月补能耗(TJ)       
+----------------------------------------------------------------------------------------------------
+方案1: 仅空间电梯                               6.12         1415.27         4.590        1061.57        
+方案2: 仅火箭（10发射场）                          7.20         4556.33         5.402        3417.63        
+方案3: 混合（电梯+10发射场）                        3.31         2858.12         2.482        2143.83        
+方案4: 混合（电梯+低纬发射场）                        4.57         2191.85         3.426        1644.07        
+
+
+====================================================================================================
+## 场景: 舒适度系数 α = 250
+====================================================================================================
+
+### 水需求参数
+- 每人每日用水: 102.5 L (生存 2.5 + 卫生 100.0)
+- 每日循环水量 (水库): 10250.0 吨
+- 每日医疗用水: 0.1 吨 (20 次)
+- 每日循环损耗 (10%): 1025.0 吨
+- 每日总补充量: 1025.1 吨
+- 每月补充量: 30753.0 吨
+- 首批运输量 (水库+30天缓冲): 41003.0 吨
+- 年补充量: 374161.5 吨
+
+### 四种运输方案对比
+
+方案                                       首批天数         首批能耗(TJ)        月补天数         月补能耗(TJ)       
+----------------------------------------------------------------------------------------------------
+方案1: 仅空间电梯                               27.87        6445.67         20.903       4834.37        
+方案2: 仅火箭（10发射场）                          32.80        20751.21        24.602       15563.79       
+方案3: 混合（电梯+10发射场）                        15.07        13016.93        11.301       9762.94        
+方案4: 混合（电梯+低纬发射场）                        20.80        9982.48         15.601       7487.04        
+
+
+====================================================================================================
+## 横向对比汇总表
+====================================================================================================
+
+### 方案1: 仅空间电梯
+
+α          首批天数         首批能耗(TJ)        月补天数         月补能耗(TJ)        年能耗(PJ)     
+--------------------------------------------------------------------------------
+1          0.79         182.82          0.593        137.24          1.6697      
+50         6.12         1415.27         4.590        1061.57         12.9158     
+250        27.87        6445.67         20.903       4834.37         58.8182     
+
+### 方案2: 仅火箭（10发射场）
+
+α          首批天数         首批能耗(TJ)        月补天数         月补能耗(TJ)        年能耗(PJ)     
+--------------------------------------------------------------------------------
+1          0.93         588.58          0.698        441.82          5.3754      
+50         7.20         4556.33         5.402        3417.63         41.5811     
+250        32.80        20751.21        24.602       15563.79        189.3595    
+
+### 方案3: 混合（电梯+10发射场）
+
+α          首批天数         首批能耗(TJ)        月补天数         月补能耗(TJ)        年能耗(PJ)     
+--------------------------------------------------------------------------------
+1          0.43         369.21          0.321        277.15          3.3719      
+50         3.31         2858.12         2.482        2143.83         26.0832     
+250        15.07        13016.93        11.301       9762.94         118.7824    
+
+### 方案4: 混合（电梯+低纬发射场）
+
+α          首批天数         首批能耗(TJ)        月补天数         月补能耗(TJ)        年能耗(PJ)     
+--------------------------------------------------------------------------------
+1          0.59         283.14          0.443        212.54          2.5859      
+50         4.57         2191.85         3.426        1644.07         20.0028     
+250        20.80        9982.48         15.601       7487.04         91.0923     
+
+
+====================================================================================================
+## 结论与分析
+====================================================================================================
+
+### 1. 舒适度系数影响
+- α从1增加到250时，年补充量从10621吨增加到374161吨 (+3423%)
+
+### 2. 运输方案比较
+- 能耗最优: 方案1 (仅电梯) - 比能量 157.2 GJ/吨
+- 时间最优: 方案3 (混合全部) - 日运力 2,721 吨/天
+- 平衡方案: 方案4 (混合低纬) - 兼顾能耗和时间
+
+### 3. 占用运力比例
+- α=1: 年补充量占电梯年运力的 1.98%
+- α=50: 年补充量占电梯年运力的 15.30%
+- α=250: 年补充量占电梯年运力的 69.68%

p3/water_transport_results.csv

@@ -0,0 +1,13 @@
+scenario_name,scenario_name_cn,description,daily_capacity_tons,annual_capacity_tons,specific_energy_gj_per_ton,initial_tons,initial_days,initial_energy_gj,initial_energy_pj,monthly_tons,monthly_days,monthly_energy_gj,monthly_energy_pj,annual_tons,annual_days,annual_energy_pj,comfort_factor
+Elevator Only,方案1: 仅空间电梯,仅使用3个太空电梯运输水资源,1471.2328767123288,537000,157.2,1162.9999999999998,0.7904934823091246,182823.59999999998,0.18282359999999998,872.9999999999998,0.5933798882681562,137235.59999999998,0.13723559999999996,10621.499999999998,7.219455307262568,1.6696997999999998,1
+Rocket Only (All Sites),方案2: 仅火箭（10发射场）,使用全部10个发射场进行火箭运输,1250.0,456250,506.0901546034429,1162.9999999999998,0.9303999999999998,588582.849803804,0.588582849803804,872.9999999999998,0.6983999999999998,441816.7049688055,0.4418167049688055,10621.499999999998,8.4972,5.375436577120468,1
+Combined (All Sites),方案3: 混合（电梯+10发射场）,电梯与全部10个发射场火箭混合运输,2721.232876712329,993250,317.4629076645565,1162.9999999999998,0.4273798137427635,369209.3616138792,0.3692093616138792,872.9999999999998,0.32081047067707014,277145.11839115777,0.2771451183911578,10621.499999999998,3.9031940599043535,3.3719322737590867,1
+Combined (Low-Lat Sites),方案4: 混合（电梯+低纬发射场）,电梯与低纬度(<30°)发射场混合运输,1971.2328767123288,719500,243.45730063987577,1162.9999999999998,0.5899861014593466,283140.8406441754,0.28314084064417544,872.9999999999998,0.44287004864489216,212538.2234586115,0.2125382234586115,10621.499999999998,5.388252258512855,2.58588171874644,1
+Elevator Only,方案1: 仅空间电梯,仅使用3个太空电梯运输水资源,1471.2328767123288,537000,157.2,9002.999999999998,6.11935754189944,1415271.5999999999,1.4152715999999999,6752.999999999998,4.590027932960893,1061571.5999999999,1.0615715999999997,82161.49999999997,55.845339851024185,12.915787799999997,50
+Rocket Only (All Sites),方案2: 仅火箭（10发射场）,使用全部10个发射场进行火箭运输,1250.0,456250,506.0901546034429,9002.999999999998,7.202399999999998,4556329.661894795,4.556329661894795,6752.999999999998,5.402399999999998,3417626.814037049,3.417626814037049,82161.49999999997,65.72919999999998,41.58112623745076,50
+Combined (All Sites),方案3: 混合（电梯+10发射场）,电梯与全部10个发射场火箭混合运输,2721.232876712329,993250,317.4629076645565,9002.999999999998,3.308426881449785,2858118.557704002,2.858118557704002,6752.999999999998,2.4815957714573362,2143827.0154587496,2.1438270154587493,82161.49999999997,30.19274855273092,26.08322868808145,50
+Combined (Low-Lat Sites),方案4: 混合（电梯+低纬发射场）,电梯与低纬度(<30°)发射场混合运输,1971.2328767123288,719500,243.45730063987577,9002.999999999998,4.5671924947880465,2191846.077660801,2.191846077660801,6752.999999999998,3.425774843641417,1644067.1512210805,1.6440671512210805,82161.49999999997,41.680260597637236,20.002817006523145,50
+Elevator Only,方案1: 仅空间电梯,仅使用3个太空电梯运输水资源,1471.2328767123288,537000,157.2,41002.999999999985,27.86982309124766,6445671.599999998,6.445671599999998,30752.99999999999,20.90287709497206,4834371.599999998,4.834371599999998,374161.4999999999,254.31833798882673,58.81818779999998,250
+Rocket Only (All Sites),方案2: 仅火箭（10发射场）,使用全部10个发射场进行火箭运输,1250.0,456250,506.0901546034429,41002.999999999985,32.80239999999999,20751214.60920496,20.75121460920496,30752.99999999999,24.602399999999992,15563790.524519674,15.563790524519675,374161.4999999999,299.3291999999999,189.35945138165604,250
+Combined (All Sites),方案3: 混合（电梯+10发射场）,电梯与全部10个发射场火箭混合运输,2721.232876712329,993250,317.4629076645565,41002.999999999985,15.067802668009055,13016931.602969807,13.016931602969805,30752.99999999999,11.301127611376788,9762936.799408104,9.762936799408104,374161.4999999999,137.49705260508426,118.78239772613192,250
+Combined (Low-Lat Sites),方案4: 混合（电梯+低纬发射场）,电梯与低纬度(<30°)发射场混合运输,1971.2328767123288,719500,243.45730063987577,41002.999999999985,20.800687977762326,9982479.698136821,9.982479698136823,30752.99999999999,15.600896455872128,7487042.366578097,7.487042366578097,374161.4999999999,189.81090687977758,91.09234879336685,250