A roadmap is outlined for determining comprehensive seismic fragility curves for (single or groups of) spherical pressure vessel structures commonly found in oil refineries. The developed modeling techniques aim to strike a balance between accuracy and computational efficiency, with a focus on capturing the most pertinent failure modes relevant to these structural types. First, a set of “partial” fragilities is determined based on each vessel’s fill ratio, as the response varies substantially depending on the amount of liquid content. Considering that a seismic assessment process invariably involves simultaneous consideration of such partial fragilities, a Monte Carlo-based approach is employed for their combination. The results naturally depend on the level of correlation employed, but can be almost perfectly matched by simpler analytical methods in the edge cases of full and zero correlation.