2025-11-25T19:34:18.315072

Reliability Sensitivity with Response Gradient

Au, Cao

Engineering risk is concerned with the likelihood of failure and the scenarios when it occurs. The sensitivity of failure probability to change in system parameters is relevant to risk-informed decision making. Computing sensitivity is at least one level more difficult than the probability itself, which is already challenged by a large number of input random variables, rare events and implicit nonlinear `black-box' response. Finite difference with Monte Carlo probability estimates is spurious, requiring the number of samples to grow with the reciprocal of step size to suppress estimation variance. Many existing works gain efficiency by exploiting a specific class of input variables, sensitivity parameters, or response in its exact or surrogate form. For general systems, this work presents a theory and associated Monte Carlo strategy for computing sensitivity using response values and gradients with respect to sensitivity parameters. It is shown that the sensitivity at a given response threshold can be expressed via the expectation of response gradient conditional on the threshold. Determining the expectation requires conditioning on the threshold that is a zero-probability event, but it can be resolved by the concept of kernel smoothing. The proposed method offers sensitivity estimates for all response thresholds generated in a single Monte Carlo run. It is investigated in a number of examples featuring sensitivity parameters of different nature. As response gradient becomes increasingly available, it is hoped that this work can provide the basis for embedding sensitivity calculations with reliability in the same Monte Carlo run.

academic

Reliability Sensitivity with Response Gradient

基本信息

论文ID: 2510.09315
标题: Reliability Sensitivity with Response Gradient
作者: Siu-Kui Au (南洋理工大学), Zi-Jun Cao (西南交通大学)
分类: stat.ME cs.LG stat.ML
发表时间: 2025年10月10日
论文链接: https://arxiv.org/abs/2510.09315

摘要

工程风险关注失效的可能性及其发生的场景。失效概率对系统参数变化的敏感性与风险决策相关。计算敏感性比计算概率本身至少困难一个层次，而后者已经面临大量输入随机变量、稀有事件和隐式非线性"黑盒"响应的挑战。基于有限差分和蒙特卡罗概率估计的方法是不可靠的，需要样本数量随步长倒数增长来抑制估计方差。许多现有工作通过利用特定类别的输入变量、敏感性参数或精确或代理形式的响应来提高效率。对于一般系统，本文提出了一种理论和相关的蒙特卡罗策略，使用响应值和相对于敏感性参数的梯度来计算敏感性。研究表明，给定响应阈值处的敏感性可以通过在该阈值条件下响应梯度的期望来表达。确定期望需要在零概率事件的阈值上进行条件化，但可以通过核平滑的概念来解决。所提出的方法在单次蒙特卡罗运行中为所有生成的响应阈值提供敏感性估计。

研究背景与动机

问题背景

工程可靠性分析中，失效概率 $F(y,α) = P(Y ≥ y)$ 对敏感性参数 $α$ 的导数 $∂F/∂α$ 的计算是一个基本但具有挑战性的问题。这里 $Y = f(X,α)$ 是响应函数， $X$ 是随机输入向量， $α$ 是确定性敏感性参数。

核心挑战

高维性: 大量输入随机变量使确定性搜索算法不可持续
稀有事件: 小概率事件使直接蒙特卡罗方法代价高昂
黑盒响应: 隐式非线性响应函数难以近似或建立代理模型
有限差分问题: 传统有限差分方法需要样本数 $N_s = O(Δα^{-1})$ 来抑制方差

现有方法局限性

分布参数方法: 仅限于影响概率密度函数但不影响响应函数的参数
弱方法: 引入增广CDF作为指示函数的平滑近似，但仍需处理表面积分
有限差分: 存在根本性的偏差-方差权衡问题

核心贡献

理论突破: 推导出敏感性的关键公式 $∂F/∂α = ∫ t p(G=t, Y=y) dt$ ，其中 $G = ∂f(X,α)/∂α$ 是响应梯度
条件期望表示: 将敏感性表达为 $∂F/∂α = p(Y=y)E[G|Y=y]$ ，提供直观的概率解释
核平滑方法: 解决零概率事件条件化问题，使蒙特卡罗估计成为可能
子集仿真集成: 将敏感性计算嵌入子集仿真(Subset Simulation)中，在单次运行中获得所有阈值的敏感性估计
通用适用性: 方法适用于任何类型的敏感性参数，不受参数性质限制

方法详解

任务定义

给定响应函数 $Y = f(X,α)$ ，其中：

$X ∈ \mathbb{R}^n$ 是随机输入向量，概率密度函数为 $q(x)$
$α$ 是标量敏感性参数
失效定义为 $Y ≥ y$ ，失效概率为 $F(y,α) = P(Y ≥ y)$

目标：计算敏感性 $∂F(y,α)/∂α$

核心理论

关键公式推导

从第一性原理出发，考虑参数扰动 $Δα$ ，响应扰动为： $Y' = Y + GΔα$ 其中 $G = ∂f(X,α)/∂α$ 是响应梯度。

扰动后的失效概率为： $F' = P(Y ≥ y - GΔα)$

通过条件化消除 $G$ 的随机性： $F' = ∫ P(Y ≥ y - tΔα|G = t) p(G = t) dt$

对每个 $t$ ，进行泰勒展开： $P(Y ≥ y - tΔα|G = t) = P(Y ≥ y|G = t) + tΔα \cdot p(Y = y|G = t) + o(Δα)$

最终得到敏感性公式： $\frac{∂F}{∂α} = ∫ t p(G = t, Y = y) dt$

条件期望形式

利用 $p(G = t, Y = y) = p(G = t|Y = y) p(Y = y)$ ，得到直观形式： $\frac{∂F}{∂α} = p(Y = y) E[G|Y = y]$

蒙特卡罗实现

核平滑技术

由于需要在零概率事件 $\{Y = y\}$ 上条件化，引入核函数 $K(τ)$ 和核宽度 $w$ ： $J(y,α) = ∫ \frac{∂F(τ,α)}{∂α} w^{-1}K\left(\frac{τ-y}{w}\right) dτ = E\left[G w^{-1}K\left(\frac{Y-y}{w}\right)\right]$