基于 Go 的免疫治疗门诊离散事件仿真与 ResusBay 挤兑建模

这些让解析解非常难用，但 DES 只需在事件逻辑里表达即可。

3. 模型工程化设计：数据结构与事件机制

3.1 仿真状态（State）要包含什么？

最小可用的仿真状态包括：

• Patients[]：每个患者的到达时间、各站点开始/结束/等待、服务时长、类别、irAE 状态
• Queues[Station]：每站点一个 FIFO 队列
• Resources[Station]：容量 Cap、当前 InUse、BusyMin（用于利用率）
• Events：事件优先队列（最小堆）

3.2 为什么事件队列用最小堆？

事件队列的操作是：

• push：安排未来事件
• pop：取最早事件执行

最小堆让每次 push/pop 都是 O(log N)，非常适合'事件驱动'的仿真。

3.3 BusyMin（利用率）怎么计算才不骗自己？

严格方法：在每次资源占用/释放的瞬间更新占用区间。 为了保持代码'最小可跑'，我们采用常见的简化方式：

• 每当开始服务，BusyMin += dur
• Severe 导致'额外占椅'，用 ReleaseInfusion 事件把额外占用时间补回 Infusion.BusyMin

这让利用率在策略对比上非常稳定且可解释；若要做精确曲线（按分钟利用率热图），可以再加时间片统计（后面会给扩展建议）。

4. 随机性建模：分布选择与参数化

4.1 为什么常用对数正态（LogNormal）？

医疗服务时长经常呈现：

• 下界 > 0（不可能负时间）
• 右偏长尾（少数患者非常久）

对数正态正好符合这个形态。 因此我们用：

• SignIn、Lab、Doctor、Pharmacy、Observation、Resus：LogNormal + clamp 边界
• Infusion：按患者类型给不同均值与方差（可视为截断正态/常数 + 噪声）

4.2 参数从哪来？

实际项目中你会做三层：

1. 经验初值（像本文）用于搭框架
2. 用历史数据拟合（LogNormal/Gamma/Weibull），并做 KS 检验
3. 用'反推校准'：使仿真输出的均值、分位数、资源利用率与真实对齐

技术提示：校准（calibration）比'分布选型'更重要。选错分布但校准对了，模型仍有决策价值；反之则可能误导。

5. 策略建模：预约模板与药房预配怎么'落到代码里'

5.1 预约模板 = 到达时间分布的控制

我们实现两个典型策略：

• Uniform：把到达均匀铺在 6 小时窗口（0~360min），加小扰动
• Staggered：短输注（TypeShort）偏上午，长/联合偏下午（180~360min），加扰动

这本质上是改变系统的'输入过程'，改变拥堵峰值。

5.2 药房预配 = 缩短 Pharmacy 服务时间（同时引入成本问题）

在最小模型中，我们把预配简化为：

• PharmacyPre=true 时，Pharmacy 服务时间视作 0

现实里你会引入：

• 预配提前量（比如提前 30min 开始配）
• 取消/改期概率 → 产生浪费成本
• 稳定性窗口 → 超过窗口失效（需要重配）

这些都可以在 DES 中继续扩展（后面给扩展蓝图）。

6. 免疫治疗独有冲击：irAE 分级与 ResusBay 挤兑

6.1 为什么 Severe irAE 必须单独建模？

轻中度 irAE 通常导致：

• 观察时间延长
• 可能额外医生评估

但 Severe irAE 引入关键耦合：

• 需要 Resus（抢救床位）
• Resus 可能无空位
• 患者在输注区'占位等待转运'
• 输注椅位不能释放 → 后续患者输注排队、等待上升
• 最后表现为尾部变厚（P90、超时率上升）

这个'占位等待'就是医疗系统里非常真实的挤兑机制。

6.2 关键实现点：什么时候释放输注椅？

如果你在 irAE 发生那一刻就释放椅位，会低估真实拥堵。 我们采用更真实逻辑：

