Новым импульсом к развитию метода потенциалов самопроизвольной поляризации (ПС) может стать разработка подходов к численному моделированию и количественной интерпретации каротажных данных, позволяющих повысить информативность и расширить возможности его практического применения. В работе рассматривается постановка задачи ПС с привлечением теории связанных потоков, реалистично отражающей физико-химические процессы формирования ПС для случая диффузии ионов солей в околоскважинном пространстве. Представлен алгоритм численного моделирования данных каротажа ПС на основе метода конечных элементов в рамках осесимметричной модели флюидонасыщенного заглинизированного пласта-коллектора. Для оценки влияния УЭС пород на форму и амплитуду сигналов ПС выполнено моделирование каротажных данных при различных соотношениях УЭС коллекторов и бурового раствора. Возможности алгоритма демонстрируются на практических материалах каротажа ПС с привлечением результатов интерпретации данных ГИС и лабораторных исследований керна из интервалов песчано-глинистых коллекторов нижнего мела Западной Сибири.