• Severe irAE 发生后：仍占椅做初步处置/转运准备（prep 15~30min）
• 准备完入 Resus 队列
• 当 Resus 真正开始服务（即拿到床位）时，再触发 EvReleaseInfusion 释放椅位

这一步是模型'像现实'的关键。

7. 完整可运行 Go 代码（含 ResusBay + irAE）

✅ 复制到 main.go
✅ go run main.go 即可
✅ 输出包含 Avg/P90/超时/各站等待/Resus 等待与利用率/Severe 比例
你可以把参数（到达量、资源容量、irAE 概率）改成你院数据做仿真对比。

package main

import (
	"container/heap"
	"fmt"
	"math"
	"math/rand"
	"sort"
	"strings"
	"time"
)

// ---------- 类型定义 ----------
type PatientType int

const (
	TypeShort PatientType = iota
	TypeLong
	TypeCombo
)

type Station int

const (
	SignIn Station = iota
	Lab
	Doctor
	Pharmacy
	Infusion
	Observation
	Resus // 新增：抢救/留观床位
	Done
)

func (s Station) String() string {
	switch s {
	case SignIn:
		return "SignIn"
	case Lab:
		return "Lab"
	case Doctor:
		return "Doctor"
	case Pharmacy:
		return "Pharmacy"
	case Infusion:
		return "Infusion"
	case Observation:
		return "Observation"
	case Resus:
		return "Resus"
	default:
		return "Done"
	}
}

type IrAEGrade int

const (
	IrAENone IrAEGrade = iota
	IrAEMild
	IrAEModerate
	IrAESevere
)

func (g IrAEGrade) String() string {
	switch g {
	case IrAEMild:
		return "Mild"
	case IrAEModerate:
		return "Moderate"
	case IrAESevere:
		return "Severe"
	default:
		return "None"
	}
}

type Patient struct {
	ID                  int
	PType               PatientType
	Arrive              float64
	Start               map[Station]float64
	End                 map[Station]float64
	Wait                map[Station]float64
	ServiceDur          map[Station]float64 // irAE IrAE IrAEGrade IrAEAt float64 IrAEFired bool // Severe 专用：椅位'占位等待转运'
	InfusionInterrupted bool
	InfusionReleased    bool
	InfusionPlannedEnd  float64      // 原计划输注结束
	InfusionHoldStartTime float64     // 原计划结束时刻（开始'额外占位'计时点）
}

type EventType int

const (
	EvArrival EventType = iota
	EvServiceDone
	EvIrAE
	EvQueueEnter // 在某时刻进入某站点队列
	EvReleaseInfusion // 在某时刻释放输注椅
)

type Event struct {
	T     float64
	Type  EventType
	PID   int
	At    Station
	index int
}

type EventPQ []*Event

func (pq EventPQ) Len() int {
	return len(pq)
}

func (pq EventPQ) Less(i, j int) bool {
	return pq[i].T < pq[j].T
}

func (pq EventPQ) Swap(i, j int) {
	pq[i], pq[j] = pq[j], pq[i]
	pq[i].index, pq[j].index = i, j
}

func (pq *EventPQ) Push(x any) {
	e := x.(*Event)
	e.index = len(*pq)
	*pq = append(*pq, e)
}

func (pq *EventPQ) Pop() any {
	old := *pq
	n := len(old)
	e := old[n-1]
	*pq = old[:n-1]
	return e
}

// ---------- 队列/资源 ----------
type Queue []int

func (q *Queue) Enq(pid int) {
	*q = append(*q, pid)
}

func (q *Queue) Deq() (int, bool) {
	if len(*q) == 0 {
		return 0, false
	}
	pid := (*q)[0]
	*q = (*q)[1:]
	return pid, true
}

type Resource struct {
	Cap     int
	InUse   int
	BusyMin float64
}

// ---------- 仿真 ----------
type Sim struct {
	DayMinutes float64
	Rand       *rand.Rand
	Patients   []*Patient
	Events     EventPQ
	Queues     map[Station]*Queue
	UrgentQueues map[Station]*Queue
	Resources  map[Station]*Reso
}